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ABSTRACT 


There is an extensive but widely scattered literature on the subject of coronas and reaction rims. 
The observations of many writers are summarized here, and it is shown that further chemical study 
and close observation of geological relations are essential to the solution of the corona problem. 
The writer supplies new chemical analyses of the corona-bearing rock of Gore Mountain, Adiron- 
dacks, and of the garnets in it, and a description and analysis of a recently discovered coronite from 
eastern Quebec. ‘The idea is developed that the fundamental factor in corcna formation may be 
the instability of olivine under thermal metamorphism, leading to the discharge of iron and mag- 
nesium ions which react with plagioclase, the conditions being those of an open system. 


INTRODUCTION 


When chemical reaction takes place between contiguous mineral grains in a rock, 
the new minerals generated often assume the form of sheaths or shells between and 
avuut the reacting grains. These are the ‘‘coronas’”’ of some writers, the “‘reaction 
rims” or “corrosion mantles” of others; and the minerals which compose them were 
called ‘“‘synantetic’” (7.e., meeting) minerals ke’ Sederholm. If the reaction was 
completed in a single step the resulting corona consists of a single shell; if the reac- 


247 














248 S. J. SHAND—CORONAS AND CORONITES 


tion involved two or more steps the corona is formed of two, three, or more shells 
and at least as many different synantetic minerals. It commonly happens that 
cores of the reacting minerals survive unchanged within successive shells of their 
reaction products; in such a case the entire series presents what Sederholm called 
“chemical formulae written in stone.” (The word formulae should be replaced 
by equations since that is clearly what Sederholm had in mind). 

To illustrate this idea, let us recall that crystals of garnet (pyrope) in serpentine 
are often surrounded by reaction rims of kelyphite, a green fibrous material which 
is essentially an amphibole mixed with grains of spinel. Here pyrope and olivine 
or serpentine are the reactants, spinel and amphibole the synantetic products, and 
the following equation seems to be literally “written in stone.” 


(Mg,Fe)sAlkSisOi. + (Mg,Fe)2SiO, = 4(Mg,Fe)SiO; + (Mg,Fe)AlO, 
(pyrope) (olivine) (amphibole)! (spinel) (1) 


Reaction rims may be formed between two solids, as above, or they may be formed 
between a solid’ and a liquid as we see in the zoned plagioclase crystals of many 
volcanic rocks. In the latter case we have to do with a continuous reaction series, 
for the shells consist of the same mineral (feldspar) as the cores; but in the kelyphite 
example, and in all the examples to which the name corona is usually applied, the 
reaction is discontinuous, and the successive shells are formed of different minerals. 
Typical examples of discontinuous reaction are: 

A. Olivine — hypersthene — hornblende plus spinel — plagioclase; 

B. Olivine — hypersthene — garnet — plagioclase. 

A rim formed by reaction is a reaction rim whether it is of the continuous or the 
discontinuous type, so no restriction should be placed upon the use of that term; 
but the special term corona might with advantage be reserved for reaction rims 
of the discontinuous type, especially those composed of two or more shells (two-ply, 
three-ply coronas). It is in this sense that the term corona will be used here, and 
we shall be concerned especially with coronas of types A and B and their modi- 
fications. 


PREVIOUS WORK ON CORONAS 


Attention was first directed to coronas by Térnebohm in 1877, and there is now 
an extensive literature on the subject. Many examples of multiple coronas have 
been found in Norway and Sweden: these have been studied in great detail. Other 
examples have come to light in Karelia, Baden, Bohemia, East Greenland, Tan- 
ganyika, India, Burma, and New South Wales. In North America, coronites have 
been described in Quebec, New Brunswick, New York, and North Carolina, but 
they have not been studied with such care as the Scandinavian occurrences. Memoirs 
of outstanding importance are those of Lacroix (1889), Sederholm (1916), and 
especially the comprehensive study by Brégger (1934-1935). Both Brégger and 
Sederholm give lengthy extracts from the work of earlier observers. 





1In this equation and in succeeding equations the simplest formula for amphibole is used, ignoring the small 
hydroxy! content. 
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CORONAS ABOUT OLIVINE AND PLAGIOCLASE 


The most remarkable examples of corona formation are those which involve as 
reactants the minerals olivine and plagioclase, occurring in such rocks as olivine 
gabbro and troctolite. The coronas are of the two main types indicated above (A 
and B), but there are also coronas of intermediate character; indeed Brégger (Figs. 
15, 23) has found both types within a single thin section. Some of the varieties 
that have been recorded are: 

olivine—hypersthene—actinolite plus spinel—garnet—plagioclase; 

olivine—hypersthene—diopside—garnet—plagioclase; 

olivine—diopside—hornblende plus spinel—plagioclase; 

olivine—anthophyllite—hornblende plus garnet—plagioclase. 
In every case involving olivine the reaction seems to have begun with the con- 
yersion of olivine into pyroxene or amphibole. This material, which constitutes 
the innermost shell of the corona, is practically colorless in thin section, and the 
prisms are so small and so closely packed that specific identification is often difficult. 
Different writers have claimed the mineral to be tremolite, anthophyllite, enstatite, 
hypersthene, diopside, or mixtures of two of these; but Brégger affirms that the 
first shell is mostly bronzite or hypersthene, rarely diopside, and never tremolite. 
Nevertheless Lacroix says that in the gabbro of Le Pallet the olivine sometimes 
passes directly into tremolite or actinolite. 

If there are only two shells in the corona the second shell is generally a symplec- 
tite of pale-greenish amphibole and bright-green spinel. Brégger described the 
amphibole as “‘a greenish-blue actinolite with an extinction angle of 13° to 16° and 
not very strong pleochroism.”” Du Rietz in his account of the Radmansé gabbro 
says 
“In all cases where it was possible to determine the optic character of the axial angle of the amphi- 


bole, it was seen to be negative and the axial angle large... .It is consequently evident that this 
amphibole is chemically lying between the tremolite- actinolite series and the gabbroidic hornblendes.”’ 


A thin shell of diopside has sometimes been observed between the hypersthene 
shell and the actinolite shell, for instance in the coronite of Buck Creek, North 
Carolina; but Pratt and Lewis report that this shell is inconstant, irregular in thick- 
ness, and sometimes almost absent. 

If garnet is present in the corona it normally constitutes the outermost shell, 
in direct contact with plagioclase. This shell may be formed entirely of inter- 
locking garnet grains or of intergrown garnet and amphibole; in the latter case the 
amphibole is generally free from spinel. Brégger showed that a garnet corona has 
often been formed by the transformation of an earlier actinolite-spinel corona. In 
agreement with this view, Holland described coronas of garnet around pyroxene 
in India and noted that the pyroxene is always amphibolized in the neighborhood 
of the garnet. Holland argued, therefore, that amphibolization seems to be “ 
preliminary change in the augite, necessary for the formation of the garnet.”” Never- 
theless Eskola found that in the coronas at Aurlandsfjord, Norway, “hornblende 
is clearly posterior to the garnet.” Chemical analyses of corona garnets from 
Norway and from the Adirondacks are given in a later section. If we may judge 
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by the number of described examples, the amphibole-spinel type of corona is much 
commoner than the garnet type. 


CHEMICAL EQUATIONS 


It has been customary to represent the reactions that produce coronas by means 
of balanced equations of the type A + B = C + D, employing exact numbers of 
entire molecules of each reactant. Thus to explain the production of garnet from 
anorthite and olivine we have the well-known “Becke equation’’: 


(Mg,Fe)2SiO, _ CaAl,SizOs = Ca(Mg,Fe)sAlSi;0;2 
(olivine) (anorthite) (garnet) (2) 


To illustrate the formation of clino-pyroxene one may write 


2(Mg,Fe)2SiO, + CaAlSi,O; = Ca(Mg,Fe)Siz0. + 2(Mg,Fe)SiO; + (Mg,Fe) Al,O, 
(olivine) (anorthite) (diopside) (hypersthene) (spinel) (3) 


For the generation of actinolite plus spinel: 


2(Mg,Fe)2SiO, + CaAlSig0s = Ca(Mg,Fe);Si,Oi2. + (Mg,Fe) Al.O, 


(olivine) (anorthite) (actinolite) (spinel) (4) 


These nicely balanced equations have a certain interest as exercises in chemical 
book-keeping, but they give an entirely false picture of what actually happened. 
The reactions that produced the corona minerals took place within an essentially 
solid medium consisting of crystals of olivine and plagioclase separated by films 
of water or water vapor, in which both olivine and anorthite are insoluble. Under 
these circumstances there can be no double decomposition involving exact numbers of 
molecules. What actually took place was a reaction between wandering ions of 
Mg, Fe, Ca, Al, and SiO, which were discharged at different rates from the olivine 
and the plagioclase crystals and travelled with different velocities through the water 
film. It is impossible to draw up a balance sheet between wandering ions in an 
open system. 

To illustrate the formation of garnet, for example, we should replace equation 
(2) by something like this: 


Mg*+ + Fet*+ + Ca++ + Al*+*+ + Fet++ + Si0z--- — (Ca,Mg,Fe)s(Al,Fe)2(SiO,)s 


(ions) (garnet) (5) 


The relative masses of olivine and plagioclase do not enter into the equation at all; 
it is only necessary that there should be enough of each to keep up the supply of 
ions, and this depends on the surface areas of the two minerals, not on their pro- 
portions by weight. 

The fallaciousness of the Becke equation (2) is shown by the composition of the 
supposed product, which is a garnet containing 33 molecules per cent of grossularite. 
No garnet of this composition has yet been found, and it seems to be impossible in 
nature, as Boeke showed by a statistical study in 1914. 
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In an attempt to show what becomes of the albite contained in the reacting 
plagioclase, Brégger wrote some very complex equations, for example 


18MgO, 8FeO, 13Si0. + Na,O, ALO, 6SiO2 + 4CaO, 4AL0;, 8SiO- 


(13 olivine) (albite) (4 anorthite) 
= 9MgO, 4FeO, 13Si02 + 6MgO, 3FeO, 4CaO, Na,O, Al,Oz, 14510. + 3MgO, FeO, 4Ai.0; 
(13 hypersthene) (2 hornblende) (4 spinel) (6) 


Because one can write an equation which balances arithmetically it does not by any 
means follow that the reaction is a possible one or that it happened in just that 
way. The hornblende in coronas, like the garnet, was formed by wandering ions, 
and its composition cannot be deduced from the relative masses of the primary 
minerals in the rock. 

It is customary to represent the conversion of olivine into hypersthene as an 
addition of silica, but there is no evidence that silica was added in the formation of 
coronas and there is good reason to doubt it. One of the most characteristic of 
corona minerals, spinel, is never found in quartz-bearing rocks; it appears to be 
unstable in presence of free silica. It is more likely that the hypersthene shell in the 
two-ply and three-ply coronas of gabbroic rocks was formed by the expulsion of iron 
ions and magnesium ions from olivine, thus 


2(Mg,Fe)2SiO, + 2H,O — 2(Mg,Fe)SiIO; + Mgt+ + Fet+* + 4(OH)- 
(olivine) (hypersthene) (ions) (7) 


That this reaction is not an imaginary one is shown by the common appearance of 
swarms of magnetite during the serpentinization of olivine. 

But conventional equations have their uses. If we combine equations (2) and 
(4) we see that 


Ca(Mg,Fe)3SisOie a (Mg,Fe) Al.O, = Ca(Mg,Fe)2Al-Si;Oi2 (Mg,Fe).SiO, 


(actinolite) (spinel) (garnet) (olivine) (8) 


That is, an actinolite-spinel corona is chemically equivalent to a garnet corona (dis- 
regarding for the moment the precise composition of the garnet). The kelyphite 
equation (1) illustrates the same point and shows that the reaction is reversible. 


CHARACTERS OF CORONA-BEARING ROCKS (CORONITES) 


The rocks that give rise to these coronas between olivine and plagioclase are all 
of one mineralogical type which different writers have called olivine gabbro, olivine 
norite, olivine hyperite, troctolite, or allivalite. Brégger says the Norwegian 
hyperites are medium-grained rocks with ophitic texture, normally holding about 14 
per cent of olivine but sometimes as much as 25 per cent. The plagioclase ranges 
in composition from Ab,;An; to Ab;Ans and is “very seldom so basic as Ab,Ane.” 
Du Rietz maintains that at Radmansé “‘it is the rocks rich in (normative) anorthite 
and olivine and poor in normative diallage and hypersthene, which have the best 
developed reaction rims.’ Most observers are in agreement about the relative 
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abundance of olivine in coronites, but there is some difference of opinion about the 
composition of the plagioclase. Du Rietz found the plagioclase at Radmansé to be 
Ango; Lundgardh found Angg at Grévstanis; and Teale said that in his Tanganyika 
example “the extinction angles observed in the feldspar indicate anorthite.” On 
the other hand Kranck states that the plagioclase in a coronite from Liverpool 
Land, East Greenland, is about Ango, and Buddington has recorded Ang, Ango, 
and An;; in Adirondack coronites. 

It is clear that the composition of the surviving plagioclase in any coronite must 
depend on the initial composition of the mineral and also on the proportion of 
plagioclase to olivine in the original rock. A very feldspathic troctolite may de- 
velop coronas without any great change in the anorthite content of the feldspar, 
whereas a feldspar-poor example must give up such a large proportion of its anor- 
thite to form coronas that the residual feldspar is likely to be very sodic. But the 
normative plagioclase will generally be rich in anorthite, as Du Rietz claims. 

In short, the kind of rock that most commonly develops coronas seems to range 
mineralogically from olivine norite to troctolite and anorthite troctolite (allivalite, 
harrisite). But the great majority even of troctolites are not coronitic, so it is 
certain that mineralogical composition is not the sole determining factor in corona 


formation. 
CHEMICAL COMPOSITION OF CORONITES 


Very few coronites seem to have been chemically analyzed. Brégger furnished 11 
analyses of olivine hyperites from the Krageré region, where some of the best ex- 
amples of corona formation are found, but there is a tantalizing uncertainty whether 
the analyzed specimens were corona-bearing or not. On page 3 of his monograph 


Brégger wrote 


‘St was necessary to procure a number of new, modern analyses of wholly fresh types of olivine- 
hyperites, from occurrences where the olivine of the rock is still preserved and not already replaced 
by rhombic pyroxene”; and (footnote 4, p. 3) “The hyperites from Odegirden, Bamle, by which 
the olivine is already mostly replaced by rhombic pyroxene, represent not wholly unaltered, fresh 
rocks, and the analyses of the hyperites from this occurrence are thus not included in the table.” 


These statements seem to mean that the hyperite specimens chosen for analysis 
were free from coronas, so one is left in doubt as to the value of the 11 analyses as 
evidence of the composition of coronites. The average of Brégger’s 11 analyses, 
and the corresponding norm, are given in Table 1. 

The norm falls in the subrang auvergnose, III.5.4.4-5. There is nothing remark- 
able about either the analysis or the norm. The figures are closely comparable with 
data for average olivine norite, olivine gabbro, and olivine diabase. The only 
point that might give rise to comment is the high proportion of FeO in relation to 
Fe,O; and the consequent richness of the dark silicates in iron. But a large excess 
of FeO over Fe,O; is a common feature in this subrang. 

There are many coronites among the gabbros and norites of the Adirondacks, 
and many chemical analyses of these rocks are now available, but only in a few 
instances do we know that the rock analyzed was a coronite. Among the older 
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descriptions there is one of the gabbro of Split Rock Mine, which was analyzed by 
W. F. Hillebrand and briefly described (with a colored plate showing garnet coronas) 
by Kemp (1899, p. 383). This analysis is given in Table 2, column 1. Three recent 


TABLE 1.—Average composition of 11 olivine hyperites from Krageré 
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TABLE 2.—A nalyses of coronites from the Adirondacks 




















at) | (2) | 3) a) 
rs 47.88 | 46.49 | 47.99 | 48.30 
RIES 18.90 | 18.93 | 16.90 18.45 
ere 1.39 | 1.08 1.49 | 1.73 
(Fe, Mn)O........ 10.61 9.55 | 9.83 | 9.97 
ae 7.10 10.04 | 9.56 | 5.49 
eee 8.36 8.91 | 8.99 | 8.30 
I his sg 2.75 | 2.61 | 2.85 | 3.37 
ee 81 ~~ 4 . a 75 
. oaaeeee 1.20 1.28 .93 | 2.22 
ee 1.02 | .90 | 45 | 1.56 
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(3) Buddington no. 59. 
(4) Buddington no. 66 





(1) Hillebrand 
(2) Buddington no. 56. 


analyses of coronites from the northwestern Adirondacks, accompanied by short 
petrographic descriptions, were given by Buddington (1939). They are included 
in Table 2, columns 2, 3, 4. 

Analysis no. 1 falls in hessose, II.5.4.4. No. 4 also falls in hessose, near auverg- 
nose; and no. 2 is just on the boundary line between these two subrangs. These 
three are garnet coronites. No. 3 is a hornblende-spinel coronite without garnet; 
it falls in auvergnose, III.5.4.4. 

The close similarity of these four analyses, and their remarkable resemblance 
to Brégger’s average for the Krageré district, may give the impression that corona 
formation is restricted to rocks within a very narrow range of chemical composition; 
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but before we draw any conclusions let us look at another group of coronite analyses, 
those from northern Karelia, on the shores of the White Sea. The garnet coronites 
of this region were described by Machkovtsev who supplied the four analyses quoted 


in Table 3. 


TABLE 3.—A nalyses of coronites from Karelia 











| (1) | (2) (3) | (4) 
Ses | 45.59 49.60 50.42 | 52.43 
eh a 7.74 12 75 10.67 | 9.36 
ee 2.06 | 47 | 1.39 | 1.98 
(Fe, Mn)O........ 11.13 9.73 9.23 12.36 
“ESSE 24.34 | 16.80 | 17.36 | 6.46 
ae 6.69 9.14 | 8.49 | 10.87 
Na:O 1.23 | 1.14 | 1.49 | 2.99 
KO... .58 67 66 | 1.11 
Tae .30 | 22 .30 1.71 
Se .12 | 13 .14 .20 

100.54 | 100.65 100.13 99.47 





Of these analyses the first falls in custerose, IV.1.4.1.2; the second in auvergnose, 
III.5.4.4.; the third just misses auvergnose and falls in palisadose, IV,1.2.2.2.; 
and the fourth, on account of its greater content of alkalis and silica, falls in kilauose, 
III.5.2.4. These rocks are much less feldspathic than the Adirondack examples; 
the normative feldspar amounts to 27, 41, 37, and 41 per cent, and the normative 
olivine to 46, 21, 18, and 2 per cent. The fourth analysis is remarkable by yielding 
a plagioclase of the (normative) composition AbsAn; together with only 2 per cent 
of normative olivine. It seems possible that this rock may have received some 
addition of soda and silica from the aplite veins which traverse the area. 

If we now combine Tables 2 and 3 we see that the range of composition has been 
considerably increased. If in addition we take account of Brégger’s 11 analyses 
(not quoted separately here) then we arrive at the following range of composition 
in 19 known coronites: 


a abate Oe a, 45.6-52.4 
caida shay sats ie wowace 7.7-21.4 
oe a Nach ks 6.9-12.3 
DO rene Pi ee Pe ot 4.0-24.4 
SSC anne oe a re 6.7-11.1 
a err 9- 4.8 


If we add the coronite from eastern Quebec which is described in a later section 
of this paper, the range of composition is still further widened, the maximum for 
FeO being increased to 17.6 and the minimum for CaO being lowered to 4.4 per cent. 

Thus it would be wrong to conclude that coronas are restricted to rocks ofa 
narrow range of composition, and equally wrong to suppose that every rock that 
falls within this range is necessarily a coronite. Washington tabulated 197 analyses 
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in auvergnose and 165 in hessose, but only 2 of the corresponding rocks are known 
to be coronites. Clearly, then, the precise chemical composition of a rock is not 
the main factor in corona formation. 

A few other analyses of coronitic rocks are on record, but their value is doubtful. 
Du Rietz quoted three analyses of coronites from Radmansé, made in the year 1872. 
They are remarkably good analyses for their day, but there is reason to suspect 
that alumina was overestimated, so we shall not make use of them. Two of the 
norms fall in corsase, IIT.5.5, and the third in kedabekase, ITI.5.5. 

Lacroix gave in 1899 an analysis of a coronitic facies of the gabbro of Le Pallet. 
It is incomplete and was judged fourth-rate by Washington. The norm falls in 
hessose, IT.5.4.4. 

Pratt and Lewis furnished an analysis of the coronitic troctolite of Buck Creek, 
North Carolina. Unfortunately it is impossible to reconcile the analytical data with 
the mineralogy of the rock. Washington rates this analysis C4 or fifth-rate. 


POSSIBLE IMPORTANCE OF FeO 


That an abundance of ferrous oxide may be an important factor in corona forma- 
tion is indicated by the common occurrence of coronas about iron ore. Térnebohm 
in his pioneer study of coronas recorded that the ilmenite in Norwegian hyperites 
is often enclosed in a shell of deep red-brown mica, and that in turn by a shell of 
green hornblende. Brégger added that the hornblende “is often symplectically 
filled with spinel.” Brégger noted that shells of garnet are sometimes formed 
directly between iron ore and plagioclase, but he said that the general order of 
succession of synantetic minerals about iron ore is (1) biotite; (2) brownish horn- 
blende with spinel; (3) dark-green hornblende with garnet but without spinel. 

Thus in the same Norwegian rocks both olivine and iron ore have reacted with 
plagioclase to give hornblende, spinel, and garnet. Between olivine and iron ore 
there is nothing in common except an abundance of ferrous oxide. But even an 
abundance of ferrous oxide in the heavy silicates does not by itself cause the forma- 
tion of coronas. The ferro-gabbros of the Skaergaard intrusion in East Greenland 
contain abundant olivine ranging in composition from Faygs to Fayos; yet Wager 
and Deer described no coronas about the olivine in these rocks and only narrow fibrous 
rims about the iron ore. 


PRESENT STATUS OF THE PROBLEM 


It has been demonstrated that neither a lime-rich plagioclase nor an abundance 
of olivine nor a high content of ferrous oxide in the silicate minerals is sufficient 
to determine the formation of coronas, although all three are favorable factors. 
It seems to follow that the factor which “pulled the trigger” that set the reaction 
going was a physical one. Before we can reach any safe conclusion about the cause 
of corona formation we must know more about the geological environment and 
geological history of coronites. Information of this character is deplorably scarce. 
Lacroix noted that the corona-bearing rock of Odegirden is intersected by dikes of 
“granulite” or granite pegmatite. Du Rietz says that the ancient rocks into which 
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the gabbro of Radmansé was intruded were afterwards invaded by granite. Lund- 
gardh found that the Grovstanis region is traversed by many dikes of aplite, peg. 
matite, and amphibolite; one of the pegmatite dikes indicated on Lundgirdh’s map 
is nearly 100 meters wide. Brégger had little to say about field relations, but he 
noted that the hyperites of Gomé and Langé are cut by dikes of porphyritic hyperite, 

The coronite of Gore Mountain, in the Adirondacks, is cut by syenite and granite 
(Krieger). Philipp reported that the coronitic gabbro of Wiesental, Baden, is cut 
by hornblende granite, and detached blocks of the gabbro are enclosed in granite. 
A similar observation was made by Kranck, who found boulders of coronitic diabase 
enclosed in ‘‘migmatitic gneiss” near Cape Greg, East Greenland. Lacroix found 
that the gabbro of Le Pallet is cut by many dikes of lamprophyre and fine-grained 
granite. In northern Karelia, Machkovtsev described a series of gabbroic rocks 
which have been invaded by great numbers of aplitic intrusions. As a result of 
contact metamorphism the gabbros and norites were converted first into coronites 
and further into amphibolites. 

Of the writers named above, Philipp, Kranck, and Machkovtsev have expressed 
themselves unequivocally in favor of thermal metamorphism as the cause of corona 
formation. Buddington favors the wider term “dynamo-thermal metamorphism.” 
Sederholm said he regarded typical coronas as “truly secondary and metamorphic” 
and added that “‘some of the phenomena in question . . . occur in rocks which have 
undergone contact action of the most unquestionable character.” But the older 
view, that coronas are due to late-magmatic reactions or to the operation of “min- 
eralizing agents,” still has its supporters. 

It should now be evident that the problem of corona formation demands for its 
solution a much more careful study than has yet been made of the geological con- 
ditions under which coronas are formed, as well as much additional information 
about the chemistry of coronas and synantetic minerals. As a contribution toward 
the eventual solution of the problem the writer offers a chemical study of the garnet 
coronite of Gore Mountain, and of the garnets in it, and a description and analysis 
of a newly discovered actinolite-spinel coronite from eastern Quebec. 


GARNET CORONITE OF GORE MOUNTAIN, ADIRONDACKS 


Coronas in Adirondack gabbros have been described by Kemp (1894), Miller 
(1914), Gillson, Callahan, and Millar (1928), and Krieger (1937), but no full descrip- 
tion worthy to compare with Brégger’s great work on the Norwegian coronites has 
yet appeared. The description of the Gore Mountain coronite by Krieger is very 
brief, but the photographs that accompany it show that the coronas consist of (1) 
a shell of pyroxene, described as “augite, hypersthene, and clino-hypersthene”; 
(2) a narrow shell of plagioclase and biotite; (3) a shell of intergrown garnet, horn- 
blende, and pyroxene; (4) a shell of pure garnet in contact with plagioclase. 

A chemical and physical study of Adirondack garnets and the rocks that contain 
them has been made by Mr. S. B. Levin in the laboratories of Columbia University 
and Hunter College. To that study (not yet published) the present writer con- 
tributes a chemical analysis of the Gore Mountain coronite (Table 4) and an analysis 
of garnets from the coronas in this rock (Table 5). 
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This analysis compares closely with some of those given by Brégger and especially 


with the average which we have already discussed (Table 1). 
subrang (auvergnose) as Brégger’s average. 


TABLE 4.—Chemical analysis of the Gore Mountain coronite 


(S. J. Shand, analyst) 


It falls in the same 
Out of 197 analyses included in this 
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TABLE 5.—Chemical analyses of corona garnets 
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subrang, in Washington’s Tables, only about 1 in 10 shows more than 20 per cent 
of normative olivine, so in this respect the Gore Mountain rock is above the average. 
The large excess of ferrous oxide over ferric oxide is noteworthy; it is reflected in the 
high iron content of the normative olivine and hypersthene. 

The minute garnets in this rock were separated by repeated fractionation with 
Clerici solution. A perfect separation could not be effected, but it is thought that 
the final product had a purity of perhaps 95 per cent. An analysis of corona garnets 
from Aurlandsfjord, Norway, by Eskola is added for comparison. 
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If we combine the figures for Al,O3, Fe2O3, Mn2O3, and TiO: (or TigO3) and write 
the sum as R.O;, then the composition of these garnets works out as follows: 
Molecules per cent 


Gore Mtn. Aurlandsfjord 
SED Ds as MOMMNRINGINE) 55 oho. 5a Sis ovine s ole Bibs Ae BSE SG Ses 41 38 
0 inns gs das eblin adliownctreees 35 44 
RGR FURAN Dg SU ROTMIERTERE) oo) Sook so soe booeG ants esas sans 19 18 


Here we have a clear demonstration of the failure of Becke’s equation (2). To 
derive a garnet of the above composition (approximately 4 molecules of iron-mag. 
nesium garnet to 1 molecule of calcium garnet) from olivine and anorthite it would be 
necessary to write 


5CaAl.Siz0; + 5(Mg,Fe)2SiO; + 4(Mg,Fe)SiO; 
(S anorthite) (S olivine) (4 hypersthene) 
= 5[{Caoe(Mg,Fe)e.JAleSisO2 + 2Ca(Mg,Fe)Si.0, 
(5 garnet) (2 diopside) (9) 


This equation gives a garnet of the desired composition, but in other respects it 
is untenable. It is certain that all equations of this type, which involve the bal- 
ancing of whole numbers of molecules, are equally fallacious under the conditions 


that give rise to coronas. 


NEW CORONITE FROM EASTERN QUEBEC 


The rock now to be described was collected by Mr. David St. Clair in the course 
of an investigation carried out for the Freeport Sulphur Company. The locality 
is on the east side of Lake Napetipi, in Pontchartrain Township, Saguenay County, 
Quebec, about 45 miles from the Labrador boundary. Latitude and longitude are 
approximately 51° 25’ and 58° 5’. Adams described a coronite from Saguenay in 
1893, but the new locality is more than 600 miles northeast of the former one. Mr. 
St. Clair reports that 


“the basic intrusive body is ovate in shape, with a long axis of approximately two-thirds of one 
mile. The structure is distinctly that of a basin, the sheeting of the old gneiss series dipping be- 
neath the intrusion on all sides at angles varying between 10° and 45°.” 

Samples collected along two lines of traverse were submitted to the writer for 
microscopic study. Of the 14 specimens received, one proved to be a very beautiful 
hornblende-spinel coronite. In view of the rarity of such rocks this example was 
thought worthy of detailed examination, and the Freeport Sulphur Company has 
kindly consented to the publication of the description that follows. The writer 
expresses his appreciation of this courtesy and his thanks for the geological data so 
willingly imparted by Mr. St. Clair. 

The floor of the intrusion is a banded quartz-plagioclase-garnet gneiss with inter- 
calated sillimanite gneiss. A specimen of the eruptive rock collected 3 feet above 
the contact is a fine-grained biotite norite with aplitic texture. Eight feet above 
the contact the rock is a coarse-grained, melanocratic biotite norite with light and 
dark bands, and it retains this character up to at least 28 feet above the contact. 
Near the top of the hillside, some 250 feet away from the exposed contact and perhaps 
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150-175 feet perpendicularly above the floor of the intrusion, a single specimen was 
collected which proved to be the corona-bearing troctolite now to be described. 
Specimens collected at random over the central part of the intrusion include a coarse- 
grained hornblende rock, a hornblende diorite, and two samples of hornblende-bearing 


harzburgite. 


TABLE 6.—A nalysis of coronite from eastern Quebec 
(S. J. Shand, analyst) 
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Sp. gr. 3.36 at 4° C. 


It is impossible without further geological data to interpret the internal structure 
of the intrusion. The troctolite is so rich in olivine as to suggest an origin by crystal 
settling, yet both troctolite and harzburgite, the densest rocks in the series, were 
found about the top of the exposed mass, overlying lighter norite. It is unlikely, 
then, that the intrusion is a sill. If it is funnel-shaped, as Mr. St. Clair suggests, 
then the olivine rocks may constitute the core and the norite the marginal facies. 
Further speculation is futile, but at least there seems no reason to doubt that all 
the rocks are related and that they were derived from a single intrusion of a highly 
ferromagnesian magma. Of special interest for our investigation of the genesis of 
coronas are the following observations: (1) The rock which developed coronas is 
the only one of the 14 samples which contains both olivine and plagioclase; and (2) 
outcrops of granite pegmatitite and aplite were observed by Mr. St. Clair among 
the ferromagnesian rocks. 


CHEMICAL COMPOSITION OF THE NEW CORONITE 


The only specimen of this rock was fortunately big enough to yield ample material 
for the analysis given in Table 6. 

This analysis falls in the subrang argeinose, IV.1.5.1.2. If we assume that there 
was little gain or loss of material in the conversion of the rock into coronite, then the 
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original rock must have been a troctolite containing about 52 per cent of olivine 
and probably less than 32 per cent of plagioclase, with some clino-pyroxene and a 
little iron ore. This is an uncommon type of rock. In Washington’s Tables there 
are only three analyses in the subrang argeinose. One of them is the original “har. 
risite” (anorthite-olivine rock with predominant olivine) from Rum, Scotland; 
another is a hornblende peridotite from the Pyrenees; and the third is a picrite from 
Nassau. All three rocks show even more (normative) olivine than the Labrador 
coronite, yet none of them is corona-bearing. It is perhaps significant that in each 
of these rocks the ratio of FeO: MgO is less than half as great as in the Quebec rock, 


MINERALOGY OF THE NEW CORONITE 


Minerals recognizable in the hand specimen are olivine, in clear green grains which 
weather red brown; biotite, in conspicuous brown flakes; glittering black grains of 
spinel; and irregular tables of greyish-green feldspar up to 5 millimeters wide. Under 
the microscope nine minerals are distinguishable which fall into three groups as 
follows: 
Pre-corona minerals: olivine, clino-pyroxene, plagioclase, ilmenite; 
Corona minerals: hypersthene, spinel, cummingtonite; 
Post-corona minerals: biotite, apatite. 
The coronas around the olivine grains are built up as follows: 

Ist shell: hypersthene (with occasional grains of cummingtonite); 

2nd shell: symplectite of hypersthene and spinel; 

3rd shell: cammingtonite with or without spinel. 
No garnet has been observed. Some clino-pyroxene may possibly be associated 
with the hypersthene, but this has not been established. Either the second shell 
or the third may be in contact with plagioclase. 


OLIVINE forms roughly half of the rock. The grains, which reach 2-3 mm. in diameter, are 
invariably rounded by corrosion and enveloped in reaction products, chiefly hypersthene but some- 
times cummingtonite (Pl. 1). Olivine grains are also enclosed in plates of clino-pyroxene of pre- 
corona age. The olivine is clear and colorless under the microscope, showing only a slight discolora- 
tion along cracks. On the universal stage several grains showed an optic axial angle of 2V(—) = 

°-86°, indicating a content of about 24 molecules per cent of fayalite or 20 per cent of FeO by 
weight. 

HyPERSTHENE forms shells around most of the olivine grains. The shells, which are commonly 
from 0.2 to 0.5 mm. thick, are formed of short overlapping rods which stand roughly at right angles 
to the surface of the olivine but do not all have the same crystallographic orientation. A single 
cleavage is often developed, but in these small grains it is hardly ever possible to see two cleavage 
directions. That the mineral is hypersthene is shown first of all by the pleochroism. In the thicker 
sections the characteristic pleochroism from rose-pink to colorless or very faint green is convincingly 
shown. For further identification a large number of grains were examined on the universal stage, 
but on account of their small size and wavy extinction few of them gave useful measurements. 
The optic axial angle, in a number of grains which showed a single optic axis, appeared to be 2V(—) = 
64°-68°, but a grain showing both optic axes gave 2V(—) = 72°. These measurements, though 
far from satisfactory, indicate a probable content of 30 molecules per cent of FeSiQs, or about 20 
per cent of FeO by weight. 

SPINEL (PLEONASTE): For a distance of a few tenths of a millimeter from the margin of a2 olivine 
grain the hypersthene is free from enclosures; then it passes abruptly into a micrographic inter- 
growth with spinel. This vivid green mineral forms angular grains, droplike particles, and thin, 
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wormlike threads which trend roughly at right angles to the wall of the hypersthene shell (PI. 1, 
fig. 2). The colorless mineral between the threads does not differ in any recognizable respect from 
the adjacent hypersthene; even the pinkish pleochroism can still be recognized here and there. 
The pleonaste threads are rounded in cross section; they are about .01 to .02 mm. thick and 10 to 
20 times as long. They are neither quite straight nor quite parallel; sometimes they radiate like the 
ribs of a fan, and occasionally they throw off branches. By rough estimation the spinel may form 
as much as half of the aggregate. This hypersthene-spinel symplectite is identical in appearance 
and texture with the actinolite-spinel symplectite which was so carefully described and illustrated 
by Brégger, but it contains no actinolite or any other amphibole. A pyroxene-spinel symplectite 
has been recorded before (by Lacroix in the Pyrenees; by Benson in New South Wales; by Barth 
in the Seiland district), but it seems to be much rarer than the actinolite-spinel type. Benson and 
Barth mention both ortho- and clino-pyroxene in the symplectite but give no exact data regarding 
them. In the present instance only hypersthene was identified. 

Spinel occurs also in the brown amphibole (cummingtonite), but only in scattered grains; the 
vermicular habit has not been observed in this case. In the plagioclase too spinel dust and larger 
grains are fairly common. The grains tend to be arranged in parallel planes following the twinning 
lamellae of the plagioclase. It is remarkable that, where hornblende and plagioclase meet, the 
lines of spinel grains often continue without change of direction from the plagioclase into the amphi- 
bole, as if the latter had replaced the plagioclase without disturbing the enclosed spinel grains (PI. 1, 
fig. 2). It is interesting to read an exactly similar observation by Lacroix (1899, p. 233) in his 
account of the gabbro of Odegirden: “II est a remarquer que ces inclusions disposées suivant l’axe 
vertical du feldspath se continuent dans l’amphibole en restant orientées comme dans le feldspath 
voisin.”” 

CiINo-PyROXENE: This mineral forms a few relatively large poikilitic plates and skeleton crystals 
which enclose many olivine grains. It was probably ophitic toward plagioclase, but hornblendiza- 
tion has obscured this relation. The mineral is not simple; it has, first, a minute lamellar structure 
due to the parallel intergrowth of two kinds of pyroxene which do not extinguish in the same position; 
in addition it is intergrown on a coarser scale with brown hornblende (cummingtonite) formed by 
its own partial transformation. The evidence of transformation is complete. The normal clino- 
pyroxene is colorless, but in some parts it develops a slight brownish tinge, accompanied by weak 
pleochroism. In a further stage one can see that the brown coloration is due to a multitude of tiny 
platelets of hornblende. Larger areas, irregularly placed, become completely trar formed into 
brown hornblende which sometimes forms graphic intergrowths with the pyroxene, and the margins 
of the crystals are generally hornblendized along most of their length. Grains of olivine enclosed 
in the pyroxene are angular to rounded, as if they had been partly resorbed during the growth of 
the pyroxene. When an olivine grain lies partly inside and partly outside of a pyroxene plate, 
a corona is developed on the outside only. This relation seems to mean that the crystallization of 
clino-pyroxene preceded corona formation. Universal-stage measurements gave for the pyroxene 
2V(+) = 60°-64°, with an extinction angle exceeding 30°. Clino-pyroxene was not identified as 
a constituent of the coronas. 

CUMMINGTONITE: The formation of a brown amphibole by transformation of clino-pyroxene has 
just been described. In addition, a pale-green amphibole appears in the third shell of the coronas. 
The green and the brown varieties are not sharply distinguishable, many grains having intermediate 
tints which might be described as greenish-brown or brownish-green. The deepest brown is seen 
where the amphibole is in contact with iron ore; the purest green appears in amphibole which borders 
plagioclase; but there is a complete gradation from green to brown. The pleochroism is: Z, deep 
brown or greenish brown; Y, yellowish brown or pale green; X, colorless. On the universal stage 
every crystal, whether brown or green, proved to be optically positive with 2V(+) = 84°-86° and 
a maximum extinction angle of 20°. Clearly this mineral is neither actinolite nor common horn- 
blende; it is cummingtonite, and it must contain something like 65 molecules per cent of griinerite, 
H,Fe7SisOx, or (very roughly) 25 per cent of FeO by weight. 

Grains of spinel may be present either in the brown or in the green amphibole, but much of the 
amphibole is free from spinel. The vermicular type of hornblende-spinel symplectite has not 
been observed in this rock. 
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PracriociasE: Every crystal of plagioclase has been intensely corroded and is surrounded by 
a shell of cummingtonite, or partly by cummingtonite and partly by the hypersthene-spinel symplec- 
tite. The crystals are full of the minute, dustlike enclosures which have been observed in the feld- 
spar of almost all coronites; larger grains of spinel are present too, arranged along the cleavage and 
twinning planes. The marginal parts of many crystals are free from inclusions, but the clear feld- 
spar has the same composition as the dusty variety beside it because the twinning lamellae continue 
from one part into the other without any change in extinction angle. Measurements on the uni- 
versal stage were generally poor. The optic axial angle proved to be fairly constant, for six measure- 
ments, five of which included both optic axes, all fell within the range 2V(+) = 78°-82°. This may 
indicate either Angz or Anges. The pole of (010) gave erratic results ranging from Ang to Ang; 
since the latter value is inconsistent with the observed 2V it may be assumed that the composition 
of the plagioclase is not far from An. This agrees with 2V(+) = 80° and also with the extinction 
angle of 25° observed in sections perpendicular te (010). 

Iron ORE (probably ilmenite mostly) forms less than 5 per cent of the rock. It appears as ir- 
regular scraps and skeleton forms which are almost always surrounded by deep-brown cummington- 
ite or red-brown biotite. The ore is almost entirely oxidic, although in one case a few infinitesimal 
specks of sulphide were observed. 

BrotireE in relatively large poikilitic plates is often prominent, although it is quantitatively in- 
significant. It is a deep red-brown variety with pleochroism to pale yellow brown, and it is nearly 
uniaxial. Of all the dark silicates biotite was probably the last to form. It encloses some grains of 
apatite which are not noticeable in other parts of the rock. 


CONCLUSIONS 


Coronas are not products of the magmatic period of crystallization. If they were, 
they should be as common as micropegmatite in diabase instead of being quite rare. 
A large proportion of the rocks in auvergnose and hessose should then be coronites, 
and so should some of those in argeinose, custerose, palisadose, and kilauose, since 
at least one coronite appears in each of these subrangs of the quantitative chemical 
classification. Late magmatic reactions cause enrichment with silica and alkalis; 
but equations (2), (3), and (4) show that coronas as a whole are no more siliceous 
than the olivine and anorthite that gave rise to them. Spinel contains no silica, 
and garnet has a lower ratio of silica to lime than anorthite has. As for alkalis, 
the pyroxene of coronas is not soda-bearing, and there is no indication that the 
amphibole is. The plagioclase in coronites shows no appearance of having been 
albitized. Micropegmatite and myrmekite are absent. 

The same objections apply to the idea of deuteric crystallization from hot aqueous 
solutions. Hornblende and biotite are often formed under deuteric conditions, 
but hypersthene and garnet are not. 

The characteristic minerals of coronas are mostly foreign to eruptive rocks (spinel, 
garnet, actinolite, anthophyllite, cummingtonite), but all of these, as well as hyper- 
sthene, are typical products of thermal metamorphism. As far as geological evi- 
dence goes (and admittedly it does not go far enough) we know that at Odegarden, 
Radmansé, Grovstanis, Wiesental, Le Pallet, northern Karelia, Cape Greg, and 
Gore Mountain, and in the new Quebec occurrence, the necessary condition for ther- 
mal metamorphism—the presence of younger intrusives—was realized. It is not 
extravagant to claim that thermal metamorphism, acting upon rocks of suitable 
mineralogical composition, seems to offer the most likely solution of the corona 
problem. 
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The reactions that gave rise to coronas took place between solid and solid in the 
presence of a minimal amount of water or water vapor, at a temperature above 
that of serpentinization or chloritization. Chemical equations that express double 
decomposition between fixed numbers of entire molecules are not valid under these 
conditions. The two solids, olivine and anorthite, served merely as sources of ions, 
and their proportions by weight had no influence on the course of the reaction. 
The fallaciousness of the “Becke equation,” fer example, has been demonstrated 
by actual analysis of the garnets from the coronas of Gore Mountain and Aurlands- 
fjord. 

The few chemical analyses of coronites that we have discussed make it clear 
that coronitic rocks may have a considerable range of chemical composition, pro- 
vided only that silica is low and that magnesia, lime, and iron are fairly abundant. 
In mineralogical terms this generally means a considerable content of olivine and 
plagioclase, although hypersthene and even iron ore may substitute for olivine. 
As regards the composition of the plagioclase, it is obvious that the original feldspar 
of the future coronite must have been richer in anorthite than the feldspar that 
remains after the coronas were formed. If the residual feldspar is labradorite, the 
original may have been in the bytownite-anorthite range. The richer a rock is in 
olivine and anorthite, the more sensitive it seems to be to metamorphism. . 

What becomes of the albite that is left over when the anorthite fraction of the 
plagioclase is used up? Brdégger suspected that it might enter the amphibole mole- 
cule, and he wrote an equation (6) to illustrate this possibility. Nobody has yet 
succeeded in isolating a corona amphibole for analysis, but Du Rietz observed that 
the amphibole in the Radmansé coronite is similar in optical characters to an amphi- 
bole that occurs in other parts of this rock, and he was able to separate and analyse 
this mineral. The content of, soda proved to be less than one-third of the amount 
indicated by Brégger’s formula, while alumina was present in more than twice the 
required proportion. It seems, then, that little albite can be accounted for in 
this way. 

Perhaps the solution of the problem is that as fast as the ions of anorthite (or 
the “anorthite molecules,” if we use the language of solid solution) were removed 
from the plagioclase, the ions of albite (or the ‘“‘albite molecules’’) redistributed 
themselves through the feldspar lattice so that the composition remained the same 
throughout. It is known that microperthite becomes homogeneous when held fer 
some time at such a comparatively low temperature as 350° (Spencer); and that 
the plagioclase in volcanic rocks tends to lose its zonal structure and become homo- 
geneous under thermal metamorphism. 

Study of the few available chemical analyses brings out that coronites seem 
to be relatively rich in ferrous silicate (orthosilicate in olivine and garnet, metasili- 
cate in hypersthene and hornblende). Together with the observation that coronas 
are often formed about iron ore, this suggests that ferrous iron may have an im- 
portant réle in the generation of coronas. We know that olivine in course of ser- 
pentinization liberates iron ions which are precipitated as magnetite, whereby 
the ratio of silica to magnesia is increased. It is suggested that the expulsion of 
iron ions and also of magnesium ions from olivine, which thereby converts itself 
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into enstatite, is perhaps the fundamental reaction in the formation of coronas, 
Since iron and magnesium ions travelled all the way from olivine to plagioclase, while 
calcium and aluminum ions failed to travel all the way in the reverse direction, it 
seems to be demonstrated that olivine was the most active agent in the formation 
of coronas. 

Briefly, the idea propounded here is that the formation of coronas is a conse- 
quence of the instability of olivine (perhaps especially of iron-rich olivine) under 
the conditions of thermal metamorphism. The iron and magnesium ions dis- 
charged from the olivine attacked the anorthite of the plagioclase, generating either 
garnet or amphibole and spinel. The albite of the plagioclase, it is conceived, re- 
distributed itself through the feldspar framework so as to keep the structure homo- 
geneous. A little albite may, however, have gone into amphibole, and some may 
have escaped from the system. Microscopic evidence suggests that some spinel 
was generated in place within the plagioclase crystals, especially along cleavage 
and twinning planes. The frequent appearance of lines of spinel grains and dusty 
particles which pass from corroded plagioclase into the adjoining hornblende grains 
without change of direction bears witness to an extremely slow replacement of 
plagioclase by hornblende. 

There are many other aspects of the subject that might be discussed here, but 
the urgent need is for facts, not words, especially for better field data about the 
occurrence of coronites, but also for more chemical data about the composition of 


corona minerals. 
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EXPLANATION OF PLATE 1 


Pirate 1.—CORONAS IN TROCTOLITE FROM EASTERN QUEBEC 


Figure 

1.—A grain of olivine (colorless, high relief) in the center of the field is surrounded by a shell 
of hypersthene (light gray), then by a shell of hypersthene and spinel (mottled). On the right 
the second shell stops, and its place is taken by cummingtonite (darker gray). Cumming. 
tonite is also seen below and to the left of the olivine, and one piece is enclosed in the olivine, 
Top and left side, plagioclase darkened by dusty inclusions. 

2.—In the upper right and lower left corners, olivine is surrounded by (1) hypersthene, (2) hyper- 
sthene with granular spinel, (3) hypersthene with vermicular spinel. Center, a triangular 
area of dusty plagioclase which has suffered extreme corrosion. Top left, a granular aggre- 
gate of cummingtonite grains through which there run lines of dusty inclusions parallel to 
and continuous with those in the plagioclase. 

3.—Confluence of four coronas, each having its olivine nucleus, (colorless, high relief), its hyper- 
sthene shell (light gray, prismatic); then a shell of hypersthene and spinel, common to all 
four coronas; and finally a narrow shell of cummingtonite next to the dusty plagioclase which 
is seen at top, bottom, and right side. 

Magnification in each case about 20 diameters. 
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ABSTRACT 


Wentworth and Powers have described four stages of glaciation on Mauna Kea, Hawaii. The 
type localities of the deposits have been studied by the writer, and only the latest or Makanaka 
deposits can be accepted as definitely glacial drift. The deposits of the first and second stages are 
paroxysmal explosion deposits. Those of the third stage are a fanglomerate that could have been 
laid down in early Wisconsin time by floods caused by hot lava melting the icecap. 


INTRODUCTION 


Two papers discussing the glaciation of Mauna Kea, Hawaii (Gregory and Went- 
worth, 1937; Wentworth and Powers, 1941), have adequately described the topog- 
raphy, climate, drainage, and vegetation. Mauna Kea is a composite volcanic 
mountain 13,784 feet high on the island of Hawaii. It is studded with cinder cones, 
many of which gave rise to andesitic lava flows that veneer a dome of primitive 
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olivine basalts and associated pyroclastics. A series of transition lavas lies between 
the primitive basalts and the andesites, and basalts are intercalated with the late 
andesites (Macdonald). The summit plateau suggests that the andesites fill , 
caldera several miles across (Fig. 1). 
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FicureE 1.—Section of Mauna Kea 
Showing the relation of the basalt dome to its hypothetical caldera, and a cap of later rocks chiefly andesite. 
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FIGURE 2.-—Sketch map of the summit of Mauna Kea 
Simplified from U. S. G. S. quadrangles, after Gregory and Wentworth. 


The whole top of the mountain is accessible on foot from forest cabins located 
high on the slopes. The principal gulches are Waikahalulu and Pohakuloa, both 
of which are cut several hundred feet into the southern slope of the mountain (Fig. 
2). Seven days’ field work in and near these gulches was done in February 1944. 
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THE SO-CALLED FOUR DRIFTS 


Wentworth and Powers (1941, p. 1206-1211) describe four drifts named from 
youngest to oldest, the Makanaka, Waihu, Pohakuloa, and pre-Pohakuloa drifts, 
and correlated tentatively with the Wisconsin, Illinoian, Kansan, and Nebraskan 
glacial stages in North America. They mapped the deposits in detail (Fig. 4), 
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FIGURE 3.—Composite stratigraphic section, Mauna Kea 
Shows rocks in the upper southwestern slope. 


extended the pre-Pohakuloa icecap to 6900 feet on the basis of drift at this altitude 
near Pohakuloa Gulch, and placed the lower limit of the Makanaka icecap at 10,550 
feet. Type localities for all four drifts are in Pohakuloa and Waikahalulu gulches. 
The writer has studied these localities and agrees that the Makanaka deposit is 
glacial, that it is presumably to be correlated with the Wisconsin stage, and that 
its lower limit is at 10,550 feet. It will be shown, however, that the Pohakuloa and 
pre-Pohakuloa drifts are explosion deposits and that the Waihu drift at the type 
locality is fanglomerate, possibly laid down as a flood deposit while the icecap was 
being melted by volcanic eruptions. 


STRATIGRAPHY 
POST-MAKANAKA LAVA FLOWS 


Four small lava flows have issued in Recent time in the summit area covéring 
parts of the Makanaka drift and outwash (Fig. 3). Several black flows on the 
lower slopes appear to be postglacial also. Thus no doubt exists that Mauna Kea 
has been intermittently active in Recent time. 


MAKANAKA DRIFT 


The Makanaka drift (Fig. 3) is characterized by large angular and subangular 
blocks of andesite, with a yellow or light-gray patina, set in an indurated light-gray 
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matrix. The drift is well exposed at 10,5U0 feet in the amphitheater head of Poha. 
kuloa Gulch. Bedding is virtually absent except in scattered lenses a few feet 
long and 1-2 feet thick. The drift is sufficiently compacted to form a sloping bluf 
50 feet high along the rim of the canyon. Blocks are readily removed from the 
matrix with a hammer. Very few of the blocks show glacial polishing and striae. 
Highly vesicular blocks of olivine basalt such as are derived from primitive-type 
pahoehoe flows are notably absent. In the upper few feet weathering has stained 
the surfaces of the blocks yellow with limonite, and the matrix has changed to 4 
yellow ochre color. The surface of the deposit is covered with loose blocks oxidized 
but rarely pitted except in a small area 2 miles northeast of Puu Keonehehee, where 
a weathered pitted aa flow supplied most of the blocks. 


MAKANAKA OUTWASH 


The Makanaka outwash (Fig. 3) is either unconsolidated or weakly consolidated 
and lies as an incomplete veneer below the terminal moraine. Its pale yellow and 
gray patined andesitic blocks reaching 4 feet in diameter distinguish it readily from 
older deposits. The matrix is sandy and light gray except where stained yellow 
with limonite. 

Large fans of unconsolidated and weakly consolidated coarse alluvium carrying 
light-colored andesitic and basaltic boulders up to several feet in diameter lie at 
the mouths of Waikahalulu and Pohakuloa gulches where they debouch on the 
plain at 6500 feet. The top seems to be chiefly Recent alluvium, but the bulk of 
the fan is probably Makanaka outwash (Jaggar, 1925). Veneering much of the 
southwestern slope between these fans and 8000 feet are deposits of unconsolidated 
and weakly consolidated poorly sorted gravel and black volcanic ash. Some of the 
debris is Recent, but much of it was probably transported by meltwater from the 
Makanaka icecap in late Wisconsin time. Some of it appears to be reworked Waihu 
fanglomerate. Wentworth and Powers (1941) mapped some of the debris as un- 
differentiated drift, but all the deposits are clearly water-laid. The few large 
boulders are no larger than those in the adjacent alluvial fans or in the beds of 
Pohakuloa and Waikahalulu streams. 


WAIHU FANGLOMERATE 


The Waihu fanglomerate (PI. 1, fig. 2) is the most easily recognizable sedimentary 
deposit on the mountain because of its great induration. It is the only deposit 
that sheds blocks of conglomerate which form boulders in the later gravel beds and 
in the present gulches. Its purplish-gray matrix consists of ash and silt, except in 
a few outcrops where its original surface has not been eroded. In these outcrops, 
the matrix is stained yellow with limonite, but the stain does not penetrate more 
than 1-2 feet into the matrix. The deposit is distinctly bedded and cross-bedded 
to the bottom in all outcrops examined. It contains subrounded, angular, and 
subangular blocks up to 4 feet in diameter, many of which are andesites similar to 
those in the Makanaka drift. Vesicular olivine pahoehoe boulders are scarce or 
absent. Most blocks are smooth to the touch. 

The deposit has been considerably eroded so that in places pinnacles of the con- 
glomerate project through lag gravel and boulders. It forms vertical canyon walls 
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50 feet high on the rim of Pohakuloa Gulch and Gulch B not far below the Makanaka 
moraine (Pl. 1, fig. 2). Blocks can seldom be removed from the matrix with a 
hammer. Blocks of the conglomerate 5 to 15 feet across tumble into the gulches 
unbroken. The large lag boulders on its surface are usually weathered brown 
and pitted. Evidently these were mistaken by Wentworth and Powers for erratics. 

The Waihu deposits form fans 100 feet thick between altitudes of 9000 and 9750 
feet between Waikahalulu and Pohakuloa gulches. A pronounced break in slope 
occurs at this level which would cause heavily laden flood waters to drop their load. 
A dark-brown fringe appears to border the lower edge of the light-colored Makanaka 
moraine. It was described by Wentworth and Powers (1941) as the Waihu drift. 
Such a fringe is clearly visible from the plain below, but it is due to a change from 
green shrubs to bare rock and does not represent the termination of the Waihu 
fanglomerate, which has been traced all the way to the plain. 

The Waihu fanglomerate thins below the 9000-foot level, and at 6900 feet it is 
usually only a few feet thick. It retains all its characteristics at this low level except 
that the pebbles and cobbles are usually only a few inches across. In a few places, 
however, lag boulders several feet across lie on its surface. 

At many places on the upper slopes the Waihu fanglomerate is overlain by pre- 
Makanaka lava flows (Fig. 5). One fills a shallow valley cut into the fanglomerate 
0.3 mile east of Pohakuloa Gulch. At no place could the fanglomerate be found 
passing under Makanaka drift, although in several places it is overlain by Makanaka 
outwash. 

The Waihu fanglomerate is older than the late Makanaka outwash, and as it 
forms the rim of both walls of Pohakuloa Gulch it antedates the gulch. Thus 
Pohakuloa Gulch was cut chiefly by meltwaters of the Makanaka glacier, probably 
during the wasting of the icecap. The firm cementation of the Waihu fanglomerate 
is not conclusive evidence of great age, as the older underlying Pohakuloa explosion 
debris is less well cemented. The cementation seems more likely to be due chiefly 
to the composition of the matrix. 

The Waihu fanglomerate with its well-developed bedding (Pl. 1, fig. 2) might 
be outwash of an earlier glacial stage, but it is certainly not drift as described by 
Wentworth and Powers. The volume of Waihu fanglomerate is vastly greater than 
that of the Makanaka outwash. If the fanglomerate originated as outwash from 
an earlier glacial stage, it would indicate a much longer period of glaciation than 
the Makanaka stage. If such a glacial stage occurred, the icecap must have been 
smaller than the Makanaka icecap, otherwise an ancient terminal moraine would be 
found below the Makanaka moraine in this area. Mauna Kea was intermittently 
active as a volcano during late Pleistocene and early Recent time as shown by the 
flows above and below the Waihu and Makanaka deposits. The hypothesis is ad- 
vanced that the Waihu fanglomerate was laid down by floods caused by eruptions 
melting the icecap in early Wisconsin time, thereby releasing great volumes of water 
loaded with glacial debris. Such floods and the associated mudflows are common 
in Iceland and the deposits commonly become firmly cemented soon after coming 
to rest. The terminal boundary of an icecap melted by a volcanic eruption does 
hot necessarily coincide with the glacier that reforms. It is therefore possible that 
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FIGURE 4.—Idealized areal geologic map by Wentworth and Powers 
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Ficure 5.—Idealized areal geologic map by Stearns 


Same area and scale as Figure 4. 


small areas of drift may be found close to the Makanaka terminal moraine that are 
Waihu in age. 

Negative evidence supporting early Wisconsin age for the Waihu deposits is the 
absence of shore lines above the 25-foot level on the pre-Waihu lavas on the dry 
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Ficure 1 
Beppep Exp.iosion Deposits 
About 50 feet thick. overlain by thin- 
bedded lavas capped with Waihu fan- 
glomerate, Pohakuloa Gulch, altitude 
10,200 feet. 
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Figure 2.—Beppep Wainu FANGLOMERATE 
60 feet thick, Guich B, altitude 9500 feet. Note lag gravel on surface. 


EXPLOSION AND FANGLOMERATE DEPOSITS, MAUNA KEA 
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Ficure 1.—DEeETAILED VIEW OF 
VesIcUuLAR BASALTIC ExPLosION 
Desnris 
Waikahalulu Gulch, altitude 9600 feet. 


. 


FicurRE 2.—ERosionAL UNCONFORMITY IN Explosion DEBRIS 
Gulch B, altitude 9000 feet. 


EXPLOSION DEPOSITS, MAUNA KEA 
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western coast of Mauna Kea near Kawaihae, in contrast to a full suite of shore lines 
extending more than 250 feet above sea level on the adjacent older lavas of Kohala 


Mountain. 


EXPLOSION DEPOSITS IN WAIKAHALULU GULCH 


The Pohakuloa “drift” in Waikahalulu Gulch, the type locality, is weakly con- 
solidated explosion debris 90 feet thick derived entirely from olivine basalts, chiefly 
of pahoehoe type. It consists of rock fragments rudely sorted by size in each bed 
in a uniform gray matrix of vitric-crystal ash. Lapilli beds are common. The 
absence of erratics, striated pebbles, and soled blocks, and the uneroded contact 
between the deposit and the underlying friable cinders all deny glacial origin and 
show clearly that the debris was laid down rapidly by a series of paroxysmal ex- 
plosions (Pl. 2, fig. 1). All the blocks are subangular, have a powdery hammered 
surface from attrition, and are rough to the touch in contrast to the smooth surfaces 
on the glacial blocks in the Makanaka and Waihu deposits. Francois Matthes, 
ina letter dated July 3, 1944, has pointed out that the rasping feel of volcanic ejecta 
in situ enables one to readily separate them from iceworn, streamworn, or wind- 
worn boulders. 


Absence of andesite blocks in this deposit only 40 feet below the Waihu and 


Makanaka deposits, both of which contain chiefly andesite blocks, proves that this 
explosion deposit antedates the andesitic eruptions. If the writer’s interpretation 
of the age of the Waihu and Makanaka deposits is correct, the andesites in the 
summit area are all late Pleistocene and Recent in age. 


EXPLOSION DEPOSITS IN POHAKULOA GULCH 


The pre-Pohakuloa drift in the type locality in Pohakuloa Gulch is explosion 
debris (Pl. 1, fig. 1) similar in most respects to the one in Waikahalulu Gulch, except 
that it contains blocks of a porphyritic basalt rich in olivine and augite phenocrysts. 
These blocks seem to indicate a different explosion. If such explosion beds are 
common in the olivine-basalt series, they would account for the rapid increase in 
slope from 5° in the lower part of the mountain to 15° or 20° in the upper part. 

Similar explosion deposits in Gulch B, near Pohakuloa Gulch, were trenched by 
a narrow valley at least 60 feet deep as shown by an unconformity in the deposits 
at about 9050 feet (Pl. 2, fig. 2). This gulch was filled with explosion debris, chiefly 
water-laid, prior to being covered with about 90 feet of lava. Some of the blocks 
in this water-laid deposit have smooth surfaces. 


ASHY HILLWASH DEPOSITS 


Lenses from 1 to 8 feet thick and a few feet to several hundred feet long are ex- 
posed between lavas of all ages in the southwestern slope of the mountain. They 
are readily identified by the thin bedding, abundant vitric ash commonly containing 
cinders, poor sorting, weak consolidation, and boulders that show little diversity in 
tock types. These deposits were formed in gullies and swales by sheetwash during 
heavy rains and are not related to glacial or explosion deposits. Such deposits are 
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commonly interbedded with the lavas of steep semiarid Hawaiian volcanoes, ¢. 
pecially with those of the dying andesitic third phase of volcanism (Stearns, 1940, 
p. 1947-1948). They perch water in the overlying lava at the Waihu Spring; 
At the lower Waihu Springs the same bed of hillwash is underlain by the Waihy 
fanglomerate. Thus none of the so-called glacial springs on Mauna Kea described 
by Wentworth and Powers (1943) appear to be perched by drift. 


SUMMARY AND CONCLUSIONS 


During Wisconsin time Mauna Kea was covered with an icecap several hundred 
feet thick that extended to an altitude of 10,500 feet and left a distinct terminal 
moraine and small outwash fans. The Waihu fanglomerate may have been deposited 
by floods caused by rapid melting of the icecap by lava flows in early Wisconsin 
time. If not, it is outwash of a pre-Makanaka glacier. All older so-called drift 
deposits are definitely explosion deposits, possibly hydromagmatic, due to hot 
lava breaking through the icecap, but they are not indicative of early glaciations, 
All so-called drift deposits and erratics, supposedly indicating that early icecaps 
extended as low as 6900 feet on Mauna Kea, are believed to be conglomerates or lag 
boulders from eroded stream-laid deposits. The absence of deposits of early Pleisto- 
cene glaciers is best explained by the hypothesis that the volcano had not at that 
time reached a height sufficient to nourish glaciers in a tropical climate. An alter- 
nate hypothesis is that the deposits lie deeply buried by subsequent lavas. 
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ABSTRACT 


The composition of the stream system of a drainage basin can be expressed quantitatively in terms 
of stream order, drainage density, bifurcation ratio, and stream-length ratio. , 

Stream orders are so chosen that the fingertip or unbranched tributaries are of the 1st order; 
streams which receive 1st order tributaries, but these only, are of the 2d order; third order streams 
receive 2d or ist and 2d order tributaries, and so on, until, finaily, the main stream is of the highest 
order and characterizes the order of the drainage basin. 

Two fundamental laws connect the numbers and lengths of streams of different orders in a drainage 


sin: 

(1) The law of stream numbers. This expresses the relation between the number of streams of a 
given order and the stream order in terms of an inverse geometric series, of which the bifurcation ratio 
ryis the base. 

(2) The law of stream lengths expresses the average length of streams of a given order in terms 
of stream order, average length of streams of the ist order, and the stream-length ratio. This law 
takes the form of a direct geometric series. These two laws extend Playfair’s law and give it a 
quantitative meaning. 

The infiltration theory of surface runoff is based on two fundamental concepts: 

(1) There is a maximum or limiting rate at which the soil, when in a given condition, can absorb 
rain as it falls. This is the infiltration-capacity. It is a volume per unit of time. 

(2) When runoff takes place from any soil surface, there is a definite functional relation between 
the depth of surface detention 62, or the quantity of water accumulated on the soil surface, and the 
rate g. of surface runoff or channel inflow. 

For a given terrain there is a minimum length x, of overland flow required to produce sufficient 
runoff volume to initiate erosion. The critical length x, depends on surface slope, runoff intensity, 
infiltration-capacity, and resistivity of the soil to erosion. This is the most important single factor 
involved in erosion phenomena and, in particular, in connection with the development of stream 
systems and their drainage basins by aqueous erosion. 

The erosive force and the rate at which erosion can take place at a distance « from the watershed 
line is directly proportional to the runoff intensity, in inches per hour, the distance x, a function of the 
slope angle, and a proportionality factor K., which represents the quantity of material which can be 
torn loose and eroded per unit of time and surface area, with unit runoff intensity, slope, and terrain. 

The rate of erosion is the quantity of material actually removed from the soil surface per unit of 
time and area, and this may be governed by either the transporting power of overland flow or the 
actual rate of erosion, whichever is smaller. If the quantity of material torn loose and carried in 
suspension in overland flow exceeds the quantity which can be transported, deposition or sedimenta- 
tion on the soil surface will take place. 

On newly exposed terrain, resulting, for example, from the recession of a coast line, sheet erosion 
occurs first where the distance from the watershed line to the coast line first exceeds the critical length 
x, and sheet erosion spreads laterally as the width of the exposed terrain increases. Erosion of such 
a newly exposed plane surface initially develops a series of shallow, close-spaced, shoestring gullies or 
rill channels. The rills flow parallel with or are consequent on the original slope. As a result of 
various causes, the divides between adjacent rill channels are broken down locally, and the flow in 
the shallower rill channels more remote from the initial rill is diverted into deeper rills more closely 
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adjacent thereto, and a new system of rill channels is developed having a direction of flow at an angle 
to ~ initial rill channels and producing a resultant slope toward the initial rill. This is called croge. 
grading. 

With progressive exposure of new terrain, streams develop first at points where the length of over. 
land flow first exceeds the critical length x., and streams starting at these points generally become the 
primary or highest-order streams of the ultimate drainage basins. The development of a rilled 
surface on each side of the main stream, followed by cross-grading, creates lateral slopes toward the 
main stream, and on these slopes tributary streams develop, usually one on either side, at points 
where the length of overland flow in the new resultant slope direction first exceeds the critical length 
Xe. 
Cross-grading and recross-grading of a given portion of the area will continue, accompanied in 
each case by the development of a new order of tributary streams, until finally the length of overland 
flow within the remaining areas is everywhere less than the critical length x. These processes fully 
account for the geometric-series laws of stream numbers and stream lengths. 

A belt of no erosion exists around the margin of each drainage basin and interior subarea while the 
development of the stream system is in progress, and this belt of no erosion finally covers the entire 
area when the stream development becomes complete. 

The development of interior divides between subordinate streams takes place as the result of com- 
petitive erosion, and such divides, as well as the exterior divide surrounding the drainage basin, are 
generally sinuous in plan and profile as a result of competitive erosion on the two sides of the divide, 
with the general result that isolated hills commonly occur along divides, particularly on cross divides, 
at their junctions with longitudinal divides. These interfluve hills are not uneroded areas, as their 
summits had been subjected to more or less repeated cross-grading previous to the development of 
the divide on which they are located. 

With increased exposure of terrain weaker streams may be absorbed by the stronger, larger streams 
by competitive erosion, and the drainage basin grows in width at the same time that it increases in 
length. There is, however, always a triangular area of direct drainage to the coast line intermediate 
between any two major streams, with the result that the final form of a drainage basin is usually ovoid 
or pear-shaped. 

The drainage basins of the first-order tributaries are the last developed on a given area, and such 
streams often have steep-sided, V-shaped, incised channels adjoined by belts of no erosion. 

The end point of stream development occurs when the tributary subareas have been so completely 
subdivided by successive orders of stream development that there nowhere remains a length of 
overland flow exceeding the critical length x-. Stream channels may, however, continue to develop 
to some extent through headward erosion, but stream channels do not, in general, extend to the 
watershed line. 

Valley and stream development occur together and are closely related. At a given cross section 
the valley cannot grade below the stream, and the valley supplies the runoff and sediment which 
together determine the valley and stream profiles. As a result of cross-grading antecedent to the 
de: elopment of new tributaries, the tributaries and their valleys are concordant with the parent 
stream and valley at the time the new streams are formed and remain concordant thereafter. 

Valley cross sections, when grading is complete, and except for first-order tributaries, are gen- 
erally S-shaped on each side of the stream, with a point of contraflexure on the upper portion of the 
slope, and downslope from this point the final form is determined by a combination of factors, includ- 
ing erosion rate, transporting power, and the relative frequencies of occurrence of storms and runoff 
of different intensities. The longitudinal profile of a valley along the stream bank and the cross 
section of the valley are closely related, and both are related to the resultant slope at a given location. 

Many areas on which meager. stream development has taken place, and which are commonly 
classified as youthful, are really mature, because the end point of stream development and erosion 
for existing conditions has already been reached. 

When the end point of stream and valley gradation has arrived in a given drainage basin, the 
remaining surface is usually concave upward, more or less remembling a segment of a parabaloid, 
ribbed by cross and longitudinal divides and containing interfluve hills and plateaus. This is called 
a “graded” surface, and it is suggested that the term “peneplain” is not appropriate, since this surface 
is neither a plane nor nearly a plane, nor does it approach a plane as an ultimate limiting form. 

The hydrophysical concepts applied to stream and valley development account for observed 
phenomena from the time of exposure of the terrain. Details of these phenomena of stream and 
valley development on a given area may be modified by geologic structures and subsequent geologic 
changes, as well as local variations of infiltration-capacity and resistance to erosion. 

In this paper stream development and drainage-basin topography are considered wholly from the 
viewpoint of the operation of hydrophysical processes. In connection with the Davis erosion cycle 
the same subject is treated largely with reference to the effects of antecedent geologic conditions and 
subsequent geologic changes. The two views bear much the same relation as two pictures of the same 
object taken in different lights, and one supplements the other. The Davis erosion cycle is, in effect, 
usually assumed to begin after the development of at least a partial stream system; the hydrophysical 
concept carries stream development back to the original newly exposed surface. 
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LIST OF SYMBOLS USED 


A = area of drainage basin in square miles. 
a = slope angle. 
¢ = distance from stream tip to watershed line. 
§ = depth of sheet flow in inches at the stream margin or at the foot of a slope length /,. 
§, = average depth of surface detention or overland flow, in inches, on a unit strip of length J,. 
5, = depth of sheet flow in inches at a distance x from the crest of the slope or watershed line. 
D; = drainage density or average length of streams per unit of area. 
¢ = energy expended by frictional resistance on soil surface, ft.-lbs. per sq. ft. per sec. 
E, = average erosion over a given strip of unit width and length /, per unit of time. 
é, = erosion rate or quantity of material, preferably expressed in terms of depth of solid material, 
removed per hour by sheet erosion. 
E, = total erosion = total solid material removed from a given strip of unit length per unit of time. 
f = infiltration-capacity at a given time ¢ from the beginning of rain, inches per hour. 
f. = minimum infiltration-capacity for a given terrain. 
fo = initial infiltration-capacity at beginning of rain. 
F, = erosive force of overland flow, Ibs. per sq. ft. 
Fy = tractive force of overland flow, Ibs. per sq. ft. of surface. 
F, = stream frequency or number of streams per unit area. 
i = rain intensity—usually inches per hour. 
I = index of turbulence or percentage of the area covered by sheet flow on which the flow is tur- 
bulent. 
K, = coefficient in the runoff equation where 6, is used instead of 6 as the depth of sheet flow. 
k, = proportionality factor required to convert the rate of performance of work in sheet erosion 
into equivalent quantity of material removed per unit of time. 
K; = a proportionality factor which determines the time /, required for infiltration-capacity to be 
reduced from its initial value fo to its constant value fe. 
K; = corresponding coefficient (to K;) in the equation for laminar overland flow. 
K, = constant or proportionality factor in equation expressing runoff intensity in terms of depth 6 
of overland flow. 
I, = average length of streams of order o. 
], = maximum length of overland flow on a given area. 
I, = length of overland flow or length of flow over the ground surface before the runoff becomes 
concentrated in definite stream channels. Ff 
h, lp etc. = average lengths of streams of 1st and 2d orders, etc. 
L’ = extended stream length measured along stream from outlet and extended to watershed line. 
Iy = total length of tributaries of order o. 
M = exponent in the equation: g. = K.6”, expressing the runoff intensity in terms of depth of sheet 
flow along the stream margin. 
n = surface roughness factor, as in the Manning formula. 
number of streams of a given order in a drainage basin. 


=“ 
i] 


4, = total number of streams in a drainage basin. 
Y,, Nz etc. = total number of streams of 1st, 2d orders, etc. 
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o = order of a given stream. 
qs = surface-runoff intensity—usually inches per hour. 
4: = runoff intensity in cubic feet per second from a unit strip 1 foot wide and with a slope length ,, 
p = stream length ratio/bifurcation ratio = 11/rp. 
rv, = bifurcation ratio or ratio of the average number of branchings or bifurcations of streams of a 
given order to that of streams of the next lower order. It is usually constant for all orders of 
streams in a given basin. 
r, = stream length ratio or ratio of average length of streams of a given order to that of streams of 
the next lower order. 
r, = stream length ratio using extended stream lengths. 
r, = ratio of channel slope to ground slope for a given stream or in a given drainage basin, = 5,/s,. 
R; = initial surface resistance to sheet erosion, lbs. per sq. ft. 
R, = subsurface resistance to sheet erosion, lbs. per sq. ft., or resistance of a lower surface or horizon 
of the soil to erosion after the surface layer of resistance R; is removed. 
o = supply rate = i — f. 
s = order of main stream in a given drainage basin. 
Sc = channel slope. 
Sg = resultant slope of ground surface of area tributary to a given parent stream. 
S = surface slope = fall/horizontal length. 
t, = duration of rainfall excess or time during which rain intensity exceeds infiltration-capacity, 
» = mean velocity of overland flow, feet per second. 
vz = mean velocity of overland flow at the distance x from the watershed line. 
Va = depth of depression storage on a given area, inches. 
w, = weight of runoff, including solids in suspension, Ibs. per cu. ft. 
w, = width of marginal belt of no erosion = x, sin A, where A is the angle between the direction of 
the divide and the direction of overland flow at the divide. 
%¢ = critical length of overland flow or distance from the watershed line, measured in the direction 
of overland flow, within which sheet erosion does not occur. 
The customary procedure of expressing rainfail, infiltration, and runoff in inches depth per hour, 
velocities, channel lengths, and discharge rates in foot-second units, and drainage areas in square 
miles, has been followed, with appropriate interconversion factors in the formulas. 


PLAYFAIR’S LAW 


More than a century ago Playfair (in Tarr and Martin, 1914, p. 177) stated: 


“Every river appears to consist of a main trunk, fed from a variety of branches, each running in 
a valley proportioned to its size, and all of them together forming a system of vallies, communicating 
with one another, and having such a nice adjustment of their declivities that none of them join the 
principal valley either on too high or too low a level.” 


This has often been interpreted as meaning merely that tributary streams and their 
valleys enter the main streams and their valleys concordantly. A caveful reading 
of Playfair’s law implies that he envisaged a great deal more than merely the con- 
cordance of stream and valley junctions. He speaks of the “nice adjustment” of 
the entire system of valleys and states that each branch runs in a valley “propor- 
tioned to its size.” 

Playfair did pioneer work based on ocular observations. There were available to 
him neither the results of measurements nor the hydrophysical laws necessary to their 
quantitative interpretation. It appears that the time has now come when such @ 
quantitative interpretation can be undertaken. 
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QUANTITATIVE PHYSIOGRAPHIC FACTORS 


GENERAL CONSIDERATIONS 


In spite of the general renaissance of science in the present century, physiography 
as related in particular to the development of land forms by erosional and gradational 
processes still remains largely qualitative. Stream basins and their drainage basins 
are described as “‘youthful,” “‘mature,” ‘‘old,” “poorly drained,” or ‘well drained,” 
without specific information as to how, how much, or why. This is probably the 
result largely of lack of adequate tools with which to work, and these tools must be 
of two kinds: measuring tools and operating tools. 

One purpose of this paper is to describe two sets of tools which permit an attack 
on the problems of the development of land forms, particularly drainage basins and 
their stream nets, along quantitative lines. 

An effort will be made to show how the problem of erosional morphology may be 
approached quantitatively, and even in this respect only the effects of surface runoff 
will be considered in detail. Drainage-basin development by ground-water erosion, 
highly important as it is, will not be considered, and the discussion of drainage de- 
velopment by surface runoff will mainly be confined to processes occurring outside 
of stream channels. The equally important phase of the subject, channel develop- 
ment—including such problems as those of the growth of channel dimensions with 
increase of size of drainage basin, stream profiles, and stream bends—will not be 
considered in detail. 


STREAM ORDERS 


In continental Europe attempts have been made to classify stream systems on the 
basis of branching or bifurcation. In this system of stream orders, the largest, most 
branched, main or stem stream is usually designated as of order 1 and smaller tribu- 
tary streams of increasingly higher orders (Gravelius, 1914). The smallest un- 
branched fingertip tributaries are given the highest order, and, although these streams 
are similar in characteristics in different drainage basins, they are designated as of 
different orders. 

Feeling that the main or stem stream should be of the highest order, and that 
unbranched fingertip tributaries should always be designated by the same ordinal, 
the author has used a system of stream orders which is the inverse of the European 
system. In this system, unbranched fingertip tributaries are always designated as 
of order 1, tributaries or streams of the 2d order receive branches or tributaries of the 
ist order, but these only; a 3d order stream must receive one or more tributaries of 
the 2d order but may also receive 1st order tributaries. A 4th order stream receives 
branches of the 3d and usually also of lower orders, and so on. Using this system 
the order of the main stream is the highest. 

To determine which is the parent and which the tributary stream upstream from 
the last bifurcation, the following rules may be used: 

(1) Starting below the junction, extend the parent stream upstream from the 
bifurcation in the same direction. The stream joining the parent stream at 
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the greatest angle is of the lower order. Exceptions may occur where geologic con- 
trols have affected the stream courses. 


u 
| 
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Ficure 1.—Well-drained basin 
(Cherry Creek, N. Y., quad., U. S. G. S.) 


(2) If both streams are at about the same angle to the parent stream at the junc- 
tion, the shorter is usually taken as of the lower order. 

On Figure 1 several streams are numbered 1, and these are 1st order tributaries. 
Streams numbered 2 are of the 2nd order throughout their length, both below and 
above the junctions of their 1st order tributaries. The main stream is apparently 
ac'b although it joins ad at nearly a right angle. It is probable that the original 
course of the stream was dcb, but the portion above dc was diverted by headwater 
erosion into stream ac’. The well-drained basin (Fig. 1) is of the 5th order, while the 
poorly drained basin (Fig. 2) is of the 2d order. Stream order therefore affords a 
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simple quantitative basis for comparison of the degree of development in the drainage 
nets of basins of comparable size. Its usefulness as a basis for such comparisons is 
limited by the fact that, other things equal, the order of a drainage basin or its stream 
system generally increases with size of the drainage area. 
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FIGURE 2.—Flat sandy area, poorly drained 
(Bridgeton, N. J., quad. U. S. G. S.) 


DRAINAGE DENSITY 


Figures 1 and 2 show two small drainage basins, both on the same scale; one well 
drained, the other poorly drained. These terms, well drained and poorly drained, 
while in common use in textbooks on physiography, are purely qualitative, and 
something better is needed to characterize the degree of drainage development 
within a basin. The simplest and most convenient tool for this purpose is drainage 
density or average length of streams within the basin per unit of area (Horton, 1932). 
Expressed as an equation 


(1) 


- 
Drainage density, Dg = —z 


where SZ is the total length of streams and A is the area, both in units of the same 
system. The poorly drained basin has a drainage density 2.74, the well-drained one, 
0.73, or one fourth as great. 

For accuracy, drainage density must, if measured directly from maps, be deter- 
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mined from maps on a sufficiently large scale to show all permanent natural stream 
channels, as do the U. S. Geological Survey topographic maps. On these maps 
perennial streams are usually shown by solid blue lines, intermittent streams by 
dotted blue lines. Both should be included. If only perennial streams were in- 
cluded, a drainage basin containing only intermittent streams would, in accordance 
with equation (1), have zero drainage density, although it may have a considerable 
degree of basin development. Most of the work of valley and stream development 
by running water is performed during floods. Intermittent and ephemeral streams 
carry flood waters, hence should be included in determining drainage density. Most 
streams which are perennial in their lower reaches or throughout most of their courses 
have an intermittent or ephemeral reach or both, near their headwaters, where the 
stream channel has not cut down to the water table. These reaches should also be 
included in drainage-density determinations. 

In textbooks on physiography, differences of drainage density are commonly 
attributed to differences of rainfall or relief, and these differences in drainage density 
are largely used to characterize physiographic age in the sense used by Davis (Davis, 
1909; Wooldridge and Morgan, 1937). In the poorly drained area (Fig. 2) the mean 
annual rainfall is about 30 per cent greater than in the well-drained area (Fig. 1), 
Therefore some other factor or factors are far more important than either rainfall or 
relief in determining drainage density. These other factors are infiltration-capacity 
of the soil or terrain and initial resistivity of the terrain to erosion. 


LENGTH OF OVERLAND FLOW 


The term “length of overland flow,” designated /,, is used to describe the length of 
flow of water over the ground ‘before it becomes concentrated in definite stream 
channels. To a large degree length of overland flow is synonomous with length of 
sheet flow as quite commonly used. The distinction between overland flow and 
channel flow is not so vague or uncertain as might at first appear. Overland flow is 
sustained by a relatively thin layer of surface detention. This disappears quickly— 
often in a few minutes—through absorption by the soil or infiltration after rain ends. 
Surface detention and surface runoff may, in fact, end before rain ends if, as is often 
the case, there is at the end of the storm an interval of residual rainfall having an 
intensity less than the infiltration-capacity. Channel flow is sustained by accumu- 
lated channel storage. This drains out slowly and lasts for hours or even days after 
channel inflow from surface runoff ends. 

In addition to its obvious value in various ways in characterizing the degree of 
development of a drainage net within a basin, drainage density is particularly useful 
because of the fact that the average length of overland flow /, is in most cases approx- 
imately half the average distance between the stream channels and hence is approxi- 
mately equal to half the reciprocal of the drainage density, or 


one _ (2) 


~ 2Da 


‘0 


Later it will be shown that length of overland flow is one of the most important inde- 
pendent variables affecting both the hydrologic and physiographic development of 
drainage basins. 
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In this paper it is frequently assumed for purposes of convenience that the average 
length of overland flow is sensibly equal to the reciprocal of twice the drainage density. 
From considerations of the geometry of streams and their drainage areas the author 
has shown (Horton, 1932) that the average length of overland flow is given by the 
equation 


ia : (3) 


2D14/1—(%Y 
So 


where s, is the channel or stream slope and s, the average ground slope in the area. 
Values of the correction factor or of the ratio /,/2Dz for different values of the slope 
ratio s-/S, are as follows: 





0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 
1.86 1.67 1.40 1.25 1.15 1.09 1.05 1.02 1.005 


Se/Sq 
1,/2Da 


The ground slope or resultant slope of the area tributary to a stream on either side 
is necessarily always greater than the channel slope since the ground surface has a 
component of slope parallel with and of the same order of magnitude as that of the 
stream, and in addition it has a component of slope at right angles to the stream. 

Table 4 shows the average channel and ground slopes of streams in the Delaware 
River and some other drainage basins, derived from topographic maps. The chan- 
nel slope in these instances is commonly from half to one fourth the ground slope. 
Often on an area which as a whole is nearly horizontal, but the surface of which is 
interspersed with hills, the ground slope may be and frequently is two or three times 
the channel slope. If the channel slope is less than one third the ground slope, the 
error resulting from the assumption that average length of overland flow is equal to 
the reciprocal of twice the drainage density may in general be neglected. 


STREAM FREQUENCY 


This is the number of streams, F’;, per unit of area, or 


F, = (4) 


be | oe 


where V = total number of streams in a drainage basin of A areal units. 

Values of drainage density and stream frequency for small and large drainage 
basins are not directly comparable because they usually vary with the size of the 
drainage area. A large basin may contain as many small or fingertip tributaries 
per unit of area as a small drainage basin, and in addition it usually contains a larger 
stream or streams. ‘This effect may be masked by the increase of drainage density 
and stream frequency on the steeper slopes generally appurtenant to smaller drainage 
basins. 


COMPOSITION OF DRAINAGE NET 


The term ‘drainage pattern” is used in rather a restricted sense in many books on 
physiography, implying little more than the manner of distribution of a given set 
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of tributary streams within the drainage basin. Thus, for example, with identically 
the same lengths and numbers of streams, the drainage pattern may be dendritic, 
rectangular, or radial. Neither the drainage pattern nor the drainage density, nor 
both, provide an adequate characterization of the stream system or drainage net ina 
given basin. There may be various combinations of stream numbers, lengths, and 
orders which will give the same drainage density, or there may be similar forms of 
drainage pattern with widely different drainage and stream densities. Something 
more is needed as a basis for quantitative morphology of drainage basins. The au- 
thor has therefore coined the expression ‘‘composition of a drainage net,”’ as distin- 
guished from ‘“‘drainage pattern.” Composition implies the numbers and lengths of 
streams and tributaries of different sizes or orders, regardless of their pattern. Com- 
position has a high degree of hydrologic significance, whereas pattern alone has but 
little hydrologic significance, although it is highly significant in relation to geologic 
control of drainage systems. 

Cotton (1935) and others have used the term “‘texture” to express composition of a 
drainage net as related both to drainage density and stream frequency. For quanti- 
tative purposes two terms are needed, since two drainage nets with the same drainage 
densities may have quite different numbers and lengths of streams. Numerical 
values of drainage density independent of other units are needed for various purposes. 


LAWS OF DRAINAGE COMPOSITION 


The numbers and lengths of tributaries of different orders were determined for 
the streams listed in Table 1. The numbers and the lengths of streams varied with 
the stream order in a manner which suggested 2 geometrical progression. Plotting 
the data on semilogarithmic paper it was found (Fig. 3) that the stream numbers fall 
close to straight lines, and (Fig. 4) the same is true of the stream lengths. From the 
manner of plotting, these lines are necessarily graphs of geometrical series, inverse 
for stream numbers of different orders and direct for stream lengths.! 

From the properties of geometric series, the equation of the lines giving the number 
N, of streams of a given order in 2 drainage basin can be written 
(5) 


N. o = . 
From the laws governing geometric series it is easily shown that the number N of 
streams of all orders is 


eh 6) 
%—1 





By definition, o is the order of a given class of tributaries, s is the order of the main 
stream, and 7; is the bifurcation ratio. 
The equation correlating the lengths of streams of different orders is, similarly, 


oA (7) 


i. = L rT; 





1 In the figure the lines of best fit were drawn by inspection. Somewhat more accurate lines could of course be ob- 
tained by the method of least squares or the correlation method. However, the agreement of the observed points with 
the lines located by inspection is so close in most cases that little would be gained in accuracy by the use of these methods. 
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These equations may appear formidable, but they are merely the statement in 
symbolic form of the simple algebraic laws of geometric series. Equations (5) and 
(7) are the most important and are readily solved by means of logarithms. 







Stream length- Miles 


Stream length- Miles 


| 


Ficure 4.—Relation of stream lengths to stream order in different drainage basins 


As an example, Table 2 shows the observed and computed numbers and lengths of 
streams of different orders, based on the following values of the variables: 


% = 3.12 
s =5 

y= 2.31 
l, = 0.994 


Actual stream numbers must of course be integers, while the computed numbers 
may be fractions. Some variation between the computed and observed stream 
numbers and lengths must be expected, for various reasons. Several drainage basins 
listed in Table 1 contain large lakes, and the drainage density is less than it would be 
if the lake did not exist, since there would then necessarily be a stream of the highest 
order traversing the lake bed. Lower Rondout Creek represents an incomplete 
drainage basin to which the geometric-series laws do not necessarily apply because 
it contains a stream or portions of streams of higher order than those originating 
within this particular area. For some other areas the data were derived from early 
editions of topographic maps which do not show all the low-order tributaries. In 
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order that the equations shall give correct results the drainage basin must be reason. 
ably homogeneous. This is true of the drainage basin of Esopus Creek above Olive 
Bridge, which is wholly mountainous. The drainage basin of lower Esopus Creek is, 
however, rolling and permeabie, with great differences in soil, vegetal cover, rainfall, 
and climate, as compared with the upper basin. Upper Esopus Creek drainage basin 


TABLE 2.—Observed and computed stream lengths and stream numbers 
Drainage basin of Esopus Creek above Olive Bridge, New York. 











Number of streams Average stream length (miles) 
Stream Order es | —— 
From ee By eq. (5) From — By eq. (7) 
| | 
1 90 94.75 | 0.994 0.994 
2 25 | 30.37 | 2.45 | 2.30 
3 9 9.73 5.64 5.31 
4 1 | 3.12 6.00 12.2 
5 1 | 1.00 | 29.00 28.3 








has a much higher bifurcation ratio and stream-length ratio than the lower basin 
(Table 1), and the composition of the drainage nets is quite different, as shown 
graphically on Figures 3 and 4. 

The importance of these equations lies both in their practical application and in the 
fact that they represent laws which evolve from physical processes which Nature 
follows rather closely in the development of stream systems under such diverse 
conditions as those of upper and lower Esopus Creek. The size of the drainage basin 
does not enter the equations directly. It is indirectly involved, since the order of the 
main stream would in general be higher in the larger drainage basin, for two homo- 
geneous drainage basins of different sizes. The order of the main stream is a factor 
in the equations, and the drainage basin in which the main stream is of the higher 
order will have, in general, more tributaries of a given order. 

The data given in Table 1 cover a wide range of conditions, from precipitous 
mountain areas, like upper Esopus Creek, and highly dissected areas, like those of 
Seneca and Owasco lakes, to moderately rolling and flat areas. They cover also 
drainage basins ranging in size from a few square miles up to several hundred square 
miles. 

The bifurcation ratio (Table 1, column 11) ranges from about 2 for flat or rolling 
drainage basins up to 3 or 4 for mountainous or highly dissected drainage basins. 
As would be expected, the bifurcation ratio is generally higher for hilly, well-dissected 
drainage basins than for rolling basins. 

The values of the length ratios (column 12) range from about 2 to about 3; the 
average is 2.32. 

In the examples given in Table 1, the stream lengths were measured to stream tips 
as shown on U. S. Geological Survey topographic maps. If stream lengths had been 
measured as extended to watershed lines, the resulting stream-length ratios would 
have been materially reduced. If c is the average length from the stream tip to the 
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watershed line, and /, and /, are actual average stream lengths of two successive or- 
ders, /2 being the higher, then, as computed in Table 1, 


eae (8) 


h 


If measured as extended to the watershed lines, 





sabes ° 
The quantity r; will always be less than r;. The average valve of r; for the streams 
listed in Table 1 is 2.32. If the stream lengths were extended to the watershed 
lines, this value would lie between 2.00 and 2.32. The theoretical value of r; for 
streams flowing into larger streams at right angles is 2.00, but ry will be greater for 
streams entering at acute angles, as do most streams on steeper slopes. The distance 
along the course of a stream from its mouth extended to the water shed line is called 
“mesh length.” The use of this quantity instead of actual stream length is prefer- 
able in physiographic studies, and its use leads to closer agreement with the theo- 


retical values. 


In Figures 3 and 4, the agreement between the mean lines for the different streams 


and the observed data is so close that the two following general laws may be stated 
regarding the composition of stream-drainage nets: 

(1) Law of Stream Numbers: The numbers of streams of different orders in a given 
drainage basin tend closely to approximate an inverse geometric series in which the 
first term is unity and the ratio is the bifurcation ratio. 

(2) Law of Stream Lengths: The average lengths of streams of each of the different 
orders in a drainage basin tend closely to approximate a direct geometric series in 
which the first term is the average length of streams of the 1st order. 

Playfair called attention to the ‘‘nice adjustment’’ between the different streams 
and valleys of a drainage basin but chiefly with reference to their declivities. These 
two laws supplement Playfair’s law and make it more definite and more quantitative. 
They also show that the nice adjustment goes far beyond the matter of declivities. 


TOTAL LENGTH OF STREAMS OF A GIVEN ORDER 


Since the total length of streams of a given order is the product of the average 
length and number of streams, equations (5) and (7) can be combined into an equa- 
tion for total stream length of a given order. 

The total length L, of tributaries of order o is: 


Lo = hret-erye™! (10) 


The total lengths of all streams of a given order is the product of the number of 
streams and length per stream. The number of streams is dependent on the bifurca- 
tion ratio r, and increases with stream order, while the length per stream is dependent 
on the stream length r; and decreases with increasing stream order. Thus the total 























292 R. E. HORTON—EROSIONAL DEVELOPMENT OF STREAMS 


lengths of streams of a given order should have either a maximum. or a minimyn 
value for some particular stream order. A maximum or minimum may not occy 
because the stream order required to give the maximum or minimum stream length 
may exceed the order of the main stream, in which case the total lengths of streams 
of a given order will either increase or decrease progressively with increasing stream 
order. An exception occurs where r, and 7; are equal. Then the total lengths of 
streams of all orders are the same and equal to hr3. The ratio of r; to rp is designated 
p and is an important factor in relation both to drainage composition and physio. 
graphic development of drainage basins. As will be shown later, the value of the 
ratio p = - is determined by precisely those factors—hydrologic, physiographic, 
cultural, and geologic—which determine the ultimate degree of drainage development 
in a given drainage basin. 

By summation of the total stream lengths for different orders, as given by equa- 
tion (10), the total stream length within a drainage basin can be expressed in terms 
of four fundamental quantities: /,, 0., 7», and 77. 


CHANNEL-STORAGE CAPACITY 


Natural channel storage is a principal factor in modulating flood-crest intensities 
as a flood proceeds down a system of stream channels. A knowledge of relative 
amounts of channel storage at different locations is required for various problems of 
flood routing and flood control. It can readily be shown that the normal channel 
storage at a given discharge rate in a given stream channel varies as a simple-power 
function of the stream length, usually less than the square of the stream length, and 
this relation can readily be determined from the stage-discharge relation at a gaging 
station. Equation (10) provides a means of determining total and average stream 
lengths for each stream order. From this the channel storage provided by each 
order of streams in the drainage basin can be determined, and by summation the 
total channel storage in the stream system becomes known. This illustrates the 
practical application of quantitative physiography to a variety of engineering 
problems. 

Different stream systems may have substantially the same drainage density and 
yet differ markedly in channel-storage capacity. The higher-order stream channels 
have larger cross sections and contain more channel storage per unit length than 
lower-order streams. If r;/r, is high, the greater length of larger stream channels 
may afford greatly increased channel storage per unit of drainage area as compared 
with a drainage basin with the same drainage density and a lower value of r;/ro. 


GENERAL EQUATION OF COMPOSITION OF STREAM SYSTEMS 
From equation (10) the total length of streams of a given order is: 
Lh = Lr,2-ere-1 


The total length of all streams in a drainage basin with the main stream of a given 
order s is the sum of the total lengths of streams of different orders, or: 


DL =; [ret + ty + ret Fr? + ++ nore (11) 


This equation is cumbersome and can easily be simplified. 
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Let: 
Tj 
p=" ft = pf, (12) 
te 
Lo = h petri", (13) 


Applying subscripts 1, 2, etc., to designate the total lengths of streams of different 
orders: 
TL=L+h+ih3+-::-L (14) 
and from (13): 
SL = hry-'(p'-1 +p) + ptt + --- pel) 
(15) 
= rt (1 + p + p? + p? + +++ pt) 


The term in parentheses is the sum of a geometric series with its first term unity and 
a ratio p and is equal to: 


go-i 





p-—1° 


Substituting this value in equation (15): 











=i 

ZL = hrs: ama? (16) 
The drainage density is, from (1): 
py a2 
__ A . 
Substituting the value of [Z from (16): 
hry * pt —1 (17) 
Da = uni 


This equation combines all the physiographic factors which determine the composi- 
tion of the drainage net of a stream system in one expression. Aside from its scien- 
tific interest in this respect, it can also be used to determine drainagedensity. Values 
of the factor (p* — 1)/(p — 1) can be obtained from Figure 5. 


RELATION OF SIZE OF DRAINAGE AREA TO STREAM ORDER 


Since equation (17) incorporates all the principal characteristics of the stream 
system of a drainage basin, it may be considered a quantitative generalization of 
Playfair’s law. It can be written in such a form as to give any one of the quantities 
h, Da, A, re, ri, and s when the other five quantities are known. If the ratio p < 1, 
then, for larger values of s, p* is small, and p* — 1 may be takenas —1.0. The equa- 
tion (17) may then be written: 


log [(1 — p) Da A/h] 





=1+ (18) 
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If p > 1, then, for large values of s, p* is sensibly the same as p* — 1, and pt — 1 js 
positive. This leads similarly to the equation: 


(0 — 1)DsA 


(s — 1) log rm + s log p = log j 
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FIGURE 5.—Diagram of factor p* — 1/p — 1 
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5 = 08 le = 1)D2A/h} + log 7 


19 
log r, + log p " 


The quantity log 7, is small relative to log (oe — 1) DzA and may be neglected, so 
that in either case the order of the main stream developed in a drainage basin of a 
given area A increases for larger values of s in proportion to the logarithm of the 
area A. If, for example, with given values of p, Da, re, ri, and 4, an area of 10,000 
square miles is required to develop a stream order s, then, under the same conditions, 
in a drainage basin of 100,000 square miles, the main stream would be of order one 
unit higher, and in a drainage basin of 1,000,000 square miles the main stream would 
be two units higher in order than in an area of 10,000 square miles. This shows at 
once why stream systems with extremely high orders do not occur—there is not room 
to accommodate the requisite drainage basins on the solid surface of the earth. The 
orders of the Mississippi, Amazon, and other large rivers have not been determined 
accurately, but the Mississippi River quite certainly does not exceed the 20th order. 

From equation (17) drainage density should vary inversely as the drainage area A, 
other things equal. Actually other things are not equal in drainage areas of different 
sizes, and, although the bifurcation ratio, stream-length ratio, and average length of 
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{st order streams may be the same in two drainage basins physiographically similar 
and of different sizes, the order s of the main stream will in general be larger for the 
larger drainage basin. As a result the drainage density may increase, decrease, or 
remain substantially unchanged in two similar drainage basins of different sizes. 





Average charrel slope ~ Ft. per Mile 


3 
Order of stream 
FicureE 6.—Law of stream slopes 


Neshaminy, Tohickon, and Perkiomen drainage basins. 
LAW OF STREAM SLOPES 


In addition to the various quantitative relationships between the different factors 
involved in drainage composition, expressed by equation (17), there are certain other 
quantitative relationships. The relation of stream slope to stream order in agiven 
drainage basin is hinted in Playfair’s law. As an illustration, the slopes of streams 
of different orders in the Neshaminy, Tohickon, and Perkiomen drainage basins have 
been plotted in terms of stream order (Fig. 6), and there is a fairly definite relation- 
ship between slope of the streams and stream order, which can be expressed by an 
inverse geometric-series law. 


DETERMINATION OF PHYSIOGRAPHIC FACTORS FOR DRAINAGE BASINS 


To determine completely the composition of a stream system it is necessary to 
know: (1) the drainage area, A, (2) the order s of the main stream, (3) the bifurca- 
tion ratio rp, (4) the stream length ratio r;, and (5) the length /, of the main stream 
or preferably the average length /, of 1st order streams. If these data are given, then 
the drainage density, stream frequency, and other characteristics of the stream 
system can be determined by calculation, using the equations which have been given. 
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From equation (5): 


\ No = rp-9), 


N, s-1 = 7d. 


This shows that the bifurcation ratio 7» is equal to the number oi streams of the next 
to the highest order for the given drainage basin. 

If the stream numbers for different stream orders are plotted on semilog paper 
(Fig. 3), the bifurcation ratio r, can be determined by simply reading from the aver- 
age line the number of streams of the second highest order. 

From equation (7): 


b, = Lywrje7h, 
If o = 2, 


-~ =n. 


h 


The stream length ratio r; can therefore be obtained by dividing the average stream 
length of any order by the average stream length of the next lower order, the values 
of stream lengths being read from the diagram of stream lengths plotted in terms of 
stream order. It is preferable to use these data rather than actual measured values, 
as the number of streams of a given order—particularly the higher order streams— 
may not be exactly the normal number for the given drainage composition. Stream 
numbers can, of course, be only integers, and there may be either two, three, or four 
streams of the second highest order in a given drainage basin where there should be 
three. 

In Table 1 and Figures 3 and 4, the stream lengths and numbers of all orders were 
determined directly from topographic maps. Where this is done the order s becomes 
known directly. 

In analyzing the drainage net of a stream system it is desirable to trace, with 
different colors for each order, the stream system from the base map. When the 
higher-order streams are determined some of the lower-order streams may prove to be 
the head-water portions of higher-order streams. Figure 7 shows the drainage basin 
of Hiwassee River above Hiwassee, Georgia, with 1st order streams shown by dotted 
lines and stream orders indicated by figures. 

The determination of stream lengths and orders by direct measurement from 
maps which are on a sufficiently large scale to show all 1st order streams is so laborious 
as to be practically prohibitive except for smaller drainage basins. 

Fortunately, all the required quantities—J,, h, rs, 71, and Dz—can be determined 
from smaller-scale maps from which the lower-order tributaries are omitted. The 
maps must show correctly the streams for several of the higher orders. The order 
of the main stream is of course unknown since it is not in general known which of the 
lower orders of streams are omitted from the map. The streams shown are assigned 
orders assuming that the main stream has an unknown order s, the next lower order 
of stream shown is designated 2,and soon. The number of streams of each assumed 
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Order No of Length Aver. length 


treams = miles miles 
1 , l4c 1Z 0.49 Drainage area =826 5q.mi 
2 32 Al 1.28 
3 2 32.8 3.65 Dy, = 2.06 
A % 24.6 12.30 
170.4 


FiGuRE 7.—Drainage net, upper Hiwassee River 


order is counted, their stream lengths measured from the map, the results tabulated 
as follows and plotted as shown by F igure 8A. 


Data for Perkiomen Creek 


Order Assumed Number of Average stream 
inverse streams lengt! 
order (miles) 

Ss 1 1 20.5 

s-—1 2 2 13.75 

s—2 3 10 3.61 

s—3 4 


32 1.39 
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If it is assumed that the main stream is of the 

4th order, then /; = 1.38 miles; 

5th order, then J, = 0.50 mile; 

6th order, then /; = 0.20 mile. 
Since /, is not far from half a mile, the main stream is of the 5th order. From line 
B (Fig. 8) the number of 2d order streams is 3.15. This is the bifurcation ratio, 
From line A the lengths of 2d and 1st order streams are, respectively, 1.38 and 0.52 
miles. This gives the stream length ratio: 


1.38 
= 052 = 2.70. 


Data for at least four stream orders are required to determine the order of the main 
stream from incomplete data by this method. Care must also be used in deter. 
mining the lines A and B accurately to secure correct results. 

The values of the stream lengths as far as known are then plotted on semilog paper 
(Fig. 8A), in terms of inverse stream orders, a line of best fit drawn to represent the 
plotted points and this line extended downward to stream length unity or less. 

To determine the order of the main streain it is necessary to know the order of 
magnitude but not the exact value of the average length of streams of the ist order. 
The length /, of streams of the ist order is rarely less than a third of a mile, a value 
which is approached as a minimum limit in mountain regions with heavy rainfall, 
as in the southern Appalachians. Also it is rarely greater than 2 or 3 miles, values 
which are approached as maximum limits under some conditions in arid and semiarid 
regions. Data from which the order of magnitude of J, can be determined are always 
available from some source. In general all that is required is to know whether 
is of the order of half a mile, 1 mile, or 2 miles or more. The point at which the 
stream length shown by the line ab (Fig. 8A) extended downward has a value about 
the same as the known value of /, for the given order indicates the order of the main 
stream. 

This method for determining the order of the main stream is of limited value in 
some drainage basins, particularly large drainage basins, such as that of the Missis- 
sippi River, which are not homogeneous, and where there may be large variations 
in the length of 1st order streams in different portions of the drainage basins, so that 
the order of magnitude of J, may be difficult to determine. A small portion of a 
drainage basin, with suitable conditions of high rainfall, steep slopes, etc., may add 
several units to the value of s for the main stream, although it has little effect on the 
weighted average value of /, for the drainage basin as a whole. For basins which are 
reasonably homogeneous the method is accurate. Proof of its validity is readily 
obtained by applying this method to a drainage basin where the values of /; and the 
drainage density Dz have been determined from measurements on a map showing 
streams of all orders, but using in the determination only the data for streams of 
higher orders. This was done in preparing Figure 8, which is of the 5th order, al- 
though only data for the first four stream orders were used in the computation, it 
being assumed that /, was of an order of magnitude between 1 and 1.5. 
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dily This determination of s gives also the average length /, of 1st order streams. The 
‘the | bifurcation ratio r, and the stream-length ratio r; are determined by the slopes of the 
ving | lines A and B on Figure 8. It is not necessary to know the order of the main stream 
s of | to determine these quantities. When ro, r:, A, s, and ) are known, the drainage 
al | density can be determined by means of equation (17). 

n, it This method of determining s has the advantage that it is at least as accurate 
when applied to large as when applied to smaller drainage basins. In general, data 
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for more stream orders will be available from a map for a large drainage basin tha 


for a small basin. 
Table 3 shows the drainage composition of Neshaminy, Tohickon, and Perkiome, 
Creek stream systems, derived in the manner described, together with the drainag 


TABLE 3.—Observed and computed drainage densities, Neshaminy, Tohickon, and Perkiomen cre; 











Item 

1 Stream Neshaminy Tohickon Perkiomen 

2 Location Below Forks Point Pleasant | Near Frederic 

3 Drainage area, square miles 139.3 102.2 152.0 

Computed values: 

4 Stream order from map 5 5 5 

5 |e 35.0 33.0 27.0 

6 fr 3.45 3.00 3.15 

7 on 2.92 2.85 2.70 

8 p=n/r 0.85 0.95 0.86 

9 f(e) 3.35 4.50 3.80 
10 rs 141.6 81.0 98.4 
11 r°1/A 1.02 0.79 | 0.65 
12 (11) X (9) 3.82 3.56 2.47 
13 i 0.50 0.53 0.52 
14 (12) X (13) = Da 1.91 1.89 1.28 
15 Drainage density from map 1.60 1.91 | 1.24 











densities as computed by equation (17) and as derived from direct measurement from 
topographic maps. 

Drainage densities computed by equation (17) will usually be somewhat higher 
than those derived directly from maps if stream lengths are measured directly and 
only to the fingertips of the stream channels, because the stream lengths and mesh 
lengths are sensibly identical for higher-order streams, whereas there may be 10 to 25 
per cent or even 50 per cent difference between stream length and mesh length for 
low-order streams. In computing drainage density from values of h, 75, and 1; ob- 
tained graphically, the computed value corresponds more nearly to drainage density 
expressed in terms of mesh length than in terms of actual stream length for lower- 
order streams. 


RELATION OF GEOLOGIC STRUCTURES TO DRAINAGE COMPOSITION 


The exainples of drainage composition shown in Table 1 and in Figures 3 and 4 in 
nearly all cases represent special or abnormal conditions, such as the presence of 
large lakes in several of the drainage basins. While this table agrees well with the 
laws of stream numbers and stream lengths even with such pronounced geologic 
control of topography as that afforded by the drumlin areas in the Ganargua Creek 
drainage basin, there are other conditions where geologic controls apparently exert a 
definite influence on drainage composition. Figure 9 shows two drainage basins the 
boundaries of which are definitely fixed by geologic structures. 
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Figure 9.—Drainage patterns of Laurel Fork and Glady Fork, Cheat River drainage basin 
(From Beverly, W. Va., quad., U. S. G. S.) 
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Data of stream lengths and stream numbers in these basins are as follows: 








Order Length Number of Average length 
(miles) streams (miles) 
Basin “‘A”—Laurel Fork 
50.8 sq. mi. 1 47.80 94 0.51 
2 23.36 26 0.90 
3 4.10 + 1.02 
4 19.62 1 19.62 
94.88 
Dz = 94.88 = 1.87 
50.8 
Basin “B”—Glady Fork 
55.44 sq. mi. 1 38.42 95 0.41 
2 34.76 37 0.94 
3 6.50 6 1.08 
4 20.80 1 20.80 
100.48 
Da = 100.48 
—— = 1.82 
55.44 


In both streams the law of stream numbers is closely obeyed. The law of stream 
lengths is approximately obeyed for lower-order streams, but in both basins it is 
necessary, in order that the area should be drained, that the main stream should have 
a length sensibly equal to that of the drainage basin; this requires that the length of 
the main stream should be much greater than it ordinarily would be for a drainage 
basin of the same order, of normal form. 

One may naturally ask whether stream systems in similar terrain and which are 
genetically similar should not have identical or nearly identical stream composition. 
Data for the drainage basins of Neshaminy, Tohickon, and Perkiomen creeks in the 
Delaware River drainage basin near Philadelphia, Pennsylvania, are given in Table 
3 and on Figure 10. The drainage patterns of these streams are shown on Figure 11. 
The physiographic characteristics of these three drainage basins are closely similar. 

Tables 4 and 5 show, respectively, drainage composition of streams in the upper 
Delaware River drainage basin and drainage composition of several small tributaries 
of Genesee River in western New York. 

The streams in the upper Delaware River drainage basin are generally similar, 
with the exception of Neversink River, in topography, geology, and climate. The 
various morphologic factors for these basins are of the same order of magnitude 
although not numerically identical. The tributaries of Genesee River listed in 
Table 5 represent areas at various locations around the margins of this basin be- 
tween Lake Ontario and the New York-Pennsylvania State line and comprise a wider 
range of geologic and topographic conditions than occurs in the Delaware River 
drainage basins, and there are correspondingly greater variations in the morphologic 
factors, particularly bifurcation ratio, length of 1st order streams, and drainage 


density. 
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Nomber of streams 


Average Lergth of stream ~Miles 





Stream order *% 


Ficure 10.—Stream numbers and stream lengths 


Neshaminy, Tohickon, and Perkiomen drainage basins. 


It is found from plotting stream lengths and stream orders, subject to the limiting 
conditions already described, that both these laws are quite closely followed. De- 
partures from the two laws will, however, be observed, and if other conditions are 
normal these departures may in general be ascribed to effects of geologic controls. 
Asa rule the law of stream numbers is more closely followed than the law of stream 
lengths; Nature develops successive orders of streams by bifurcation quite generally 
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in a uniform manner, regardless of geologic controls. Stream lengths, on the other 
hand, may be definitely limited by geologic controls, such as fixed boundaries of the 
outline of the drainage basin. 


TABLE 5.—Physiographic factors—tributaries of Genesee River, western New York 











Item Sloder Gates Rus? Red Spring Stony 
No. Iter Creek Creek Creek Creek Creek | Creek 
@ (2 ™ (a) ‘sy © (wy (@ 
1 Stream Siader Crk. | Gates Crk. |Rusk Crk. | Red Crk. |Spring Crk.| Stony Crk 
2 | Location Kanaseraga | at movth | Fillmore |at mouth | at mouth | at mouth 
3 Drainage Arca Sq. Mi is.6 20.2 43.3 25.3 20.3 22.2 
4 | Width ratio 112 126 0.39 acs 144 089 
5 | Order main stream © 3 4 3 3 3 < 
GS Stream numbers 
First order 24 3z 338 14 is zs 
Second - 4 10 10 4 4 & 
Third . i 3 { 1 ! 3 
Fourth + ! f 
7 | Average stream ‘length ~miles 
First order 0.6) 0.49 0.48 429 122 oa! 
Second - 2.12 oss us 150 2.00 128 
Third - 3.15 2.15 10.50 €25 31s 325 
Fourth . 4.00 S500 
8 |Average stream slope ~ Ft/mi 
First order 271 89 17! 26 37 287 
Second + iz) 41 80 iT) 39 196 
Third . 89 39 49 7 it 169 
Fourth «+ 2s 120 
8 | Bifurcation rahe rb 4.40 310 B10 390 330 3.10 
10 | Stream length ratio TY 221 2.00 459 257 1718 Les 
1! | Total stream length mi 21.00 37.50 €025 30.25 30.00 3875 
12 | Drainage density Da 171 1.86 ial 1.20 147 173 
13 | Ratio © Wr, 0.50 064 os OGG 046 060 
14 | Latitude ~rorth 42-25 42-5 42-25 43-05 43-05 42-30 
is | Longitude ~west 71-45 77-30 7-05 T1-40 78-05 71-40 





























INFILTRATION THEORY OF SURFACE RUNOFF 
GENERAL STATEMENT 


The factors and formulas given serve as measuring tools for the quantitative 
comparison of upland features of drainage basins. Quantitative science develops 
by the correlation of observed relationships through scientific laws and principles, 
which may therefore be described as operating tools. Two principal kinds of operat- 
ing tools are needed in connection with upland erosion: (1) the laws governing the 
sheet flow of surface runoff, and (2) the laws governing (a) soil resistivity to erosion, 
(b) erosive force, (c) erosive power of sheet flow, and (d) transporting power of sheet 
flow. The first of these tools is supplied by the infiltration theory of surface runoff, 
developed by the author (Horton, 1935; 1937; 1938). Only a few salient features 
of this theory are pertinent to the present discussion. 

The infiltration theory of surface runoff is based on two fundamental concepts: 

(1) There is a maximum limiting rate at which the soil when in a given condition 
can absorb rain as it falls. This is the infiltration-capacity (Horton, 1933). 

(2) When runoff takes place from any soil surface, large or small, there is a definite 
functional relation between the depth of surface detention 6, or the quantity of water 
which accumulates on the soil surface, and the rate of surface runoff or channel 
inflow 4. 
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These two concepts, in connection with the equation of continuity or storage 
equation, form the basis of the infiltration theory. It has hitherto been assumed 
that surface runoff was some definite fraction of rain. If that were true, then all 
rains, however low their intensity, should produce runoff. This is not an observed fact. 


INFILTRATION-CAPACITY 


Infiltration-capacity, f, is governed by physical laws and processes which involve 
the simultaneous downward flow of water and the upward flow of displaced air 
through the same system of soil pores (Horton, 1940; Duley and Kelly, 1941) and 
is used in the sense of a limiting rate of flow, like the capacity of a water pipe. 

The infiltration capacity of a given terrain, including soil and cover, is controlled 
chiefly by (1) soil texture, (2) soil structure, (3) vegetal cover, (4) biologic structures 
in the soil, especially at and near the surface, including plant roots and root perfora- 
tions, earthworm, insect, and rodent perforations, humus, and vegetal debris, (5) 
moisture content of the soil, and (6) condition of the soil surface, whether newly 
cultivated, baked, or sun-cracked. Temperature is probably also a factor, although 
its effect is often masked by biologic factors, which also vary with temperature and 
season. 


The infiltration-capacity of a given area is not usually constant during rain but, — 


starting with an initial value fo, it decreases rapidly at first, then after about half an 
hour to 2 or 3 hours attains a constant value f.. The relation of the infiltration- 
capacity to duration of rain can be expressed accurately by the following equation, 
with f, fo, and f. in inches per hour: 


i =fet+(h— feye*s* (20 


where ¢ is the base of Naperian logarithms, ¢ is time from beginning of rain, in hours, 
and Ky is a proportionality factor (Horton, 1939; 1940). This equation can easily be 
derived on the assumption that infiltration-capacity is governed chiefly by the con- 
dition of the soil surface and is reduced at the beginning of rain by effects which 
result from the energy of falling rain and which operate after the manner of exhaustion 
phenomena. These effects include packing of the soil surface, breaking down of the 
crumb structure of the soil, swelling of colloids, and the washing of fine material 
into the larger pores in the soil surface. 

As an example typical of many experimental determinations of the change of 
infiltration-capacity during rain, the values of f have been computed at different 
times, ¢, from the beginning of rain for a soil with initial infiltration-capacity fy = 
2.14 in. per hour and which drops to a constant value f- = 0.26 in. per hour in 2 
hours. The quantity Ky determines the rate of change of infiltration-capacity 
during rain for a given rain intensity and in this case has the value 3.70. 


t; 0.0 0.2 0.4 0.6 0.8 1.0 1.5 2.0 
f,in. per hr.; 2.14 1.16: 0:09 0.46 0.36 0.31 0.28 0.26 


These values of f show that the infiltration-capacity drops off rapidly at first, then 
more slowly as it approaches f.. Between rains drying out of the soil and restoration 
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of crumb structure leads to restoration of f toward or to its initial value. The geo. 
morphic significance of this decrease of infiltration-capacity during rain is illustrated 
by the fact that if infiltration-capacity remained constant at the high value it usually 
has at the beginning of rain, then there would be little surface runoff or soil erosion 

It is the minimum value f. of infiltration-capacity which predominates during 
most of long or heavy rains which are chiefly effective in producing floods and sheet 
erosion. Infiltration-capacity is highest and least variable for pure coarse sands, 
For terrain with such soils, the infiltration-capacity may always exceed the rain 
intensity, so that even in winter, when the ground is frozen, no surface runoff occurs, 
The infiltration-capacity concept accounts for the absence of sheet erosion and the 
meager development of drainage in many sandy regions, even with abundant rainfall, 

“Transmission-capacity” is the volume of flow per unit of time through acolumn 
of soil of unit cross section, with a hydraulic gradient unity or with a hydraulic head 
equal to the length of the soil column. Infiltration-capacity and transmission- 
capacity are related, but they are not identical; the infiltration-capacity is usually 
less than the transmission-capacity. Under conditions where transmission-capacity 
prevails, the soil column is fully saturated, and the entire cross section of the void 
space participates in hydraulic flow. When infiltration-capacity prevails, the soil 
column is not usually fully saturated, and air must escape upward through the soil 
as fast as water flows downward into the soil. This occupies a fraction, though 
usually only a small fraction, of the pore space. Also, because of the surface reduc- 
tion of infiltration-capacity, already described, the soil surface cannot usually absorb 
water as fasi as water can flow downward through the interior of the soil mass. Asa 
result the soil is not saturated appreciably above its capillary-capacity, even during 
the heaviest rains. This fact can readily be verified during a rain by picking upa 
handful of garden soil at a depth of a few inches. Asa rule, no water can be squeezed 
out of the moist soil although, if saturated, water could readily be squeezed out. 
A common fallacy is the statement that during a storm a severe flood was produced 
by the soil becoming saturated. This happens only in the case of heavy clay soils. 
Even if the water table rises to the soil surface, as it sometimes does in swampy areas, 
the soil can still absorb water, which appears later, somewhere, as ground-water or 
wet-weather seepage. As an example, in the central New York flood of July 1935, 
analyses showed that, in the Cayuga and Seneca drainage basins, the terrain main- 
tained an infiltration-capacity of 0.2 to 0.3 inch per hour over the areas where the 
most intense rainfall and runoff occurred. 

Infiltration-capacity of a given terrain can be determined in several ways, either 
with fully controlled conditions or for a drainage basin as a whole, under natural 
conditions. Thousands of determinations of infiltration-capacity have been made. 
These and other similar data, when more fully analyzed and classified, will, it is 
believed, form one of the most important tools for a quantitative study of drainage- 
basin morphology. 

Rain falling at an intensity i which is less than f will be absorbed by the soil surface 
as fast as it falls and will produce no surface runoff. The rate of infiltration is then 
less than the infiltration-capacity and should not be designated “‘infiltration-capac- 
ity.” If the rain intensity i is greater than the infiltration-capacity f, rain will be 
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absorbed at the capacity rate f; the remaining rain is called “rainfall excess.” This 
accumulates on the ground surface and for the most part produces runoff, and the 
difference between rain intensity and infiltration-capacity in such a case is denoted 
by « and designated the supply rate (s = 7 — f). For a constant rain intensity i, 
in inches per hour, the runoff intensity q., in inches per hour, approaches the supply 
rate o asympototically as a maximum or limiting value as the rain duration increases 
(Horton, 1939; Beutner, Gaebe, and Horton, 1940). The total surface runoff is ap- 
proximately equal to the total supply o/., where /, is the duration of rainfall excess. 


OVERLAND OR SHEET FLOW 


In the minds of most persons the term sheet flow probably implies a greater depth 
of flow than usually occurs. Sheet erosion is used in contradistinction to channel 
erosion, and the use of sheet flow to describe overland flow not concentrated in 
channels larger than rills is appropriate, but it may not imply flow to depths mea- 
sured in feet or even in inches but rather in fractions of an inch. 

Since 1 inch per hour equals approximately 1 second-foot per acre or 640 c.s.m.?, 
and an acre is 208 feet square, the surface-runoff intensity g, in cubic feet per second 
from a unit strip 1 foot wide and a slope length /, will be: 


g1 = 0.000023 logs (21) 


where g, is the runoff intensity in inches per hour. Discharge = Depth X Velocity» 
or if 6 is the depth of sheet flow, in inches, and » the velocity in feet per second, g; = 
05 
1D It follows that the depth of sheet flow at any point on a slope where the slope 
length is 7, will be: 
5 = 0000277104. 
ree 


(22) 


Ona gently sloping lawn, with a length of overland flow of 100 feet and a velocity of a 
quarter of a foot per second, a depth of surface detention of 0.11 inch will produce 1 
inch runoff per hour. Walking over such a lawn while this runoff intensity is occur- 
ring, one may not notice that surface runoff is taking place; yet it is this same un- 
obtrusive and almost imperceptible overland flow which, with greater depths and 
larger volumes and on longer slopes, is largely responsible for carving the landscape 
of drainage basins into observed forms. 


LAW OF OVERLAND FLOW 


The velocity of turbulent hydraulic flow is expressed in terms of the Manning 
formula: 


sa = ef} (23) 
n 


where 9 is the mean velocity in feet per second, » is the roughness factor, having the 
same general meaning for sheet flow as for channel flow, R is the hydraulic radius or 





*¢.s.m. = cubic feet per second per square mile. 
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ratio of area of cross section to wetted perimeter. For sheet or overland flow, R 
becomes identical with the depth 6. S is the slope. Since discharge or volume of 
flow per time unit per unit width equals the product: Velocity X Depth X Width, 
the runoff intensity in inches per hour from a strip of unit width, for turbulent flow, 


can be expressed by: 


qe = Ka (24) 


where K, is a constant for a given strip of unit width, having a given slope, roughness, 


and slope length. 
A similar equation: 


qs = K,6*S (25) 


can be derived from Poiseuille’s law for nonturbulent or laminar flow. 

Overland flow may be either wholly turbulent, wholly laminar, or partly turbulent 
and partly laminar—patches of laminar flow being interspersed with turbulent flow 
or vice versa. Since the equations for turbulent and laminar flow are of the same 
form, it follows that for either laminar or turbulent flow, or for mixed flow, the rela- 
tion between depth of surface detention and runoff intensity, in inches per hour, 
should be a simple power function of the depth of surface detention or: 


aaa K.6™” (26) 


where gq, is the runoff intensity in inches per hour, 6 is the depth of surface detention 
at the lower end of the slope, in inches, K, is a coefficient involving slope, length of 
overland flow, surface roughness, and character of flow, and the exponent M has a 
value of 5/3 for fully turbulent flow. 

Except for very slight depths of surface detention, this simple law of surface runoff 
is remarkably well verified by plot experiments (Fig. 12). The circles (Fig. 12) 
indicate points derived directly from the hydrograph, and the solid lines the resulting 
relation curves plotted logarithmically. The points fall almost precisely on the 
relation lines, indicating an accurate functional relationship between 6, and q,. 

Except on steep slopes there are always depressions, often small but numerous, 
on a natural soil surface. If the derived points were plotted for smaller depths than 
those shown on Figure 12, the corresponding relation lines would curve off to the 
left, indicating that the power-function relation of g, to detention depth changes for 
very slight depths of surface detention. This represents the effect of depression 
storage. When runoff is taking place, flow through the depressions also occurs, 
although usually slowly, and hence the full cross sections of the depressions partic- 
ipate in determining the law of overland flow. The runoff becomes zero, however, 
when the depth of detention is reduced to the depth of depression storage Va, although 
water still remains in the depressions. Consequently for slight depths the relation 
lines curve to the left. 





3 Expressed hydraulically as the ratio: fall/horizontal length. For steeper slopes the sine of the slope angle should be 
used in place of S. 
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INDEX OF TURBULENCE 


The exponent in the surface-runoff equation is not the same for all the lines (Fig. 
12). This exponent, as indicated by the slopes of the lines, came out 2.0 for the line 
at the left, 1.48 and 1.47 for the next two lines, and 1.09 and 1.00 for the lines at the 
right. This suggests that in these cases the flow was not wholly ordinary turbulent 
flow. 

For fully turbulent flow the exponent M should be 5/3, while for fully laminar flow 
the expenent M should be 3.0. 

Where the exponent M > 5/3, the larger exponent can readily be explained if 
part of the overland flow is laminar flow, and this is quite certain to occur where the 
flow is alternately in thin films on ridges in the form of laminar flow, and through 
depressions wholly or partly as turbulent flow. An index of turbulence applicable 
to such a case is expressed by the equation: 


I=}(3— M). (27) 


If the flow is fully turbulent and M = 5/3, this gives J = 1.0. If the flow is fully 
laminar and M = 3.0, this gives J = 0. 
By transposition, the runoff exponent M can be expressed in terms of the index 
of turbulence: 
M =3-41. (28) 


Turbulence in overland flow increases downslope from the watershed line. In 
laboratory and field-plot experiments with plot lengths /, usually less than 25 feet, 
the flow over surfaces without vegetal cover is usually-partially turbulent. On long 
natural slopes, with /, much greater—frequently 1000 feet or more—the flow is fully 
turbulent except for extremely slight depths or close to the head of the slope. 


TYPES OF OVERLAND FLOW 


The study of overland flow in accordance with the infiltration theory has revealed 
various phenomena of microhydraulics not commonly present in ordinary channel 
flow. 

Partially turbulent flow described may be considered “mixed flow.” In general it 
consists of turbulent flow interspersed with laminar flow. 

If the area on which flow occurs is covered with grass or other close-spaced vege- 
tation, the flow may be “subdivided.” Part of the energy available for overcoming 
resistance is expended on the grass blades and stems, reducing the amount of energy 
available for expenditure on the soil surface. For the limiting condition of complete 
subdivision of the flow, all the resistance to flow would be due to the vegetation, and 
the law of overland flow would be: 


a = Kaba. (29) 


The velocity of overland flow would be sensibly constant regardless of the depth of 
surface detention. Some experiments show substantially this condition. Because 
of the increased resistance, the depth of surface detention required to carry a given 
rate of runoff is very greatly increased, and the velocity of overland flow is corre- 
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spondingly decreased where there is dense cover of grass, grain, or similar vegetation. 
This is an important fact in relation both to surface runoff and soil erosion. 

It is often found in runoff-plot experiments that the hydrograph does not have a 
smooth surface but is broken into irregular waves or surges. This type of flow may 
be designated “‘surge” flow and may be due to several causes: 

(1) Under certain hydraulic conditions steady flow cannot occur even on a smooth, 
unchanging surface (Jeffreys, 1925). 

(2) Plant debris, especially of the sand-burr type, may be loosened and carried 
along with the flow, forming debris dams, behind which the water piles up, and these 
hold back the water temporarily and then release it in relatively large volumes, 
producing irregular waves (Beutner, Gaebe, and Horton, 1940). 

(3) Active surface erosion may produce a succession of irregular waves due either 
to mud or mud-and-debris dams of the type last described, to the breaking down of 
divides between natural depressions, or to the lateral incaving of the walls of gullies 
(Horton, 1939). Erosion may produce traveling mud dams or mud flows similar to 
those sometimes produced on a larger scale in mountain canyons by cloudburst 
storms. Wherever a mud or debris dam is formed, water accumulates behind it 
until presently the dam moves down the slope with the accumulated water behind it. 
In case of surge flow or traveling back-water due to debris dams, there is often no 
consistent relation between depth of overland flow and runoff intensity. 


RAIN-WAVE TRAINS 


On slopes which are not too flat, shallow flow in the form of a uniform sheet may be 
hydraulically impossible. The flow then takes the form of wave trains or series of 
uniformly spaced waves in which nearly all the runoff is concentrated. The author 
has twice observed such rain-wave trains in intense storms. They occur most com- 
monly in rains of high intensity, particularly those of the cloudburst type, char- 
acterized by their ability to tear up sod on slopes and carry fences and other large 
debris into stream channels. Rain-wave trains occur only under suitable hydraulic 
conditions and have been described in another paper (Horton, 1939). Observation 
strongly indicates that rain-wave trains may be important in initiating erosion on 
sloping lands. If, for example, the waves are 6 feet apart, then each wave contains 
as much water as would be contained in a length of slope of 6 feet with uniform flow. 
The successive waves, with their concentration of runoff and energy, can strike 
sledge-hammer blows on obstructions. They may initiate erosion where it would 
never occur from the same runoff intensity with steady flow. The difference between 
the two cases is like that of breaking a rock with a few sledge-hammer blows, when a 
million taps with a pencil tip would expend the same amount of energy but produce 
no effect. 

The following types of sheet or overland flow take place: 


Type of flow M I 
RR ee bai ia Naas 3 0 
Mixed laminar and turbulent....................... 3t05/3 Otol 
Turbulent.......... Nd dia Cite aiGilcak wees 5/3 1.0 
Subdivided or superturbulent....................... 5/3 to1 1 to 14 
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PROFILE OF OVERLAND FLOW 


It can readily be shown from the infiltration theory that the profile of sheet or 
overland flow, or the relation of depth 6 of surface detention to the distance « down- 
slope from the watershed line, is expressed by a simple parabolic or power function 


Outlet 


Ficure 13.—Half section of a small drainage basin 
Vertical scale greatly exaggerated. 
A similar power function 





Illustrating runoff phenomena. 
(Horton, 1938). This relation is illustrated on Figure 13. 
expresses the relation of velocity of overland flow in terms of distance from the water- 
For turbulent flow: 
(30) 


oe x 


3/5 
a & i) 
(31) 
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shed line (Horton, 1937). 


and, in general, for 6, in inches, for any type of flow: 
(32) 
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where ¢ = supply rate, in inches per hour; /, = total length of slope on which over- 
land flow occurs, in feet; x = distance, in feet, downslope from the watershed line; 
K, is a coefficient derived from the Manning formula for turbulent flow, and approxi- 
mately applicable to other types of flow. Its value is: 





1020+/S 
adele (34) 
and the exponent: 
2 2 P 
3M@—«9. — 4 (35) 


in which S is the slope, J the index of turbulence, /, the length of overland flow, and 
n is a roughness factor, of the same type as the roughness factor in the Manning 
formula. The equations for turbulent flow are derived directly from the Manning 
formula and the law of continuity and are rational. The equations for other types 
of flow are closely approximate. The depth 6, as given by these equations is the 
total depth of surface detention, including depression storage. The equations for 
5, fail at points close to the watershed line if, as is often the case, depression storage 
persists to the watershed line. In the equations both for turbulent and other types 
of flow, it is assumed that the velocity varies as ~/S, as for turbulent flow. This 
may not be entirely correct although numerous experiments indicate that, in mixed 
flow, most of the resistance is that due to turbulence. Theoretically, for laminar 
flow the velocity should vary directly as the slope, not as +~/S. 

These equations apply primarily to steady flow. Experiments show that they are, 
however, closely approximate during the early stages of runoff, while surface deten- 
tion is building up to its maximum value. While detailed discussion of the effect 
of the various factors on surface-runoff phenomena cannot be undertaken here, com- 
parison of the equations shows that the velocity of overland flow increases, while the 
depth at a given point x decreases, as the slope increases. Increasing roughness 
of the surface decreases the velocity but increases the depth of surface detention. 

The equations for depth and velocity profiles, in conjunction with that for K;, are 
of fundamental importance in relation to erosional conditions, since they express the 
two factors, 6. and vz, which control the eroding and transporting power of sheet 
flow, in terms of the independent variables which govern surface-runoff phenomena. 
There are six variables: (1) rain intensity, 7; (2) infiltration-capacity, f; (3) length of 
overland flow, /.; (4) slope, S; (5) surface-roughness factor, ; (6) index of turbulence 
or type of overland flow, 7. To apply these equations to erosion and gradational 
problems one must also have Jaws governing the relation of velocity and depth of 
overland flow to the eroding and transporting power of overland flow. 


SURFACE EROSION BY OVERLAND FLOW 


SOIL-EROSION PROCESSES 


There are always two and sometimes three distinct but closely related processes 
involved in surface erosion of the soil: (1) tearing loose of soil material; (2) transport 
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or removal of the eroded material by sheet flow; (3) deposition of the material in 
transport or sedimentation. . If (3) does not occur, the eroded material will be car- 
ried into a stream. 

Every farmer has noticed that the spots most vulnerable to erosion are the steeper 
portions of the hill or valley slopes, neither at the crest nor at the bottom of the hill 
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Ficure 14.—Half profile of a valley slope 


Illustrating soil-erosion processes. 


but intermediate. All soils possess a certain resistivity to erosion, and this resistivity 
may be increased greatly by a vegetal cover, especially a good grass sod. The under- 
lying soil may have a much smaller resistivity to erosion, and, if the surface conditions 
are changed by cultivation or otherwise so as to destroy the surface resistance, erosion 
will begin on land which has not hitherto been subject to erosion. 

Figure 14 shows a half profile of a typical stream valley slope, with the vertical 
scale greatly exaggerated. The line oabc represents the soil-surface profile—flat in 
the region 9, near the crest, steepest in the region ab, about mid-length of the slope, 
and relatively flat at the foot of the slope, in the region bc. The line odef represents 
the surface of sheet or overland flow in an intense rain, the depth of overland flow 
increasing downslope from 0 toward f. In the region oa no erosion occurs throughout 
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a distance x, from the crest of the slope, and this is called the belt of no erosion. Here 
the energy of the sheet or overland flow is not sufficient to overcome the initial re- 
sistance of the soil surface to erosion, even in the most intense storm. In the belt 
ab, mid-length of the slope and where the slope is steepest, active erosion occurs. 

Beginning at a the amount of material carried in suspension by the overland flow is 
proportional to the ordinate between the dotted line ab’c’ and abc. At a it is zero; 
at b it is represented by the vertical intercept 5b’. Beginning at a a given volume 
of water, for example, the water flowing over 1 square foot of soil surface, picks up a 
certain amount of eroded soil matter and carries it in suspension. Passing over the 
next adjacent square foot of area the same water picks up another increment of soil 
matter and holds it in suspension, and so on, the amount of material in suspension 
increasing until at some point 0 the overland flow is fully charged with material in 
suspension and can carry no more. Between the point 6 and the stream channel, no 
material is carried away because any material picked up must be replaced by an 
equal quantity of material deposited from that already in suspension. If the slope 
decreases as shown on the diagram, then the ability of the overland flow to carry 
away imaterial may decrease, in which case deposition of material or sedimentation 
on the surface will occur instead of erosion. 


RESISTANCE TO EROSION 


The physical factors governing soil erosion are: (1) initial resistivity, R;; rain 
intensity, 7; infiltration-capacity, f; velocity and energy of overland flow or eroding 
force, F. The breaking down of the soil structure, tearing the soil apart and lifting 
or rolling soil particles or aggregates, requires the expenditure of energy. Erosion 
can occur at a given location only where the amount of energy expended as frictional 
resistance on the soil surface exceeds the amount of energy required to overcome the 
initial resistance of the soil to erosion. An exception occurs in some cases where the 
soil is churned up inte a semifluid mass by intense rain before surface runoff begins, 
producing high initial erosion rate. Sustained erosion can occur only where the 
condition above described is fulfilled. 

The term “soil” as related to surface erosion includes not only the soil substance 
but also the vegetal cover and the structures—physical and biologic—in the surface 
layers of soil. Soils are of two general classes: (1) indigenous, or those formed by 
weathering of underlying parent rock, either igneous or sedimentary. Such soils 
generally prevail outside of glaciated and loess-covered regions. For some types of 
rock the formation of soil in situ is extremely slow. After a shallow surface layer of 
soil is formed, the formation of additional soil is restricted by the previously formed 
soil cover, but even in full exposure the rate of soil formation from many types of 
consolidated and igneous rocks is so slow that when the soil cover has been removed 
the land becomes worthless. (2) Preformed and transported soils. These consist 
of rock material comminuted by glacial or aeolian action and transported and de- 
posited. Such soil is often a mixture of transported and indigenous soil material and 
includes sedimentary soils deposited on lake or ocean floors and afterward exposed. 
Transported soils, particularly those of giacial origin, are often highly fertile at the 
time they are laid down, as is evidenced by the growth of thrifty forest vegetation 
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within a few years on soils recently exposed by glacial retreat. Erosion of trans- 
ported soils, while a serious menace, is not in general so completely destructive as in 
case of indigenous soils unless the land is gullied and scarified to such an extent as to 
make cultivation impracticable. With equal runoff intensity the resistance of soil 
material to erosion generally increases with the fineness of the soil particles or soil 
texture, the resistance being small for fine uncemented sands but so high for cemented 
hardpan and tough clay that erosion rarely if ever occurs even on bare soil. 

Resistance to erosion is, however, governed more largely by vegetal cover, biologic 
structures, and physical structure of the soil in the surface layers than by soil struc- 
ture. A soil which forms a hard crust on drying may be highly resistant to erosion 
although the same soil when newly cultivated erodes easily. The coherence of soil 
particles and consequently the resistance to erosion is generally increased by the 
presence of colloidal matter, particularly that of vegetal origin. Vegetal cover is the 
most important factor in relation to initial resistance to soil erosion. Its effects on 
the resistivity of the soil to erosion are complex but include: 

(1) Vegetal cover breaks the force of raindrops, thereby reducing the effect of the 
energy of falling rain in breaking down the crumb structure of the soil and packing the 
soil surface. For some soils with little coherence, breaking down of the crumb struc- 
ture by rain impact reduces the soil to a fluid condition, readily susceptible to erosion, 
while for other soils packing of the soi! surface tends to increase the resistance to 
erosion. 

(2) A grass sod operates somewhat like a carpet covering the underlying soil and 
tends strongly to inhibit erosion. 

(3) Fine soil particles adhere to root hairs and plant roots near the soil surface and 
act strongly as a soil binder. In a forest similar effects are produced largely by the 
grass cover but are accentuated by differences in soil structure as between natural 
or undisturbed and cultivated soils and by the presence of an undisturbed humus 
layer near the soil surface. In addition there is often a dense matting of roots of 
trees, herbaceous vegetation, and litter within a forest. Some of the runoff may be 
subsurface runoff and pass through this mat of litter and roots but at so greatly 
reduced velocity as to inhibit erosion. Factors have been devised which stress the 
resistivity of soil material to erosion in terms of the chemical and physical composition 
of the soil. Such factors are, however, inadequate to express the resistivity of a given 
terrain to erosion because of the predominant effect of vegetation and soil structure 
and condition, which are not reflected in indexes of the erodibility of the soil material 
itself. 

The resistivity of a given terrain to surface erosion can be expressed quantitatively 
in terms of the force in pounds per square foot required to institute erosion. This 
quantity can readily be determined on a given soil surface from measurements of the 
distance from the watershed line downslope to a point where erosion begins. 

Nearly all the factors which control resistance of a soil to erosion also controi 
infiltration-capacity of the soil. At a given point on a given slope and with a given 
rain intensity the erosion rate is governed by various factors, one of the most impor- 
tant of which is the infiltration-capacity of the soil. In many instances factors which 
tend to promote a high resistance to erosion also tend to restrict or reduce the in- 
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filtration-capacity and vice versa. Consequently open-textured, coarse, sandy soils, 
such as soils of sand dunes, with little vegetal cover, may never be subject to erosion 
even in the most intense rains although the soil has little resistance to erosion, be- 
cause the infiltratiou-capacity is so high that little or no surface runoff ever occurs. 


ERODING FORCE 


Erosion by aqueous agencies involves three processes: (1) dislodgment or tearing 
loose of soil material and setting it in motion. This is called “entrainment.” (2) 
transport of material by fluid motion. (3) sedimentation or deposition of the trans- 
ported material. 

Let x = distance from divide or watershed line, measured on and along the slope 
(not horizontally); 6. = depth of overland flow at x, in inches; w, = weight per cubic 
foot of water in runoff, including solids in suspension; a = slope angle; » = velocity 
of overland flow at x, in feet per second. The energy expended in frictional re- 
sistance per foot of slope length on a strip 1 foot wide running down the slope, per 
unit of time, for steady flow, will be equal to the energy of the volume of water passing 
over a unit of area per unit of time. This is the product of the weight, fall, and veloc- 
ity, or: 

~~ 
e = m7? sin a. (36) 
The energy equals the force times the distance moved. Hence the force exerted 
parallel with the soil surface per unit of slope length and width is: 
e éz , 
F, = ead laa a. (37) 
Equation (37) is known as the DuBoys formula and is a rational expression which 
may be used as a basis for determining the eroding force per square foot of soil surface. 
F, is the force available to dislodge or tear loose soil material. Sometimes not all 
this force is expended in tearing loose soil material; on a grass-covered surface a large 
portion of it may be expended in frictional resistance on grass blades and stems and 
have little effect in tearing loose soil material unless or until the grass is broken 
over by the impact of the overland flow. However, these factors are best included 
in the resistivity R; of the soil to erosion. 

As shown in connection with surface runoff, for turbulent flow the depth of over- 
land flow at the distance x from the watershed line can be expressed in terms of slope, 
runoff intensity, and surface roughness: 


ao nx\sl5 4 
6, = (<=) . yea ‘ (38) 
The slope S is ordinarily expressed as the tangent of the slope angle a or as tan a. 


Also, for steady overland flow, « = qs, in inches per hour. Substituting these values 
in (38): 





qnx\3/5 1 
6: = BIS SOY ° . 
1020 tan®? a 
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Substituting this value of 6, in (37) gives as the total eroding force at x: 


Wy [ Wanx 36 sina 
i= — ;, 
rs (3) tan®? @ (39) 








CRITICAL LENGTH x,—-BELT OF NO EROSION 


As indicated, erosion will not occur on a slope unless the available eroding force 
exceeds the resistance R; of the soil to erosion. The eroding force increases down- 
slope from the watershed line (Equation 39). The distance from the watershed line 
to the point at which the eroding force becomes equal to the resistance R; is called 
the “critical distance” and is designated x,. Between this point and the watershed 
line no erosion occurs. This strip adjacent to the watershed line, and immune to 
erosion, is designated the “‘belt of no erosion.” An expression for the width of the 
belt of no erosion can readily be obtained from equation (39) by substituting R; for 
F,, making x = x,, and solving the equation for x,. The runoff is free from sediment 
where erosion begins, and w, = 62.4 lbs. per cu. ft. 


12 \s/s 
020{——} = 65.0. 
1020 (2) 65.0 


Substituting this constant in (39) gives: 


65 ( R; \s8 
qn i) me: 


where f(S) is a function expressing the effect of slope on the critical length x, and is 
given by the equation: 


sin a 


tan”? @ ° (41) 


F(S) = 


Numerical values of the slope function f(S) are given on Figure 15. For slopes less 
than 20°, f(S) increases nearly in proportion to the slope. The critical length x, 
varies inversely as the runoff intensity g, in inches per hour, inversely as the rough- 
ness factor m, and directly as the 5/3 power of the resistance R; (equation 40). 

Table 6 gives numerical values of x, for R; = 0.01, 0.05, 0.10, 0.20, and 0.50 lb. 
per square foot, for slope angles of 5°, 10°, and 20°, and for four different runoff in- 
tensities. These are computed for the roughness factor m = 0.10 but can easily be 
applied to other roughness factors, since the value of x, is the same if the product gn 
is the same. 

Renner (1936) observed the percentages of areas having different slope angles 
which were subject to erosion in the Boise River drainage basin, Idaho; his results 
are shown on Figure 16. The extent to which erosion occurred on a given slope 
increased to a maximum on a 40-degree slope and thereafter decreased to zero ap- 
proaching a 90-degree slope angle. 

A comparison of the results of Renner’s observations with the value of the slope 
function f(S) is also given on Figure 15, and the two curves are in close agreement. 
Equation (41) for the slope function is rational in that it is based directly on funda- 
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mental physical laws and principles with respect to small and moderate slopes. For 
such slopes both Renner’s observations and observations by Fletcher and Beutner 
(1941) show that for slopes of less than about 20 degrees the amount of erosion 
generally increases about in proportion to the slope. For steeper slopes the slope 





80 


Horton 
Ea ? 








Erosion ~ Percent 
































49 GO 
Gradient — Percent 


Ficure 15.—Horion slope function for surface erosion 
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sina, 
J(s) = ——,- , where a = slope angle. 
tan-“a 


function f(S) must be considered as empirical, but its validity is confirmed by com- 
parison with Renner’s observations. This function is predicated on uniform turbu- 
lent flow. For very steep slopes such flow cannot occur. 

The critical length x, is the most important factor in relation to the physiographic 
development of drainage basins by erosion processes and also in relation to erosion 
control. The value of x, (Table 6) is highly sensitive to changes in the variables by 
which it is controlled, in particular the resistance R; and the runoff intensity q,.. 
With a newly cultivated bare soil, with R; small, 0.05 Ib. per square foot, for example, 
with a 10-degree slope and a runoff intensity of 1 in. per hour, the width of the belt of 
no erosion would be 35.1 feet, whereas on the same terrain, but with a good, well- 
developed grass sod to protect the soil, and R; increased to 0.5 Ib. per square foot, the 
belt of no erosion would be 1573 feet wide. The width of the belt of no erosion 
varies with the rain intensity, and, consequently, regions with frequent storms of high 
rain intensity are much more subject to erosion, other things equal, than regions 
with lower rain intensities and less frequent heavy storms. Furthermore, in a given 
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TABLE 6.—Critical length x- for various values of R; 



























































Slope angle a, and runoff intensities g,, with roughness factor nm = 0.10: x, = = (ar 3 
qn F 
| @s, inches per hour 
| 0.5 1.0 1.5 | 2.0 
Ri a (degrees) | 
1/qan 
20 10 | 6.67 5.00 
01 5 10.66 i. ae 3.55 | 2.67 
.01 10 | 4.80 | 2.40 | 1.60 | 1,20 
01 20 2.28 | 1.14 | 0.76 | 0.57 
.05 5 | 153.4 | 76.7 | 51.16 | 38.3 
.05 10 70.2 a 23.41 | 17.6 
.05 20 | 32.6 | 16.3 | 10.87 | 8.2 
| | 
| 
.10 5 | 487.6 | 243.8 | 162.61 | 121.9 
.10 10 | 218.4 | 109.2 | 72.84 | 54.6 
.10 20 | 402.8 | 51.4 | 34.28 | 25.7 
| 
.20 5 | 1535.2 | 767.6 512.0 | 383.8 
.20 10 | 689.0 | 344.5 229.8 | 172.2 
.20 20 | 325.0 | 162.5 108.4 | 81.3 
| | | 
.50 5 | 7046.0 | 3523.0 2349.8 | 1762.0 
.50 10 | 3146.0 | 1573.0 1049.2 | 787.0 
.50 20 | 1478.0 | 739.0 4929.0 | 369.0 
f 
w 
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Ficure 16.—Gradient and degree of erosion 
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locality, long intervals may elapse during which no erosion occurs on a given slope, 
than a rain of sufficiently high intensity may cause erosion on a part of the slope. 
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For example, with a moderately good grass-covered slope, with R; = 0.20 and a slope 
angie of 5°, a slope 700 feet in length would not be subject to erosion with runoff 
intensities less than 1 in. per hour, whereas with 2 in. per hour runoff intensity nearly 
the entire lower half of the slope would be subject to erosion. 

Most slopes do not have a uniform gradient from the watershed line to a stream 
but are flattest near the summit, steepest in the middle portion, and again fiat adjoin- 
ing the stream. For such a slope the belt of no erosion will usually comprise all the 
upper, flatter portion. If, for example, the slope length is 2000 feet, g, = 1.5 in. per 
hour, and the mid-portion of the slope has a gradient of 10° and a resistivity of 0.5 
lb. per square foot, erosion will begin 1049 feet from the watershed line. If the 
lower 250 feet of the slope is flatter (its gradient being 5°), then the length of over- 
land flow required to produce erosion with this slope would be 2350 feet. Conse- 
quently, no erosion would occur on the lower or flatter portion of the slope. This 
example shows why erosion is generally confined to the steeper, middle portion of a 
given slope (Fig. 14). 

As another example, on a slope which has been moderately well protected by grass 
cover, with R; = 0.20, slope angle 5°, and with the limiting maximum value of runoff 
intensity for the given terrain 2 in. per hour, the width of the belt of no erosion in the 
maximum storm would be 384 feet. If the slope length was 400 feet, erosion would 
occur only at the foot of the slope and at rare intervals. If the resistivity of the soil 
was reduced, for example, by overgrazing, to 0.1 lb. per square foot or half its former 
value, then the width of the belt of no erosion in a maximum storm wotld be reduced 
to 122 feet, and some erosion would occur on the lower part of the slope, while for a 
runoff intensity of 1 in. per hour the width of the belt of no erosion would be 243 feet. 
Under the changed conditions some erosion would occur in all storms with runoff 
intensities exceeding 0.5 in. per hour. An area having a low resistance to erosion 
and on which erosion occurs over nearly the entire area in the more intense storms 
becomes partially immune to erosion in lighter storms. Between storms of maximum 
intensity, resistivity to erosion may be built up by the growth of grass or trees so 
that when the next succeeding maximum storm occurs the surface resistivity- is 
increased, and the areal extent of erosion greatly diminished. In this manner Nature 
tends to correct the deleterious effects of surface erosion. Another result of im- 
portance is the fact that, on an area subject to erosion only in maximum storms, the 
total amount of erosion over a given time interval—a century, for example—may be 
relatively small, while, if the area is subject to erosion in storms with moderate as 
well as maximum runoff intensities, then because of the much greater frequency of 
storms of lower runoff intensities, the total erosion will be enormously increased. 
Not uncommonly the entire surface of the soil is removed in a century or less. 

Another factor of importance in relation to erosion is that the soil surface, if pro- 
tected by vegetation, has commonly a resistance to erosion many times greater than 
the underlying, unprotected soil. If the surface protection is removed and a maxi- 
mum storm occurs, erosion will then take place at a rate governed by the lower re- 
sistance of the underlying soil. The soil once exposed to direct erosion may then be 
rapidly removed. Soil removed in the belt of erosion may either be carried away or 
deposited farther downslope. The manner in which these combined effects develop 
and control the forms of valley cross sections is considered later. 
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EROSION RATE 


If a factor & be introduced in equation (37) to reduce the erosive force to terms of 
quantity of soil removed, as, for example, the depth in inches of soil materia’ pr. 
moved per hour, then the erosion rate at the point x would be, making F, = ent 


i: 


Re oc 

Cr = ket, p= (42) 
This equation is rational in form and in fact if the rate at which soil material is tom 
loose is proportional to the force available from frictional resistance on the soil surface, 
It relates, however, only to the rate at which soil material can be torn loose and does 
not take into account the ability of overland flow to transport material in suspension. 
Equation (42) is limited in its applicability to cases where the erosion rate is less than 
the transporting power. 


TOTAL EROSION AND EROSION DEPTH 


Beyond the critical distance «, and where the overland flow is not loaded to capac- 
ity with solid matter in transport, the erosion rate at a given point x may be assumed 
to be proportional to the net eroding force. Introducing a proportionality factor 
k. in equation (39) to reduce erosion force to equivalent depth of solid soil material 
removed from the surface per unit of time, and making F; equal the erosion rate ¢, 
and subtracting the value of F; at x, gives: 








= st (2 iz F(S) (a8 — x88), (43) 
The total erosion between x, and «x is: 
E; = / ° €, dx 
or, letting: ‘ J’ 
; 3/5 
“ir © C 
and integrating equation (43): 
E, = B(x8/s — xlsx), (45) 
Substituting the slope length /, for x, 
E, = Bij!’ — 1x3!) (46) 
= Bi,(i3/8 — 3/5), (47) 


The quantity ordinarily measured or measurable in the field is the total erosion per 
storm. Equation (47) can be used to work back from measured total erosion to the 
physical characteristics of the terrain which govern erosion rate. 

The average depth of erosion between x, and x is: 


E; 
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The average depth of erosion is commonly expressed in terms of depth on the 
entire area, not merely the depth on the part of the area where erosion occurs. When 
so expressed the average depth of erosion is: 


as : = BR! — x2!"). (48) 


The coefficient B in this equation is 5/8 of the coefficient of the term containing x 
in equation (43). Consequently, the average erosion depth over a given area is, for 
turbulent flow, 5/8 of the erosion depth for the same time interval at the point x. 

If the value of x, is determined from field observation, together with the average 
erosion depth, the slope length and runoff intensity being known, it becomes possible 
to determine for a given field or area the erosion force F; and the constant k,. The 
latter is: 


__ Erosion depth 
*  Eroding force ° 
These equations form a practical working basis for determining the erosion constants 
R; and k, and for comparing the erosion conditions on a quantitative basis for differ- 
ent areas. 


RELATION OF EROSION TO SLOPE LENGTH 


The average depth of erosion on a given slope increases as the 3/5 power of the 
slope length minus the 3/5 power of the quantity ,, which is constant for a given 
slope and storm (equation 48). For example, the relative erosion rates for different 
slope lengths, with x, = 100 ft. and B taken at unity, are as follows: 


lL = 100 200 500 1000 2000 5000 
Bis = 15.85 24.00 40.3 63.10 95.0 166.2 
Es = 0.00 8.15 24.45 47.25 79.15 140.35 


For a soil with R; and x, each zero, the relative erosion rates would be as shown in the 
second line of the table—i.e., proportional to i3*. The actual relation of erosion 
to slope length is, however, not quite so simple. In equation (48) and in computing 
the figures given abeve, surface-runoff intensity g, has been assumed constant for 
all slope lengths. Other things equal, runoff intensity in a given storm decreases 
somewhat as the slope length /, increases. Also, if erosion rate is determined as an 
average for a year or for several storms, the width of the belt of no erosion will vary 
in different storms. There will in general be more storms producing erosion on the 
middle and lower than on the upper portions of the slope. 

Comparable determinations of soil loss by erosion over a period of 4 to 7 years 
have been made by the U. S. Soil Conservation Service at several stations, using 
slope lengths of 145.2, 72.6, and 36.3 feet. Some of the reported results are shown 
on Figure 17 (Bennett, 1939, p. 152) Differences of soil type, slope and resistivity, 
rainfall, and infiltration-capacity account for differences in soil loss for a given slope 
length at the different stations. Variations of these conditions in different portions 
of the same slope length account for small differences in the relation of slope length to 
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erosion at a given station. in all these experiments the soils were under cultivation, 
producing corn or cotton, and the values of x, were small, particularly in summer, 
The results of these experiments cannot be directly compared with equation (48) 
because x, is unknown, and x, was probably much greater in winter than in summer, 
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if Figure 17.—Relation of erosion to slope length 


From field experiments of U. S. Soil Conservation Service. 


Allowing for differences, the results are entirely consistent with equation (48). As 
shown by equation (48) and also by these experiments, the greater the rain intensity 
: and erosion, the greater, in general, will be the variation of erosion rate with slope 
f i length. Rains of low intensity may produce erosion on scattered patches of soil as a 
: result of local variations of infiltration-capacity and resistivity and with little rela- 
tion to slope length. 


See 


RAIN INTENSITY AND EROSION 


Maximum rain intensities in a given locality generally occur in storms of the 
summer thunderstorm type. For such storms the highest rain intensity usually 
occurs before the middle of the storm and frequently within a few minutes after the 
beginning of rain. Some soils are easily pulverized when excessively dry but possess 
coherence through the operation of capillary force when partially dried after a gradual 
wetting. If an abrupt intense rain occurs on such a noncoherent soil the soil may be 
beaten into a pasty semifluid mass by rain of high intensity before runoff begins and 
before the soil surface becomes protected by surface detention. Such a semifluid 
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mass of soil may accumulate on the surface until it becomes sufficiently fluid, or the 
runoff intensity becomes sufficiently great, to overcome its plastic resistance to flow. 
It is then carried into the stream by surface runo‘f en masse. The combination of 
sudden high rain intensity on previously dried soil of low coherence is not uncommon 
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—— with high initial infiltration-capacity 
FiGurE 18.—Relation of initial infiltration-capacity to erosion 


in summer storms in semiarid regions. In addition the initial infiltration-capacity 
of the dry soil is likely to be abnormally high, and this intensifies the effect of rain 
impact by increasing the time during which the soil is directly exposed without a 
protective cover of surface detention. : 

The combination of the conditions described frequently produces what is referred 
to as a “cloudburst flood.” The term “cloudburst flood” is used because of the 
characteristics of the flood rather than those of the rain which produces it. Meteoro- 
logical conditions are, however, involved. Measurements show that, while such 
floods may carry large volumes of solids, they often carry surprisingly little water as 
runoff, : 

In Figure 18 with a rain pattern opg, conditions for initially high and initially low 
infiltration-capacity are shown by solid and dotted lines, respectively. With a high 
initial infiltration-capacity og there is no surface detention or runoff during the in- 
terval ob during which rain intensity has risen nearly to the maximum. With a 
previously wet and packed soil and low initial infiltration-capacity oj, surface deten- 
tion and runoff begin earlier at a, while the rain intensity is still low. Conditions 
such as those first described occur both on the upland and in stream channels pre- 
viously dry. These conditions in both cases usually produce a cloudburst flood, 
characterized by a wall of turbid water or fluid mixed with debris traveling down the 
Stream channel. 
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The solid material carried along in cloudburst floods is popularly believed to be 
derived from stream banks and channels. This may be true, but most of it may be 
derived from the upland, as is evidenced by the presence of fences, trees, logs, sheep, 
and other objects derived from the upland enmeshed in the semifluid mass. When 


Ficure 19.—Erosion of sodded area initiated by the breaking down of grass cover in intense rains 


the flood wave debouches from a mountain canyon, the water in and behind it may 
escape laterally or by infiltration. The mud flow then slows down, sometimes travel- 
ing so slowly that a person could keep pace with it. Finally the mud stops flowing, 
dries out, and becomes a “fossil cloudburst.” Many such “dry floods” may be 
observed below the mouths of certain canyons in the Wasatch Mountains in Utah. 
Sometimes one or more fossil floods are superposed pick-a-back fashion. An accumu- 
lation of such mud flows may form a debris cone (Horton, 1938; Bailey, 1935; Bailey 
et al., 1934). 

Hydraulic conditions do not permit the occurrence of shallow steady flow on steep 
slopes (Horton, 1938; Jeffreys, 1925); the runoff water is concentrated in a succession 
or train of more or less uniformly spaced waves. These waves concentrate the 
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impact of overland flow on irregularities or obstructions on the surface and greatly 
accentuate surface erosion. While grass or close-growing crops are in general an 
excellent preventive of erosion, they may contribute to active upland erosion in 
cloudburst floods. As long as the grass remains standing it decreases the velocity 
but increases the depth of surface detention, and the resistance to runoff is exerted 
at right angles to the grass stems and has little tendency to pull the sod loose (Fig. 19). 
When a certain depth and velocity of overland flow is attained the grass begins to 
flatten down, like the tipping down of a row of dominoes standing on end by pushing 
the first domino against its neighbor. Then (Fig. 19B) the pull of frictional re- 
sistance is exerted parallel with the surface. A bit of sod projecting more, or which 
is relatively less firmly held, than its neighbors is torn loose by the impact of the 
wave train, and erosion begins. A dense grass sod may be torn up and rolled down 
the slope like a snowball (Fig. 19C). Often several parallel strips of soa are torn off 
from the same slope, each strip a few feet in width. The less resistant underlying 
soil is thus exposed to erosion, and a cycle of erosion may begin on a slope which has 
been immune to erosion for centuries. 

Thus the conditions which initiate soil erosion and govern its rate of occurrence are 
simply and readily expressible in quantitative form in terms of known independent 
variables under some conditions, while under other conditions the factors are so 
complicated that neither the necessary and requisite conditions to cause erosion nor 
the rate of its occurrence can be predicted or expressed accurately in quantitive 
terms. 


TRANSPORTATION AND SEDIMENTATION 


The transportation or carrying forward of eroded material by overland flow or in 
stream channels takes place in various ways: (1) as bed load or material moved 
chiefly by being rolled or pushed along the ground surface or along the bottom of the 
stream channel; (2) material called ‘‘suspended load,” which is held more or less 
continuously in suspension by upward currents due to turbulence; (3) material held 
permanently in suspension by molecular agitation, called the “Brownian movement”; 
(4) solids in chemical solution in the water. The last mode of transport, while it is 
the most important process in connection with ground-water erosion, is relatively 
low in order of importance in connection with surface runoff. 

Bed load comprises materials of sizes ranging from large boulders down through 
cobblestones, shingle, pebbles, gravel, and coarse sand. Mud flows are an important 
process in the transport of material in overland flow. This is, properly speaking, bed 
load. As the terms are ordinarily applied in connection with the dynamics of streams, 
there is no sharp line of demarcation between bed load and suspended load, the former 
term applying to material carried along, on, or near, the solid boundary. Flat 
fragments, such as shingle, are transported by sliding or as bed load. Round parti- 
cles, such as fine gravel, may be transported by a combination of rolling, sliding, and 
jumping. This is also counted as bed load. The process of transport of particles 
by hopping from point to point in semielliptic arcs has been described (Gilbert, 1914) 
as “saltation.” . 

In turbulent flow, eddies thrown off at the solid boundary surfaces have an upward 
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component of velocity. At the same time there is a gradual settling of the water 
between eddies, creating a downward current equal in volume to that of the upwan 
current in the eddies. The upward velocity and the magnitude and frequency of 
eddies increase with the velocity of flow and with the roughness of the boundary 
surface. The eddies are slowly dissipated by viscous resistance as they proceed 
upward. Both experiment and theory indicate that an eddy originates near the 
boundary of the fluid as a vortex ring system consisting of the vortex ring and its 
surrounding field. The fluid comprised in the ring retains its identity throughout 
the life of the ring, as in the case of the familiar smoke ring. The fluid comprising the 
field is continually changing, like the water in a wave. 

Saltation in its simplest form involves the picking up of solid particles by ascending 
eddies. Those lifted and transported by the field are usually carried only a short 
distance; those entrapped near the center of the section of a vortex-ring may be 
carried much farther, until they are thrown out of the ring by centrifugal force. 
These two processes are more or less distinct although both depend on the laws of 
vortex motion. For this reason mathematical analyses of bed load and suspension 
transport without taking vortex motion into account are likely to prove inadequate 
and unsatisfactory. A given particle may be thrown out of one eddy at a certain 
level above the channel bottom, picked up by another, and carried forward, and so 
buffeted about, like a player in a football scrimmage. It may remain a long time 
in suspension, finally reaching the bottom, only to rest for a moment and then embark 
on another wild escapade. 

Various attempts have been made to derive empirical expressions for rate of erosion 
from runoff-plot experiments. In general the conditions of the experiments have 
been such that it cannot be determined certainly whether the results represent the 
rate at which the given surface runoff could erode soil material or the ability of the 
sheet flow to transport such eroded material. 

Most of the work done on sediment transport has been in connection with stream 
channels. Turbulent flow consists of laminar flow on which is superposed the effect 
of the transverse motion of eddies. If the flow is turbulent, then only a minute 
fraction of the energy consumed would be required to provide an equal mean velocity 
of laminar flow. The remaining energy becomes, in effect, latent at the boundary 
by conversion into rotational energy of vortex motion. 

Two principal results follow: (1) The mean velocity is reduced from that for laminar 
to that for turbulent flow; (2) the velocity distribution is changed from that for 
laminar to that for turbulent flow. 

For the usual slight depths of sheet flow the energy actually used in translational 
motion of the fluid is a much larger fraction of the total energy available than for 
types of flow commonly occurring in stream channels. The relative roughness is 
usually much greater for overland flow than for channel flow. Sand particles 0.001 
foot in diameter with overland flow 0.01 foot in depth correspond in relative rough- 
ness to boulders 1 foot in diameter in a stream channel 10 feet in depth. Because 
of these and other differences the extent to which experiments and analyses for chan- 
nels are applicable to sediment transport in sheet or overland flow is an open question. 

Channel flow is ordinarily turbulent except in lakes, while overland flow, except 
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on steeper slopes, usually comprises a mixture of areas on which the flow is turbulent 
and depressions through which the flow may be laminar. 

Much more experimental and analytical work is needed on this problem. How- 
ever, the following facts appear to be well established and suffice for present purposes: 

(1) The transporting power of sheet flow increases with the amount of eddy 
energy due to surface resistance. 

(2) Kinetic energy varies as the square of velocity, and transporting power of sheet 
flow must vary at least as the square—perhaps as some higher power—of the velocity. 

(3) There is a maximum or limiting volume of eroded material which can be trans- 
ported in suspension by a unit volume of overland flow at a given velocity. 


ORIGIN AND DEVELOPMENT OF STREAM SYSTEMS AND THEIR VALLEYS BY 
AQUEOUS EROSION 


RILL CHANNELS AND RILLED SURFACE 


The first step toward the gradation of newly exposed sloping terrain is the develop- 
ment of shallow parallel gullies wherever the length of overland flow is greater than 
the limiting critical distance x,. These are “rill channels,” and a surface covered 
with such channels is a “rilled surface.” Rill channels are usually relatively uni- 
form, closely spaced, and nearly parallel channels of small dimensions which are 
initially developed by sheet erosion on a uniform, sloping, homogeneous surface. 
They are sometimes described as “shoestring gullies,” but the term “gully” as or- 
dinarily understood connotes larger and less regular channels developed by sheet 
erosion at a later stage. A rilled surface presents a striated appearance in plan and a 
finely serrated appearance in cross section. 

Excellent examples of rilled surfaces may be found in newly made road cuts, on the 
slopes of highway and reservoir embankments, spoil banks, and mine dumps. In 
road cuts, water often drains onto the newly made slope from above the edge of the 
cut, in which case the rill channels extend the full length of the slope. If such drain- 
age does not occur then the rill channels invariably begin at a little distance below 
the top edge of the cut. Actually the value of the critical length x, may be very 
small, a few inches to a few feet, on newly exposed steep slopes, and the fact that the 
rills do not extend to the top of the slope would not ordinarily be noticcd. Where a 
rilled surface develops on a newly exposed slope the usual result is the development 
of a deep central master rill or gully, with more or less parallel, shallower, shoestring 
rills, decreasing in depth and frequency, on both sides of the master gully. These 
shoestring rills do not generally survive. The deeper ones close to the master gully 
are absorbed by the master rill by bank caving or are destroyed by the breaking. 
down of the narrow ridges between them. Those more remote are later obliterated 
when lateral slope has developed sufficiently to permit cross flow. 

On some newly exposed lands, with high -infiltration-capacity and high resistivity 
to erosion, the length of slope from the major divide to the downslope edge of the 
area may never exceed x,. Under these conditions a rilled surface may not develop. 
This condition often occurs on sand-dune areas and in some glaciated areas with deep 
permeable soils, especially where grass or other vegetal cover develops soon after the 
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disappearance of glacial ice. In the latter case x, on the newly formed surface may 
exceed the values of /, pertaining to the drainage of melt-water from the ice sheet, g 
rilled surface may not develop, and the topography will remain much the same as 
when the ice disappeared, except that gradation by solution may take place. Jy 
desert regions, with suitable relations between the rain intensity, infiltration-capac. 
ity, surface resistivity, and the slope, a rilled surface may develop with little or no 
cross-grading, so that surface gradation may never extend beyond the rill stage. 


ORIGIN OF RILL CHANNELS 


“Sheet erosion” implies the formation of either a rilled or gullied surface. From 
the discussion thus far it would appear that overland flow downstream from the 
critical point x, on a smooth uniform surface should remove a uniform layer of soil 
instead of producing a rilled surface. 

The question may fairly be asked: Why does a drainage basin contain a stream 
system? Surface runoff starts at the watershed line as true sheet flow, without 
channels. Even below the critical distance x, it should apparently continue as such 
sheet flow combined with sheet erosion. Why, then, do rill channels develop? The 
answer is that channels start to develop where there is an accidental concentration 
of sheet flow. Accidental variations of configuration may provide the requisite 
initial conditions where a local area has a lateral slope joining a longitudinal slope 
or where two lateral slopes join and form a trough. 

Most cases of active erosion observed at present represent conditions where there 
is or has been a protective vegetal cover and the initial resistivity of the soil surface 
is greater than that of the immediately underlying subsoil. 

Consider a point upslope from x,. If, as a result of change in cover conditions, 
either the resistance R; or the infiltration-capacity is reduced, the point «, may move 
upslope from the given point, which will then be susceptible to erosion. When the 
remaining protective cover is broken through at a given point, a channel or gully will 
form which will proceed rapidly upslope, chiefly by headward erosion, because of the 
lower resistivity of the underlying soil. 

This, however, is not the mode of origin of rill channels, which, it must be pre- 
sumed, often form on new ter:ain without vegetal cover and with a value of R; 
sensibly the same at and to some depth below the soil surface. Slight accidental 
variations of topography may produce a sag in which the depth of sheet flow is a 
maximum at the point a (Fig. 20), the line 6b’ representing the water surface at 
maximum runoff intensity. As a result of the greater depth at a, erosion will be 
most rapid at that point, and increased channel capacity will be provided at a, and 

‘part of the water which originally flowed in shallower depths on the adjacent area 
will be diverted into this enlarged channel. This may accelerate the process until 
the entire flow is concentrated in the rill channel (Fig. 20). This does not involve 
headward erosion in the ordinary sense. However, when a rill channel has once 
formed, sheet flow coming down the slope upstream from the head of the rill will be 
deflected toward and diverted into the rill channel, thus providing a means of rapid 


headward extension of the rill. 
This process of rill formation can often be observed on a cultivated slope during 
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heavy runoff. The size and spacing of the rill channels vary with the slope, runoft 
intensity, and length of overland flow, ranging from a few inches apart on a culti- 
yated slope to many feet or yards apart on long slopes with low runoff intensity and 
higher erosive resistance. In some areas in abandoned lake beds or exposed coastal 























FiGurE 20.—Successive stages of rill-channel development 


belts in arid regions stream development has never progressed beyond the rill stage. 
Such an example is given on the Moon Mountain, Arizona-California, quadrangle, 
U.S. G. S. topographic map (Fig. 21). Later stages of stream-channel development 
belong to the domain of channel dynamics and involve velocity distribution, silt 
equilibrium, and other factors which cannot be considered fully here. 

The ultimate dimensions of a stream channel are, as indicated by Playfair’s law, 
such that it is adapted to the area which it drains. Stream channels tend to acquire 
ultimate dimensions such as to carry all or most of the flood waters of the stream. 
This is largely because most surface erosion and channel erosion occur during floods. 


CROSS-GRADING AND MICROPIRACY 


A system of parallel shoestring gullies is transformed to a dendritic drainage net 
as the result of the tendency of the water to flow along the resultant slope lines and 
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Ficure 21.—Portion of Moon Mountain, Ariz.-Calif., quadrangle, U. S.G. S. 


is a direct consequence of the overtopping and breaking down of intermediate ridges 
between gullies by overland flow during heavier storms. 

The deepest and widest rill develops where the net length /, — x, in which erosion 
can occur is greatest. If x, varies, this may not occur where the total length /, of 
overland flow is greatest. The longest, deepest, and strongest rill channel will be 
called the “master” rill. Owing to smaller values of /, — x, proceeding away from 
the master rill on each side, the rills will be shallower, or, considering two adjacent 
rills, the bottom of the one farther from the master rill will be higher. 
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When a storm occurs exceeding in intensity preceding storms on the newly exposed 
areas, the divide between two rills may be broken down at its weakest point by (1) 
caving in of the divide between two rills, diverting the higher into the lower rill; (2) 
erosion of the divide by the deeper or lower rill, thus diverting the higher rill; (3) 
overtopping of the divide at the low point by the higher rill, again diverting it into 
the lower rill. This breaking down of divides between adjacent rill channels and 
diverting the higher into the lower rills is described as micropiracy. Micropiracy 
much resembles stream capture by lateral corrasion, but micropiracy results chiefly 
from water overtopping a low spot in the narrow ridge between two rills. Micro- 
piracy obliterates the original system of rills and their intermediate ridges on a 
uniform newly exposed surface. The process of erosion, in the course of development 
of a stream system and its accompanying valleys destroys most of the record of their 
origin. Ultimately the original slope parallel with the stream is replaced on each 
side of the stream by a new slope deflected toward the stream. This process is 
described as ‘‘cross-grading.”’ 

The initiation of cross-grading is illustrated on Figure 22, which shows a plan of a 
small area of newly exposed land, aa’bb’. The line cc’ marks the downslope limit 


of the belt of no erosion, aa’cc’.. The critical distance x, is assumed to vary with . 


slope, infiltration-capacity, and initial resistance to erosion. Rills develop down- 
slope from cc’, and their development is followed by cross-grading, as shown on the 
cross sections taken on the line dd’ and numbered 1 to 4, inclusive. The line of 
resultant slope in each case is in the direction shown by the arrows, and the rills 
increase in depth and degree of gradation at a given time proceeding away from the 
lateral boundaries toward the initial or master rill. In section 2 (Fig. 22) the rilled 
surface has developed, but flat “lands” still persist between rills, and the resultant 
slope is still parallel with the original slope. In section 3 some rills have combined 
by cross-grading, creating slight cross slopes, but the overland flow is still carried 
chiefly by rills parallel with the original slope. In section 4, with increased cross 
slopes and perhaps a heavier storm, active cross-grading has taken place, especially 
mid-length of and in the lower portion of the original slope; the direction of overland 
flow is no longer parallel with the original slope, but overland flow takes place partly 
in the rills and partly across the intervening ridges, somewhat as shown on a larger 
scale on Figure 23, and in detail for a single pair of rills on Figure 24. This process 
can sometimes be observed in heavy storms on lands cultivated nearly but not quite 
parallel with the contours, the tillage marks corresponding to the original rills above 
described. 

In Figure 22, section 5, the original serrated rilled surface has been obliterated by 
cross-grading and is replaced by an irregularly roughened surface on which a new 
tilled surface tends to develop, with flow lines parallel with the resultant slope. This 
represents the end point of the first stage of valley gradation and stream develop- 
ment. At this stage there exists only one stream in the area, and this follows the 
course of the initial or master rill. This idealized picture of the cross-grading of a 
tilled surface is based in part on field observations of eroding slopes and in part on 
the observed manner of erosion of experimental plots, using artificial rainfall. Under 
natural conditions the results are rarely so uniform as those shown on Figure 22. 

A break across a rill divide may result from numerous causes, such as a rock or 
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FIGURE 22.—Development of a valley by cross-grading 
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Ficure 23.—Successive stages of rill obliteration 


obstruction in the path of the rill, causing back-water upstream therefrom, or the : 
caving in of the bank of the rill, thus obstructing its flow and producing a side outlet. | 
A very common cause is an accidentally low divide between two adjacent rills. In 
most cases the channel of the diverted rill will be obstructed just downstream from f 
the point of diversion. If the break supplies a free outlet from the diverted rill, 
then immediately downstream from the break, as at x (Fig. 24, B), there will be little 
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or no flow, and erosion will cease, while erosion just above the break will continue. 
The channel above x will quickly become deeper than the abandoned rill channel 
below x. There will be thrust, as indicated by the downslope arrows on Figure 24, 
at the point of diversion, resulting in a tendency toward the formation of a rounded 
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Figure 24.—Development of a stream from a rill system by cross-grading 


(A) Showing angular plan of initial stream; (B) Development of bends by thrust at angles of an initial stream. 


bend by impact erosion. Similar thrust across the slope will occur at the second 
angle of each cross rill diversion, with similar results. 

Micropiracy tends to give the resulting stream at first a more or less angular 
course (Fig. 24) which, as the angular bends are rounded by erosion, will finally 
develop a more or less tortuous stream course following generally the line of resultant 
slope. On steeper slopes the resulting stream tends strongly to maintain a straight 
course and has eroding power sufficient to do so. On flatter slopes, with greatly 
reduced eroding power, centrifugal force around the initial angular stream bends 
tends to enlarge their radii. Weak bends merge with stronger ones until ultimately 
a system of stream meanders is developed. Since the master rill either initially or 
by cross-grading, on a given slope, ultimately becomes a permanent stream, it appears 
that conditions favoring the formation of stream bends on flatter slopes are inherent 
with the origin of the streams. 

The elimination of rill systems on gentle slopes is very different from the develop- 
ment of mature meander belts such as those observed in natural streams. Even 
where stream bends originate in this way the evidence of the intermediate steps is 
eliminated in the course of the process. 





tov 
wil 
alo 
cat 
oul 
thi 
nec 
anc 
wil 


to | 
tow 
lon; 
dev 
fort 


Act 
slor 
incr 


furt 
5a 


and 
side 





nnel 
24, 
ded 








ORIGIN AND DEVELOPMENT OF STREAM SYSTEMS 339 


This is not the sole explanation of the origin of stream bends, for individual bends 
arise from other and accidental causes. The fact that bends occur generally in 
rather definite systems on flatter stream slopes and in material that is almost per- 
fectly homogeneous seems to require something more than purely fortuitous causes 
to explain their origin. It is at least significant to find that a stream in its initial 
stage, as it merges from a rill system, provides the necessary conditions for the 
development of stream bends on flat slopes. 


HYDROPHYSICAL BASIS OF GEOMETRIC-SERIES LAWS OF STREAM NUMBERS 
AND STREAM LENGTHS 


General statement.—A conventionalized illustration is given of the main steps 
involved in the development of a drainage net, showing the hydrophysical basis of 
the geometric-series laws of stream numbers and stream lengths. A square area 
with its diagonal parallel with the direction of slope will be assumed. This roughly 
approximates in form a typical ovoid drainage basin. 

First stage —On Figure 25, oabc is a uniform surface sloping toward 0. The maxi- 
mum length of overland flow is the diagonal length J, = ob. A storm produces 


surface runoff of intensity sufficient to reduce the critical distance x, to some value . 


bg < J,. Sheet erosion can then occur over the area ofgh. The runoff intensity 
q, in cubic feet per second per unit width, will increase proceeding from a and ¢ 
toward the center line nearly in proportion to the length of overland flow. There 
will be no erosion at f or 4. The erosion intensity will increase from these points 
along the line fz toward the center line. The maximum intensity of surface runoff 
cannot occur until surface detention is built up to a point where the inflow to and 
outflow from surface storage are equal. Since the unit runoff intensity q,, other 
things equal, increases nearly as the length of overland flow, the critical intensity 
necessary to induce erosion will be exceeded first along the center line of the area 
and then progressively later proceeding toward f and # (Fig. 25,a). Sheet erosion 
will begin first along the diagonal line og, spreading more or less rapidly toward the 
points fand #. Ina typical storm the duration of surface-runoff intensity adequate 
to produce erosion will also be greatest near the center line of the area, decreasing 
toward the sides. Consequently, erosion will begin first, be most intense, and last 
longer near the center line of the area. A series of more or less parallel gullies will 
develop (Fig. 25, a), decreasing in depth and frequency, proceeding away from og, 
forming a rilled surface. 

In Figure 25a the lines of overland flow are parallel with the central or master rill. 
Actually as the valley develops, the lines of overland flow will follow the resultant 
slope lines converging toward the central rill first at acute angles and, as the slope 
increases, at increasingly larger angles, approaching right angles as a limiting case. 

This is the first stage of stream and valley development. It may later produce 
further upslope erosion in more intense storms. At the end of the first stage (Fig. 
25a), there is a single central stream, with a V-shaped valley. 

Second stage—The development of the first stream (Fig. 25a) has divided the area 
and reduced the maximum length of overland flow on the remaining areas on each 
side to 3 /, or to half the initial maximum length of overland flow. Lateral tribu- 
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taries have not yet developed because (1) Lateral flow cannot occur unless there is a 
lateral component of slope, and (2) until the lateral slope extends far enough from 
the main stream so that in a given storm the value of the critical length x, is reduced 












































FiGuRE 25.—Development of a drainage net in a stream basin 
(Schematic). 


to less than the width of the lateral slope, or in case of complete development of the 
initial valley, the value x, must be reduced to less than $ /,. 

The second stage involves the development of lateral tributaries and their valleys. 
During a storm of sufficient intensity to meet the prerequisites described, a pair of 
master rills or lateral tributaries with accompanying shoestring gullies will develop, 
one from each of the equal areas on opposite sides of stream 1 (Fig. 25b). These 
lateral streams are designated “2,” and their development and that of their lateral 
slopes will follow the same course as in the case of stream 1. At the end of the second 
stage the area will appear as shown by Figure 25b. In the meantime, further head- 
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ward development or downward gradation, or both, of the valley of stream 1 may 
have occurred. 

Subsequent stages.—At the end of the second stage the maximum remaining length 
of overland flow in the area has been again halved or reduced to  /,, and a still more 
intense flood is required to reduce the critical length x, to a value less than this and 
permit the development of a third group of lateral tributaries. When such a flood 
occurs, each of the No. 2 tributaries will develop a pair of lateral tributaries, desig- 
nated “3” on Figure 25c. 

The development of lateral slopes adjacent to a given stream brings in overland 
flow from additional areas, increasing downstream, and accelerates grading of the 
main stream and its immediate valley. The stream system at the end of the fourth 
stage is shown by Figure 25d. The number of the stage of stream development 
corresponds to the order of the main stream. The main stream is of the 1st order 
at the end of the 1st stage, and of the 4th order at the end of the 4th stage. 

Development of lateral tributaries and the manner in which they develop is the 
direct consequence of (1) the existence of a critical length x, of overland flow required 
to institute erosion; (2) the operation of cross-grading. Drainage patterns, while 


invariably following the two fundamental geometric-series laws as to stream length 


and stream numbers, can still develop in an infinite variety of ways. 

Two questions naturally arise: (1) What would happen if the newly exposed area 
was a continuous belt along the coast line, with no lateral boundaries? (2) Why 
and how are the boundaries of drainage basins developed? This case will be con- 


sidered later. 
ADVENTITIOUS STREAMS 


Differences occur between the development of streams under natural conditions 
and those assumed in the example because: 

(1) The drainage area is usually not rectangular but ovoid. 

(2) Newly formed tributaries follow in general the resultant slope of the cross- 
graded areas on which they develop and hence enter the parent stream at more or less 
acute angles; the steeper the slope, the more acute is the angle of stream entrance. 
Tributaries thus tend to be longer than if they entered at right angles, and the stream- 
length ratio is consequently usually greater than 2.0. 

(3) There are nearly always variations—sometimes large variations—of infiltra- 
tion-capacity in different parts of the area. 

(4) There are also variations—sometimes extreme—in the initial resistance of the 
terrain to erosion, as, for example, where part of the area is in consolidated and part 
in unconsolidated material, or part covered with vegetation and part bare. As a 
result of these departures from hypothetical conditions, the following results often 
occur: 

In certain parts of the area the length of overland flow to the parent stream or its 
larger tributaries may be less relative to x, than on the remaining areas tributary 
to the last group of streams developed. Then some 1st, and perhaps 2nd, order 
tributaries develop, entering the main stream or larger streams directly and not 
through higher-order tributaries. These streams, which result from accidental 


Paar ref els eet mer omeeepran Aleareereeh ahem aren bres optenpetosioneresieperecernereeeoaeneeteceere te tee 


Reacsbart Slanederaaaenee 


RES 








342 R. E. HORTON-—EROSIONAL DEVELOPMENT OF STREAMS 


variations of conditions within the area, may appropriately be designated “‘adventj. 
tious” streams. The development of adventitious streams increases the number of 
streams of lower orders and tends to make the bifurcation ratio greater than 2.0, as 
it usually is for natural streams. Lateral slope also increases the length of tributaries 
and makes the stream-length ratio also greater than 2.0. An increase of 1; is also 
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FiGureE 26.—Beginning of erosion on newly exposed land 


produced because streams do not extend to the watershed line although they may 
extend by headward erosion to a distance less than the critical minimum value of x, 
from the watershed line. Adventitious streams usually increase the bifurcation 
ratio more than the stream-length ratio is increased by the conditions described, 
. %, oe ‘ 
with the result that the ratio — is in general a fraction, and the total length of streams 
6 
of a given order is not constant but decreases proceeding from the lowest to the 
highest stream orders. Adventitious streams do not in general develop simulta- 
neously with larger streams in the basin but are developed later as the development 


of the stream system approaches maturity. 
STREAM DEVELOPMENT WITH PROGRESSIVELY INCREASING LAND-EXPOSURE COMPETITION 


For illustration the exposure of coast marginal lands will be assumed to be nearly a 
uniform homogeneous sloping plane, extending from a divide line aa’ (Fig. 26) to the 
new coast line cc’. It is assumed that the soil surface in the newly exposed belt 
aa’cc’ is initially bare and has a certain infiltration-capacity f and a surface resistance 
to erosion R, such that the critical length x, required to permit surface erosion to 
occur in the most intense rain is as shown on the diagram. The dashed line bd’ 
is at a distance x, from the watershed line aa’. As long as the coast line is within the 
belt aa’bb’, no streams will develop, runoff will be in the form of direct sheet flow to 
the new coast line, and no erosion takes place. There will be irregularities in the 
watershed line aa’ and in the coast line cc’, and, when the length of overland flow 
exceeds x, at some point d, erosion will begin at that point. When the coast line has 
reached the position cc’ there will be a small area, as outlined by a dashed line, within 
which /, > «., and within this area sheet flow will produce erosion and a series of rill 
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channels parallel with the direction of the initial slope surface. The first rill channel 
will be at dd’, where the length of overland flow first exceeds the critical length x. 


As the coast line recedes the belt in which erosion can occur will increase laterally 
and longitudinally, and the system of rill channels will be extended correspondingly 
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FicureE 27.—Development of first pair of tributaries on new stream system 
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in both directions. The rill at dd’ was first formed and has been longest subject to 
erosion and will become the master rill. Cross-grading will take place, producing 
new components of slope toward the rill dd’ and obliterating the original rill channels. 
New rill channels will develop following the new direction of slope, on each side of the 
original stream dd’ (Fig. 27). In general the lengths of these new rill channels will 
increase proceeding down the slope from the line cc’. At some point o (Fig. 27) 
a new rill channel will have a greater length og and greater runoff than rills between 
oandd. It will have developed earlier than rill channels entering the parent stream 
dd’ between o and d’. It therefore has greater runoff and a longer duration of runoff 
in which to cut its channel than rills formed farther down the slope. Such rill chan- 
nels will survive as a tributary stream. Such a rill channel occurs on each side of the 
parent stream in the vicinity of 0, ai cross-grading toward these tributary streams 
will also occur. Cross-grading of the areas adjacent to these two tributaries will 
produce cross-graded slopes on either side of each tributary (Fig. 28), until there is 
again a location on each of these areas favorable for the development of tributaries, 
and new tributaries will develop, usual!y one on each of the two preceding tributaries, 
as at m, n, and p (Fig. 28). This process will continue until finally there is no land 
surface above the mouths of the original tributaries where the length of overland 
flow exceeds the critical length «x. 
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F If the coast line recedes farther (Fig. 29), the area upslope from 00’ on the right 
hand side of the stream dd’ is tributary to the stream og. The original rill channels 
parallel with dd’ upslope from oo’ have been obliterated, and the runoff from the area 
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Ficure 28.—Lines of flow after cross-grading of first pair of tributary areas 
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oo’o”’ now enters stream og. As the coast line recedes a new system of rill channels 
parallel with dd’ develops downslope from oo’, and the area oo0’d’d’’ will become cross- 
graded toward dd’. When the length of overland flow within the area oo’d’d” be- 
comes sufficiently great at some point g’, a new tributary q’r will develop along the 
line of the resultant slope, and its basin will in turn be developed by cross-grading. 
There must be a certain minimum space or intercept between tributaries of the main 
stream to provide adequate length of overland flow to permit a lower tributary to 
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‘ 
the main stream to develop. Furthermore, the tributary g’r, having developed 
much later than the tributary og, will extend its drainage area laterally more slowly 
than the latter, with the result that the drainage basin will tend to have an ovoid 


outline (Fig. 29). 
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Figure 29.—Development of lower pairs of main tributaries 


Another stream may also develop at zz’ in the same manner as the stream at dd’. 
The development of this stream may have begun either a little earlier or a little later 
than the stream dd’, and the final location of the lateral divide between the two 
drainage basins will be determined by the conditions of competition. The older 
stream will absorb the greater part of the area between the two streams. A mar- 
ginal area of direct drainage d’d’’z’ is left between the two major drainage basins. 
If the length of overland flow here becomes sufficient, an intermediate subordinate 
stream will develop. 

The appearance of the final stream systems in the two drainage basins will be 
somewhat as shown by Figure 30. 

Two major factors control the development not only of the drainage basin of a 
given stream but the systems of drainage basins tributary to a new coast line: 

(1) Streams develop successively at points where the length of overland flow be- 
comes greater than the critical length -. 

(2) Competition results in the survival of those streams which have the earliest 
start or had the greatest length of overland flow, or both, and which are therefore 
able to absorb their competitors by cross-grading. 
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END POINT OF STREAM DEVELOPMENT 


Stream development on a newly exposed slope continues until the greatest remain- 
ing length of overland flow is less than the critical distance x, required to institute 
erosion. 


£. 




















a b 
Ficure 30.—Final development of two adjacent drainage basins on newly exposed land 


At a certain stage of gradation (Fig. 31) the stream oa has developed with a drain- 
age basin ecd. Before cross-grading of this area the critical length x. is, for example, 
equal to that shown by the line mm’ on the insert, and this is less than oa. After 
cross-grading of the area ocd this criticai length is somewhat reduced by increased 
resultant slope and is now mn. The greatest lengths of overland flow on the areas 
oca and oad are now along the slope lines de and ce, but these are both less than mn. 
Hence no additional streams will develop in the area ocd. 

The upper ends of the streams in a drainage basin will extend at least to the dis- 
tance x, from their watershed line, measured in the direction of slope. They may be 
extended closer to the watershed line by headward erosion, under suitable conditions. 

For streams to be perennial at their sources there must be ground-water flow at the 
head of the stream channel. In regions where there is a permanent ground-water 
horizon under the drainage basin the most common condition is that the stream is 
intermittent for a distance downstream from the point where its channel begins. 
Figure 32 shows the profile at the head of a stream. The watershed line is at a, and 
a definite channel begins at b. There is a water table underneath the headwater 
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FiGure 31.—End point of stream development 
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Ficure 32.—End point of a definite stream channel 
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belt of no erosion ab, the maximum ground-water table is at cc’, and the minimum 
at dd’. Between cc’ and dd’ the stream is intermittent. At c’ part of the infiltration 
on the upper drainage area enters the stream. At times of maximum surface runoff 
the ground-water flow may represent a considerable fraction of the total flow. If, 
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Ficure 33.—Drainage basin of Pennypack Creek 


Above Valley Falls, Pa., showing subareas from which surface runoff is derived. 


for example, the ground-water flow is one fourth of the total flow at c’, then, if the 
channel extended a little farther upslope to e, the maximum runoff would be reduced 
one fourth by elimination of ground water. There is therefore an abrupt and some- 
times considerable change in the total runoff at about the point where the maximum 
level of the water table intersects the stream channel. Surface runoff plus ground- 
water flow can generally extend the channel upstream farther by headward erosion 
than could surface runoff alone. Hence the channel usually ends near the point 
where ground-water flow is no longer effective. Ground-water flow at c’ is inter- 
mittent, but it usually continues much longer than surface runoff and by maintaining 
the soil at the head of the stream channel moist and soft it promotes extension of the 
channel by headward erosion and bank caving. 

The final results of stream development under natural conditions are illustrated 
by Figure 33. Some of the streams in the lower part of the basin are clearly adventi- 
tious. There are several drainage basins, such as A and B, where tributaries have 
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developed only on one side of the parent stream, leaving, in this case, an isolated 
plateau in the interfluve area, although the drainage development of the basin is 
evidently mature. 


STREAM-ENTRANCE ANGLES 


From geometrical considerations the following equation has been obtained for the 
entrance angle between a tributary and the higher-order stream which it enters 
(Horton, 1932)*: 


where z, is the entrance angle between the two streams; 5, is the channel slope of the 
parent or receiving stream; 5s, is the ground slope or resultant slope, which is here 
assumed to be the same as the slope of the tributary stream. 

Values of the entrance angle computed by this equation for different values of the 
ratio S-/s, are as follows: 


Se/So = 9.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 
ie om BES SES. AES? 37.0?) GOO? -66.2°, 7253": | 2838": , SR? 


As shown by Table 4, stream slopes are always less than the adjacent ground 
slope, and tributaries should enter the confluent stream at acute angles when the 
slopes of the channels of the tributary and confluent streams are nearly the same. 
The equation takes on the indeterminate form 0/0 if the two slopes s, and s, are equal. 
This means that the two streams will be parallel and will not join. Three cases will 
be considered for purposes of illustration. 

CasE 1—FLat STREAMS DEVELOPED ON A FLAT AREA: When the parent stream 
has developed and cross-grading has proceeded to a point where a pair of tributaries 
develop, the parent stream will in general have cut into the initial surface to some 
depth, and its stream slope in the vicinity of the debouchure of the tributaries will-be 
materially less steep than the original slope, while the slopes of the tributaries as 
they approach the parent stream will be materially steeper than the original slope. 
As a consequence, instead of the ratio s./s, being close to unity, this ratio will seldom 
have a value greater than 1/2 or 1/3, and the tributaries will not enter the main 
stream at acute angles, as would be the case if s, and s, were nearly equal, but will 
more generally enter the parent stream at angles of 60° to 80°. On extremely flat 
surfaces in humid regions a swampy condition often prevails, and stream-entrance 
angles are but little subject to control by erosion conditions. On semiarid plains 
where little erosion occurs, acute entrance angles of tributaries to the parent stream 
may sometimes be observed. 

CasE 2—Fiat VALLEY SLOPE WITH MODERATE TO STEEP ADJACENT GROUND 
SLOPE: Under these conditions the ratio s./s, is nearly always low, and the stream- 
entrance angles to the main or parent stream are commonly 60° or greater. As the 





‘Derivation of this equation is given correctly in the reference cited. Interpretation of the equation as there given 
is incomplete and not wholly correct. 
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stream system develops, the slope of the main stream steepens proceeding upstream, 
and the lateral ground slopes also steepen proceeding upstream. The ratio 5,/s, 
may remain sensibly constant, or it may either increase or decrease. Most commonly 
it decreases to some extent. Quite generally the entrance angles of tributaries to the 
main or initial stream are quite uniform and range from 60° upward, decreasing some- 
what upstream. 

Case 3—TRIBUTARIES ON A STEEP SLOPE: Tributaries developed on the same slope 
generally run nearly parallel, and if the main valley is relatively flat they will enter 
the parent stream at an angle of 90°, representing a limiting condition which is ap. 
proached but not often attained. Tributaries developed on the same lateral slope 
may of course join and are especially likely to join where drainage development is 
incipient, as on steep, rocky slopes and in semiarid regions where tributary develop- 
ment has been arrested at the end of the rill stage. Parallel tributaries which join 
on a steep slope under these conditions commonly have an acute angle of juncture. 
In this case the ratio s./s, is close to unity. 


DRAINAGE PATTERNS 


Much has been written regarding the forms of drainage patterns. They are usually 
classified as dendritic (treelike), rectangular or trellised, radial, and centripetal. 
The terms radial and centripetal commonly refer to the arrangement of a group of 
drainage patterns originating at or converging to a common point and do not refer 
in general to the pattern in an individual drainage basin. All drainage patterns of 
individual drainage basins are treelike, but different patterns resemble the branch- 
ings of different kinds of trees and range from those with branches entering the 
parent stream nearly at right angles, to those with tributaries nearly parallel and 
entering their parent streams at small angles. The form of the drainage pattern 
depends to a large extent on the relation of the slope of the parent stream to the 
resultant ground slope after cross-grading. If this ratio increases with successive 
cross-gradings, stream-entrance angles of successive tributaries are somewhat more 
acute for successively lower-order streams, affording the most usual type of dendritic 
drainage pattern. 

On a relatively flat surface the directions of resultant overland flow after the 
first cross-grading are nearly at right angles to the initial stream, and the second 
series of streams developed enter the parent stream nearly at right angles. Cross- 
grading of the areas tributary to these streams produces but a slight change in the 
slope ratio s./s,,so that the next order of streams also enters the parent streams 
more or less nearly at right angles. In this way a rectangular drainage pattern is 
developed. 

If, on a steep, sloping, original surface, the headwater divide forms roughly an arc 
of a circle, then the first two tributaries developed will enter the parent stream from 
opposite sides at nearly the same point (Fig. 34). These streams will develop long 
tributaries nearly parallel with the initial stream, giving rise to a centripetal drainage 
pattern (Fig. 34). 

On flat slopes each successive cross-grading of a given subarea changes the direction 
of the next stream to develop on the area through an angle approaching 90° as a 
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| Ficure 34.—Centripetal drainage pattern 
Payne Creek, Ga., Mulky Gap quad., U. S. G.$.—T. V. A. 
limit and changes the direction of overland flow through a corresponding angle. 
The direction of resultant cross-graded slope at the end of a given stage becomes the 
direction of the stream of the next succeeding stage. The directions of streams and 
of resultant slopes will change through nearly a right angle with each successive 
stage of stream development and cross-grading, and the directions of streams and of 
resultant slopes tend generally to be the same in any two stages of stream develop- 
ment which are either both even numbered or both odd numbered. 
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ASYMMETRICAL DRAINAGE PATTERNS 


Because newly developed tributaries enter their parent streams at acute angles, 
they divide their tributary areas into two parts such that the remaining upslope 
tributary areas are larger than those on the downslope side, using the terms “up. 
slope” and “downslope” with reference to the two sides of the tributary. Because 
of inequality of area, width, and slope on the two sides of a tributary, the next lower 
order of tributaries may develop with two or three tributaries on the upslope side and 
fewer or none on the downslope side, a common phenomenon, particularly in moun- 
tain areas. Since the average elevation of the upslope area is greater than that of 
the downslope area, this phenomenon is sometimes attributed to increase of rainfall 
with elevation. It may occur, however, as the result of differences of tributary area 
and length of overland flow on the upslope and downslope sides of the parent streams, 
independently of variation of rainfall or runoff on the drainage basin. Burch Creek 
and Reels Creek drainage basins (Utica, New York, quad., U. S. Geological Survey) 
afford examples of asymmetrical drainage-basins. 


PERCHED OR SIDEHILL STREAMS 


In general, streams follow the bottoms of the valleys in which they are located. 
Small—usually 1st order—streams are occasionally perched precariously on the side 
slopes of graded valleys of higher-order streams. The course of such a stream is 
often more nearly parallel with the antecedent slope than with the cross-graded slope. 
At the foot of the slope the stream often turns abruptly and debouches into the parent 
stream at nearly a right angle (Fig. 35). Evidently gradation of the valley of the 
parent stream cd reached the stage shown in the figure before the slope became steep 
enough to reduce the critical length x, below the maximum length /, of overland flow 
on the right-hand side, and J, became greater than x, only when gradation of the 
valley slope had reached the end point. A weak stream, ad, then developed bv 
micropiracy and cross-grading, but owing to some local cause, such as increased 
resistivity of the soil to erosion at increased depth below the original surface, this 
stream was unable to develop a valley of its own by further cross-grading and so 
remained high above the parent stream on the antecedent rilled surface, until, with 
increasing volume and slope, it turned nearly a right angle as it entered the parent 


stream. 
REJUVENATED STREAMS; EPICYCLES OF EROSION 


In the preceding sections it has been assumed that: (1) Uplift or exposure of new 
terrain took place continuously though not necessarily at a uniform rate, the region 
finally becoming stable; (2) the initial resistance R; of the soil surface to erosion re- 
mained constant. The effect of subsequent further elevation or subsequent sub- 


. sidence of an area on which a stream system has already developed has been exten- 


sively discussed in connection with the Davis erosion cycle (Wooldridge and Morgan, 
1937) and will not be considered further here. Before leaving the general subject of 
stream development and valley gradation consideration will be given to the effect 
of (1) differences between surface and subsurface resistivity to erosion, (2) changes 
in the surface resistivity to erosion. 
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The term “rejuvenated stream” is applied to a stream system in which a renewed 
cycle of erosion begins and which may extend the drainage net after it has reached 
maturity. Rejuvenation may result from several causes, although in the Davis 





Ficure 35.—Perched or hillside stream 


sense the term is applied chiefly where it results from widespread geologic changes 
such as renewed uplift, folding, and tilting. 

Accelerated or decreased erosion may result without any such geologic changes 
if the original terrain varies in erosional resistivity or infiltration-capacity proceed- 
ing downward from the surface. Then, as erosional gradation takes place, changes 
in the critical length of overland flow x, will occur, and if these changes are abrupt 
they may result in important effects, either (1) marked increase in drainage density 
and extension and number of minor tributaries, if R; and «, decrease downward 
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from the surface, or (2) abandonment and fossilization of pre-existing streams and 
tributaries, if R; and x, increase with increased gradation. 

A third condition may also bring about changes in erosion rate and stream de. 
velopment which is more common than rejuvenation due to strictly geologic causes, 
This occurs where, as the result chiefly of climatic or cultural changes, there ig q 
change in the surface-erosional resistivity or infiltration-capacity of the terrain 
which brings about changes in the critical length «, and in the consequent develop. 
ment of drainage. 

Accelerated erosion due to the removal or replacement of an initially resistant 
surface by a less resistant surface has been appropriately described by Bailey (1935) 
as an epicycle of erosion. This term is appropriate since it implies a marked changed 
in erosional and gradational activity, superposed on the normal erosional conditions, 
Changes in erosional conditions brought about by dust storms and the formation of 
loess veneer on soil surfaces, and changes in erosional activity resulting from im- 
proper cultivation of the soil, deforestation, fires, or overgrazing of range lands, afford 
excellent examples of epicycles of erosion. 

Where a less permeable and more resistant surface layer of soil or sod overlies 
weaker or more permeable subsoil, there will be in effect two different values of x,, 
one pertaining to the surface layer, the other to the underlying material. This occurs 
where well-established grass or other vegetal cover overlies a noncohesive sandy soil 
or where there is a layer of loess or similar fine-textured material, with moderate or 
high cohesiveness, overlying more permeable and less cohesive material, such as 
sand. 

If the overlying resistive material is broken through, the value of x, pertaining to 
the underlying material governs subsequent stream development. In such cases 
the development of a drainage net is likely to be erratic and sporadic. On much of 
the area there may be but few streams. This will be true where the larger or surficial 
value of erosive resistance R; and critical distance x, are effective. At other locations 
where the smaller subsurface values of R; and x, have become effective, active and 
extensive stream development may take place. Extensive plains, for the most part 
undisturbed by erosion, may be dissected by rapidly growing and irregularly branch- 
ing systems of gullylike channels. This condition exists in the Pontotoc Ridge region 
of the Little Tallahatchie, Mississippi, drainage basin, where deep incoherent sand 
is overlain with a thin veneer of fine uniform loessal silt. In this region x, for the 
underlying sand is practically zero, and stream development may extend far above 
the x, limit for the surface material as a result of headward erosion. The author has 
observed gullies in the Pontotoc Ridge region which in some cases have extended 
not only to but somewhat beyond the topographic boundaries of their drainage basins 
(Happ et al., 1940). This has resulted from the slumping of masses of earth from 
the nearly vertical and sometimes undermined scarp formed by the erosion of the 
deep, incoherent sand. 

The destruction of vegetation by smelter fumes early in the present century in the 
vicinity of Ducktown and Copper Hill, Tennessee, brought about a new erosion 
cycle. Glenn’s early report (1911) and the author’s later observations show that 
forest and hills sometimes protected the sod locally even where the trees were killed, — 
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and where the sod was protected no erosion occurred. As described by Glenn (1911, 
p. 78): 


“The erosion starts near the bottom of a slope, and where the soil is porous rapidly cuts a steep- 
sided gully to a depth of 5 to 12 feet below the surface, where the underlying schist is as a rule still 
measurably firm. After a gully has reached its limit in depth it widens until its walls coalesce with 
the walls of adjacent gullies, by which time most of the soil has been removed.” 


Over much of the denuded area erosion has not been as complete as that above 
described. Narrow flat lands still persist between the parallel gullies, and uneroded, 
nearly flat summits of the hills are conspicuous. In some cases the gullies afford 
excellent examples of cross-grading in progress, with remnants of the antecedent rill 
surface still visible. 

The erosional topography of this region was essentially mature before denudation 
took place wherever there was a well-established sod cover. The resistivity of the 
underlying soil to erosion is, however, so small that, lacking protection, the critical 
distance x- is reduced nearly but not quite to zero. Consequently the walls between 
initial parallel ridges on steep slopes have sometimes coalesced, as described by Glenn. 
Within a few years after destruction of the vegetation the drainage density was in- 


creased locally from ten to one hundred fold, and where this occurred the end point 


of the new erosion cycle was quickly attained. 

In the gully formation in the Pontotoc Ridge region in Mississippi and in the 
vicinity of Ducktown, Tennessee, surface and subsurface resistance R; differed, the 
surface resistance being initially greater and the terrain initially stable against ero- 
sion. Reduction of surface resistance resulted from improper cultivation in the 
Pontotoc Ridge region and from partial destruction of vegetal cover by smelter fumes 
in the Ducktown region, and an active epicycle followed in each case. The forma- 
tion of arroyos on overgrazed land affords another example of an epicycle of erosion 
where the value of x; is less for underlying soil than for undisturbed surface cover. 


DRAINAGE-BASIN TOPOGRAPHY 


MARGINAL BELT OF NO EROSION; GRADATION OF DIVIDES 


In addition to controlling the drainage density and the composition of the drainage 
pattern and fixing the end point of development of a stream system on a given area, 
the critical distance x, and the belt of no erosion which it produces govern the degree 
of gradation which can occur on a given area and the extent of gradation along and 
adjacent to both exterior and interior watershed lines or divides. 

If the angle between the watershed line aa’ (Fig. 36A) and the direction of overland 
flow is A, then for a given critical length «, there will be a belt of no erosion on the 
given side of the watershed line having a width 

W, = x, sin A. 
This marginal belt of no erosion aa’cc’ is relatively permanent. It is widest, other 
things equal, where the direction of overland flow is most nearly normal to the water- 


shed line; this is usually around the headwaters of an exterior divide. The width 
of the marginal belt of no erosion decreases for a given x, as the direction of overland 
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flow becomes more nearly parallel with the direction of the divide, a condition which 
commonly occurs along lateral segments of the main divide surrounding a drainage 


basin. 





(A) Width of belt of no erosion. 
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Ficure 36.—Belts of no erosion 
(B) Longitudinal belt of no erosion. 


If aa’ (Fig. 36B) represents the exterior divide at the head of a newly exposed area, 
then, with sufficient newly exposed surface, streams will develop, starting at d and ¢. 
The entire slope from cc’ to the outlet is subject to sheet erosion. Cross-grading 
begins adjacent to the streams and spreads laterally until there remains a narrow 
belt ff’ge’ not yet cross-graded. Dashed arrows indicate directions of overland flow 
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antecedent to, and solid arrows the corresponding directions with, cross-grading. 
This belt has, however, been previously subject to sheet erosion since it lies down- 
slope from the headward belt of no erosion, and the direction of overland flow is 
parallel with the slope. The profile of the belt ff’gg’ is concave, and it lies, except 
at its ends, considerably below the original slope. The narrow belt ff’gg’ is still 
subject to cross-grading. Slight variations in surface conditions will divert most of 
the surface runoff at a given location, as at #, into one stream or the other. The 
divide between the streams will move away from the stream into which the diversion 
occurs. The direction of overland flow on the diverted area will swing around until 
it is more or less parallel with that on the adjacent cross-graded slope, and a belt of 
no erosion will develop on the side of the divide on which diversion occurs. This 
belt will have a width x, cos A, where A is the runoff angle between the diverted 
surface runoff and the antecedent slope. This angle will vary from zero to mm’n, 
and the width of the belt of no erosion on the given side will vary accordingly. 

At some other location, 7, the stream ee’ will gain the advantage in competition 
with gg’, and the watershed line will be deflected toward ff’. As a result the water- 
shed line will become sinuous, as shown by the dashed line on Figure 36B. Inter- 
mediate between # and 7 the streams will divide the runoff more or less equally. The 
watershed line will cross the center of the belt ff’gg’, but at this location most of the 
runoff will have been diverted at h, and there will be less erosion than at either / or j. 
As a consequence of the competitive development of divides the width of the belt of 
no erosion will vary from point to point, governed locally by the slope, the direction of 
overland flow, and the amount of previously undiverted surface runoff originating 
within the belt of no erosion. The watershed line will be sinuous in plan and profile, 
and the watershed ridge will be broken up into a series of irregularly spaced hills, 
often with flat crestal plateaus, and adjacent hill crests will be at about the same 
elevations. The hills will be separated by saddles, and both will be rounded not 
only as a result of the manner of their development by aqueous erosion but also by 
secondary processes, such as earth slips and rain-impact erosion. 

A favorable location for flat-top, interfluve hills is at the junction of a longitudinal 
and a cross divide. Such junctions commonly occur where there is an angle or bend 
in the parent divide. Under these conditions the flat-top hill usually has an arm 
extending out onto the interior divide. Flat-top hills and plateaus may also occur at 
intermediate locations where there is a relatively wide belt of no erosion. 

On Figure 37, aa’ and 66’ are adjacent tributary streams which developed more 
or less simultaneously on the same side of the parent stream and which flow nearly 
parallel, and crosswise of the original slope. When these streams have developed 
on an antecedent slope, cross-grading will occur, spreading laterally on both sides of 
each stream. Dashed arrows show directions of overland flow on the antecedent 
surface, and solid arrows show the corresponding directions after cross-grading by the 
streams aa’ and bb’. Most of the surface runoff on the area aa’m will be diverted 
from the parent stream into stream aa’ by cross-grading. 

Downslope from aa’ this stream can divert only the runoff from the area aa’c. 
Overland flow on the area upslope from the watershed line ac will be parallel to this 
line, while on the area acde the antecedent direction of overland flow will still persist, 
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and a belt of no erosion will develop. The stream 66’ will receive the runoff from the 
area bb’f on the upslope side and from the area bb’g on the downslope side. The 
areas aca’ and dbgb’ will have been subject to at least two cross-gradings, and asa 
result the direction of overland flow on these areas will have been turned nearly 


2 
“4 


MO}s 





3008 PUNSIS [Sra —————__ 


> 

! 

t 

oO 

3 

S 

4 

= 

3 | 
= | 
9 | 
re) | 
= | 
: : 

x 





= 
re 


_———. 
Se es cee ces eee 
_—— 
—J ee 








Ficure 37.—Belt of no erosion at a cross divide 


through a right angle. The direction of overland flow on such areas downslope from 
streams running crosswise of the original slope may, of course, have either a down- 
slope component (Fig. 37), or it may have its direction of flow reversed with respect 
to the original slope. This happens if the direction of overland flow is deflected 
through more than 90°. 

A longitudinal section along the line xx’ is shown on Figure 37A. The initial 
surface is shown by a dashed line. In spite of the fact that e is higher than f, the 
resultant slope is not materially different on the wide and narrow sides of the valley, 
fact often noticed on topographic maps. 

It has been shown that the occurrence of a belt of no erosion along an interior 
divide between streams parallel with the original slope is contingent on the develop- 
ment of components of flow across the divide by micropiracy. A belt of no erosion 





— - 2 tiie 





ya 


rior 


3100 











DRAINAGE-BASIN TOPOGRAPHY 359 


usually occurs on each side of the divide but is relatively narrow in relation to x,, and 
in the vicinity of saddles between crestal hills it may have been subject to cross- 
grading during its development. If a divide runs crosswise of the drainage basin 
(acde, Fig. 37), the belt of no erosion will temporarily be subject to longitudinal 
erosion, but presently, as a result of erosional competition, hills and saddles will 
develop, breaking up the longitudinal components of overland flow into elements 
each less than x,, as in case of a longitudinal divide. A crosswise belt of no erosion 
will usually be wider, and the interfluve hills and plateaus developed thereon will 
usually be larger and with flatter tops, than in case of a divide running parallel with 
the original slope. 

As a drainage system develops, additional belts of no erosion are introduced along 
the new interior divides, thereby reducing the portion of the total area over which 
sheet erosion can occur, other things equal. These later divides have been longer 
subject to gradation than those developed earlier, and they are generally at lower 
elevations relative to the original surface. Streams that ultimately become the 
higher-order streams of the drainage basin usually develop early in the erosion cycle, 
and their divides are usually higher relative to the original surface than those of 
lower-order tributaries. 

A belt of no erosion once developed persists throughout subsequent stages of 
gradation although subject to variations in width with subsequent cross-grading of 
the adjacent terrain. 

If the drainage basin of a tributary is narrow and steep on one side, with overland 
flow at right angles to the divide, the belt of no erosion may extend from the water- 
shed line to the stream on that side, while a flatter slope or overland flow at an acute 
angle on the opposite side may permit erosion over all or a part of the area on that 
side. 

Discussion thus far has related chiefly to the earlier stages of gradation of a drain- 
age basin where the length of overland flow is generally much greater than the critical 
length x. At later stages of stream development the critical length x, is decreased 
by successive subdivisions, with the birth of new tributaries, until finally there re- 
mains little or no intermediate length of overland flow between the belts of noerosion 
and the streams. 

A practical illustration of a belt of no erosion is afforded by slope terracing. This 
introduces a system of artificial watershed lines or cross divides on the terraced slope 
such that the remaining lengths of overland flow are everywhere less than the critical 
distance x,, and the area between a terrace divide and the next one downslope con- 
stitutes a belt of no erosion. 

Nature accomplishes a similar result in the development of a drainage-basin system 
by the successive development of tributaries of lower orders, thereby cutting down 
the length of overland flow until it does not exceed the critical length x. anywhere 
within the drainage basin. The drainage net is then complete. Since the develop- 
ment of tributaries of successively lower orders in a stream system does not go on 
indefinitely, the drainage density approaches a finite limit. Drainage density 
seldom exceeds 3.0 and commonly is 1.0 to 2.0 in humid regions where soil erosion is 
active, as shown in column 9 of Table 1. 
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INTERFLUVE HILLS AND PLATEAUS 


Within almost any drainage basin approaching maturity, especially with steeper 
slopes, relatively flat-topped interfluve hills and plateaus are scattered over the 
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FiGure 38.—Topography of an interior cross divide 
(Coosa Bald, Ga., quad., U. S. G. S.—T. V. A.) 


area. These hills and plateaus are not, as sometimes described, ungraded areas. 
They represent remnants of antecedent slopes and are therefore areas on which 
gradation was arrested when the adjacent streams and valleys developed. 

: Interior divides running crosswise of the drainage basin have broader belts of no 
erosion and are more permanent than either interior divides or lateral segments of the 
main divide running parallel with the original slope. Interior, flat-top, residual hills 
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are most common along the lines of transverse divides and at their junctions with 
interior longitudinal divides. Figure 38 shows the topography of a transverse divide 
in the drainage basin of Nottely River between Wold and Crumby creeks (Coosa 
Bald, Georgia, quad., U. S. G. S——T. V. A. map), with interfluve hills (knobs and 
balds) and saddles (gaps). 


TABLE 7.—Flat tops in upper Esopus Creek drainage basin 

















istance above . | : Average elevation, . 
“dive Bs dee —S “? tops | Apercuimate seen major as | — stream 
| 
(1) (2) (3) (4) | (5) 
a 1500 0.09 1600 600 
a | 2000 0.61 1650 610 
5.9 2500 0.11 2200 640 
6.4 2500 0.03 2350 660 
7.35 3500 0.06 3000 700 
4.30 2500 0.51 3000 700 
9.2 3500 0.05 3150 740 
9.2 2500 0.23 3150 740 
9.2 2500 0.16 3150 740 
9.7 2500 0.24 3400 760 
10.6 3000 0.15 3300 800 
10.6 2500 1.00 3300 800 
11.9 3000 1.85 3100 900 
13.4 2500 0.14 3200 | 950 
14.3 | 3500 0.05 3450 1000 
18.3 | 3000 0.04 3200 1500 











The topographic map of Esopus Creek drainage basin above Olive Bridge in south- 
eastern New York shows several residual flat-top hills and plateaus, all located along 
the lines of interior cross divides (Phoenicia, Kaaterskill, Margaretville, and Slide 
Mountain quads., U. S. G. S. topographic maps). Column 2 of Table 7 shows the 
elevations of the highest closed 500-foot contours. They are given in order of occur- 
rence proceeding upslope from Olive Bridge. Column 3 shows the area within the 
contour. Above these contours the summits are relatively flat. Column 4 gives the 
average elevation of the main divide at the same cross section. Flat tops on the same 
cross divide are usually at nearly the same elevations, but those on adjacent divides, 
even where at nearly the same distance upslope, may be of quite different elevations; 
the one located on the divide which developed later is usually lower. 

At the head of the drainage basin the marginal belt of no erosion along the main 
divide represents a portion of the original surface. While the elevations of the sum- 
mits of flat tops usually increase upslope, those on lateral divides near the head of the 
drainage basin have been graded somewhat and are consequently usually at a lower 
elevation than the divide at the head of the basin. 

Interior interfluve hills and plateaus ave, however, not always lower than the 
adjacent peripheral divide because of conditions of exposure of the original surface. 
If the original surface was warped upward or is domelike, interfluve hills may rise 
above the adjacent main divide. The occurrence of flat-top, interfluve hills and 
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plateaus results from the development of a divide between two streams under com. 
petitive conditions and requires only aqueous erosion. The locations and size 
of such hills may, however, be governed by secondary causes. 

There may be, on a newly exposed region, local areas where high infiltration. 
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Ficure 39.—Origin of ungraded or partially graded interfluve hills and plateaus 





capacity /, large initial resistance R;, or local flatness of slope, or these combined, 
make the critical distance x, abnormally large—for illustration, 100 times as great 
as for the region in general. On Figure 39, A, B, and C represent such local areas, 
the lines of overland flow in the first stage of gradation being as indicated. The 
area A receives such intense runoff at its upslope edge that it may be eroded actively 
in the first stage of gradation. Area B may be little eroded or eroded only on its 
downslope portion, while area C will be immune to erosion. In the second and sub- 
sequent stages of stream development, none of these areas may be subject to ero- 
sional gradation, but they will remain as flat-topped interfluve hills, some of them 
at the elevations of the original surface, others at somewhat lower elevations. 


CONCORDANT STREAM AND VALLEY JUNCTIONS 


A new stream develops on a pregraded slope extending away from the parent 
stream. Thus the new stream enters the parent stream concordantly. A tributary 
valley has, in general, steeper side slopes and a smaller value of x, than its parent 
stream valley. Hence a tributary valley usually grades faster than the coincident 
gradation of its parent valley, and, although younger, if its stream does not initially 
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enter the parent stream concordantly, it ultimately reaches the grade of the parent 
stream and debouches into it concordantly. If unrestricted, the tributary stream 
would cut below the level of its junction with the parent stream, but, since it cannot 
discharge below the grade of the parent stream at the junction and can easily main- 
tain its grade at the level of the latter, it continues to discharge into the parent 
stream concordantly, as stated by Playfair’s law. 








Ficure 40.—Gradation of stream valley 


STREAM-VALLEY GRADATION 


Stream and valley gradation are closely related. The stream supplies a means of © 
disposal of eroded material from the valley and fixes the minimum level of valley 
gradation. The valley tributary to the stream supplies the runoff that grades the 
stream. Stream and valley gradation proceed together, but valley gradation tends 
to lag behind stream gradation. 

In the development of the stream and valley of a first-order tributary, runoff 
volume adequate to produce erosion close to the stream will be most frequent, and 
this part of the valley will be cut down rapidly to stream level, progressing backward 
toward the divide. At any point in the drainage basin downstream from the critical 
point «-, the overland flow is increasingly charged with material in transport. This 
cuts down the transporting power of the overland flow, and near the foot of the slope 
the eroding and transporting power may be in equilibrium, and further erosion may 
thereby be inhibited. 

Figure 40 represents the cross section typical of a mature tributary stream valley 
produced by aqueous erosion. As a result of cross-grading and re-cross-grading, the 
initial surface aa’ was cut down to bfh when the stream developed. At each side is a 
belt of no erosion. For homogeneous material, valley side-slope erosion will not 
stop at the line bc. If erosion continues until the profile on the left-hand side is bd, 
the sheet flow, charged with eroded material, arriving at d must be disposed of. 
Slope is required to carry the water from d to c, and the segment dc will not be graded 
below this minimum slope. A steeper slope may be maintained from d to c because 
of sedimentation if the sheet flow from above d is overcharged with sediment with 
respect to its transporting power at the reduced slope dc. 

For turbulent flow the critical length x. varies inversely as the surface-runoff 
intensity g,. Rainstorms range from those with intensities less than infiltration- 
capacity, and which produce no surface runoff, up to the maximum intensity pos- 
sible in the given locality. In lighter storms -, will extend to the stream at c, and no 
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erosion will occur, although there may be runoff. In moderate storms -, will exten 
to some point between 0 and c, and erosion will occur only on the lower portion of 
the slope. Only in maximum storms will x, be limited to ab. with consequent erosion 
or sedimentation throughout the length bc. 

As a result of the combination of (1) decreasing frequency of higher rain intengj. 
ties, (2) the existence of a marginal belt of no erosion, and (3) inhibition of valley. 
bottom erosion by limited transporting power of overland flow, the valley cro 
section takes on an ogive or S-shaped form, with a point of contraflexure at some 
point g and a valley cross section below this point commonly nearly parabolic in form, 

A maximum storm produces erosion increasing downslope from x, and simulta- 
neously provides material in suspension proceeding downslope and thereby increases 
the width of the belt of sedimentation adjoining the stream. The combined effect 
is in general to concentrate the steepest-slope angle in the middle portion of the slope 
extending from the stream to the watershed line. 

The slope angle in the belt of steepest slope does not continue to increase indef- 
nitely as erosion continues but approaches a value corresponding to the slope angle 
for maximum erosion. Obviously if the slope became steeper, the erosion would be 
decreased instead of increased. This limiting angle of erosion may be either greater 
than or less than the angle of repose of the soil material when wet. If it is less than 
the angle of repose, then the slope will ultimately become stable at about the angle 
of maximum erosion. If the angle of slope for maximum erosion is greater than the 
angle of repose, the slope will become subject to earth slips as well as to erosion, and 
each earth slip exposes new material, thus favoring continued erosion. Earth slips 
do not occur on all mature slopes—in fact, their occurrence is relatively uncommon, 
and this fact is explained in the manner above described but cannot readily be ex- 
plained if it is assumed, as has frequently been done heretofore, that erosion continues 
to increase indefinitely with increase of slope. Land slips, if they occur, generally 
tend to prolong the process of valley-slope gradation and often decrease the width 
of the belt of no erosion. 

When a first-order tributary develops, the length of overland flow along the course 
of the rill channels on the antecedent surface is generally less than x, over most 
of the area tributary to the stream. Hence only limited subsequent cross-grading 
can occur. Such cross-grading as does occur is often local and feeble and is usually 
confined mostly to belts closely adjoining the stream on either side. This often 
results in a relatively deeply incised V-shaped stream channel adjoined on either side 
by uneroded upland. In this case the stream and its valley are not truly concordant 
in the upper reaches of the stream although they become concordant at the mouth 
of the stream, in accordance with Playfair’s law. 

Valley gradation, while it may lag behind stream development, is usually close to 
completion when the latter is completed. As a result the drainage basins of Ist 
order streams are likely to be born mature or nearly so and yet have the appearance 
of youth in the physiographic sense, in that their slopes are steep and their immediate 
stream valleys narrow and V-shaped in cross section. They remain perpetually 
youthful in appearance not because there is no remaining gradation which could be 
accomplished but because there are no adequate tools available by which it can be 
accomplished. 
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Transverse and longitudinal valley profiles both result from the same hydrophys- 
ical processes operating on the same terrain, and they are closely related. The re- 
sultant slope is the result of cross-grading at a given point. The transverse slope is 
the component of resultant slope at a right angle to the stream. Longitudinal slope 
is the component of the resultant slope parallel with the stream. 

Many attempts have been made, beginning with Unwin (1898), to derive an equa- 
tion for a stream profile. The stream profile at bankful stage is essentially the same 
as the longitudinal profile of the valley bottom. None of the existing equations is 
wholly satisfactory. A rational equation for stream profiles must take into account 
not only the laws of hydraulics and channel dynamics but those also of sheet erosion 


and storm frequency. 
TYPICAL OVOID FORMS OF DRAINAGE-BASINS 


Whether a drainage basin develops from the head downstream or from the mouth 
upstream depends on the manner in which the land becomes exposed. On lands 
exposed by vertical uplift of oceanic marginal lands, the streams would obviously 
develop from the headwaters downstream, while on lands exposed by retreat of glacial 
ice, the drainage basin and stream system may develop from the mouth upstream. ~ 
In either case the points of first emergence ultimately become either the heads or 
mouths of principal streams because they provide the greatest lengths of overland 
flow and hence the points at which /, first becomes greater than x,. There is an angle 
or bend at the point of first emergence which remains as a landward-directed bend 
at the head of the main divide. Hence the head of a main drainage basin is usually 
somewhat rounded. . 

As a result of change of direction of overland flow by cross-grading (Fig. 30), the 
area acb between two parallel streams never becomes tributary to these streams as 
long as the coast line remains at ab. Areas tributary to rilled surfaces which initially 
drained part of the area aa’bb’ directly to ab will be captured by competitive cross- 
grading, and in this way the drainage basin of the streams aa’ and 6d’ will expand 
laterally until a permanent competitive divide cc’ is established. 

Similar competitive divides will be established to the right of stream 5d’ and to the 
left of stream aa’. An area of direct drainage similar to abc always remains between 
each two streams, with the result that the final forms of drainage basins are usually 
pear-shaped or ovoid. 

Elsewhere (Horton, 1941) the author has shown that the average of the forms of 
drainage basins of several of the larger rivers of the world is a slightly asymmetrical, 
pear-shaped figure with the apex at the outlet end. Also, assuming valley cross 
sections to be approximately parabolic, it is shown (Horton, 1941) that the trace of 
the intersection of a parabaloid with its apex upstream and a plane inclined to the 
axis of the parabaloid is a similar ovoid figure. Hence the hydrophysical laws and 
processes, as a consequence of which valley sections are generally similar to, though 
not actually parabolas, lead directly to the development of drainage basins of typical 
form. 

In spite of modifications of individual drainage basins by geologic structures, the 
semblance to the typical ovoid form is still preserved in most drainage basins. 
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DEVELOPMENT OF LARGE DRAINAGE BASINS 


In this paper, lands exposed by changes of a coast line have been used for purposes 
of illustration. If an entire continent was exposed progressively the principles 
outlined would still apply and would lead to the formation of stream systems de. 
pendent on the manner of exposure and the position of the continental divide. [y 
North America, for example, upwarping of the northern portion of the continent and 
a gentle tilt of the remainder to the southeast would provide requisite conditions for 
the existing stream systems. 

Waters of the earth are probably chiefly of volcanic origin. If at an early stage 
in geologic history a continent was exposed at a time when rainfall was light and 
oceans shallow or wanting, and there followed a period of gradually increasing rain- 
fall, then a stream system would develop with a tendency to greater concentration 
in large drainage basins than occurs where exposure of the area is gradual. 

If an entire continent was exposed before erosion begins, the maximum length of 
overland flow for the whole continent would occur at the start. Smaller coastal 
basins would be more limited in their extension inland than with a gradual exposure 
of a coastal margin. For-such a continental exposure a few major drainage basins 
would generally develop, draining the greater part of the area. The development of 
each drainage basin would take place in much the same manner as one of the larger 
basins on a progressively exposed coastal margin. Smaller coastal streams would be 
restricted in tributary area by competition with the stronger major streams in much 
the same way in the two cases. 

There would be certain differences in the operation of the processes of erosional 
gradation for the continental area: 

(1) Even an approximation to homogeneity of the area with respect to the control- 
ling factors: rain intensity, infiltration-capacity and surface resistivity to erosion, 
would not be likely to occur. 

(2) A single maximum storm covering the whole of the continent is improbable. A 
major continental drainage basin would apparently develop in sections, and the 
sections would combine one after another when erosion attained requisite stages. 
There would be opportunities for large-scale stream piracy, such as appears to have 
occurred in the upper part of the area originally tributary to Colorado River. 

The development of great drainage basins, such as those of the Mississippi, Amazon, 
and Nile rivers, must be treated on an individual basis. The recognized persistence 
of major stream systems, barring such catastrophic events as lava overflows or glacial 
submergence, suggests that a re-examination of the early erosional history of the 
earth in the light of hydrophysical processes is well worth while. 


DAVIS STREAM-EROSION CYCLE 


In discussions of the Davis erosion cycle the cycle is in effect assumed to begin 
after the development of at least a partial stream system. The initial development 
of streams is considered to be either fortuitous or governed by antecedent geologic 
structures and is largely taken for granted (Davis, 1909). The author has consid- 
ered stream development and drainage-basin topography wholly from the viewpoint 
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of the operation of hydrophysical processes. In the Davis theory the same subject 
is treated largely with reference to the effects of antecedent geologic conditions and 
subsequent geologic changes. The two views bear much the same relation as two 
pictures of the same object taken in different lights—the results are not necessarily 
in conflict; each supplements the other. The hydrophysical concept appears to be 
more fundamental because it carries back to the original, newly exposed surface. 
In comparing the two viewpoints much depends on the meanings given to the words 
“youth,” ‘‘youthful,” “maturity,” and ‘“‘peneplain.” 

As commonly applied in connection with the Davis erosion cycle, “young” and 
“youthful” relate to an area where there has been but little erosional gradation. 
“Mature” means that all or nearly all the gradation which can result from the opera- 
tion of existing agencies has been accomplished. On an ares where infiltration- 
capacity and surface resistance are sufficiently high, little erosional gradation may 
be possible under existing conditions. Such an area is actually mature although 
from its surface appearance it would be classified as youthful. 

Incipient drainage and low drainage density are accepted as prima facie evidence 
of youth. Extensive dissection and high drainage density on a given area are ac- 
cepted as necessary and sufficient proof of maturity. These are usually sufficient, - 
but they are not necessary as conditions precedent to maturity. An area of low 
drainage density and with little dissection may have been born mature if the original 
infiltration-capacity, initial resistance, and rainfall appurtenant to the drainage area 
were such that the length of overland flow was not sufficient to induce erosion on any 
part of the area. Such an area represents erosionally mature topography. It is 
mature not because much has been accomplished but because nothing more can be 
accomplished without rejuvenation by geologic agencies or as a result of climatic 
and cultural changes. 

Conversely, but less frequently, an area occurs with moderate dissection and a fairly 
high drainage density which is in a youthful stage because of a small value of 2. 
Such areas are most likely to be found where an epicycle of erosion is in progress. 
That they are not common elsewhere indicates that gradational development of 
drainage basins is far more generally complete with respect to pre-existing conditions 
than is usually supposed or assumed. The error has resulted largely from a careless 
and, as now appears, unjustifiable application of the Davis concepts rather than from 
error in the concepts themselves if correctly interpreted. For example, areas on the 
Highwood, Illinois; Rochester, New York; Dunlop, Illinois; Oak Orchard, New 
York; Parmelee, North Carolina; Anson, Texas; Fargo, North Dakota-Minnesota; 
and Oberlin, Ohio, topographic sheets of the U. S. Geological Survey have all been 
referred to as representing “‘stream erosion—youth” (Salisbury and Atwood, 1908). 
In spite of the low drainage density and lack of dissection of most of the areas re- 
ferred to, the stream systems are all or nearly all relatively mature because of high 
infiltration-capacity, flat slopes, or high surface resistance to erosion or a combination 
of these factors. Erosion on many of these areas, if it occurs, must be ascribed to 
conditions which have induced either a local or general epicycle of erosion over part 
or all of the area. 

There is often need to refer to actual age. Tributary streams and their basins 
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are usually younger, in this sense, than higher-order streams or basins. The ters 
“little dissected” and “well dissected” seem preferable to “youthful” and “mature” 
in describing erosional status. There may then be areas which are (1) young, littk 
dissected; (2) young, well dissected; (3) old, little dissected; (4) old, fully dissected, 

The term “‘peneplain” seems inappropriate. The ultimate surface of erosion 
within a main basin boundary is neither “almost a plane,” as the prefix “pene” 
implies, nor is it usually as close to being a plane as was the original surface are 
from which it has been derived. It seems better to call it a “base surface” in com. 
parison with the original surface below which the stream basin is developed. The 
base surface at its downslope end is at “‘base level’ in the usual sense. The base 
surface is, however, generally concave upward except along divides, and its margins 
intersect the initial surface around the upstream portion of the drainage area. The 
ultimate base surface is, under ideal conditions, closely similar to a segment of a 
parabaloid cut by a plane which is not parallel with the axis of the parabaloid. The 
parabaloidal surface is ribbed with ridges which represent the divides between 
streams. 

Wooldridge and Morgan (1937) use “the invaluable concept of the cycle of erosion 
as initiated by W. M. Davis.” They state (p. 184): 


“Some writers have argued that the cycle of erosion can never have run its full course and that the 
peneplain is an unrealized and unrealizable abstraction.” 


They invoke geologic factors, as others have done, to complete the erosion cycle— 
in other words, the Davis erosion cycle is not completed by erosion per se. The hy- 
drophysical concept neither denies the effects nor invokes the operation of uplift. 
It carries the matter of basin development only so far as it can be carried by purely 
erosional processes, and it shows that there is a definite end point to the development 
of streams and valleys by aqueous erosion and leaves room for the survival of hills 
and plateaus between valleys and which are not subject to further erosion or pene- 
planation. It provides a better foundation than has heretofore existed for the inter- 
pretation of the effects of changes of geologic conditions in relation both to the subse- 
quent march of the erosion cycle and in relation to changes of drainage patterns and 
drainage composition occasioned thereby. 

In accordance with the hydrophysical concept most of the observed gradation of 
divides takes place before the streams which are separated by the given divide are 
developed—in other words, the terrain where the divide is located is graded in ad- 
vance at a time when sheet erosion is taking place along or across the line which sub- 
sequently becomes the divide. Interfluve hills and plateaus are remnants of this 
pregraded surface, and when once formed they are permanent features of the topogra- 
phy. The whole concept of ultimate development of a peneplain appears to be 
founded on the idea that grading of interior divides continues indefinitely and is 
accomplished by the streams they separate, whereas in the case of the hydrophysical 
concept the gradation is already, for the most part, accomplished when the adjacent 
streams originate. 

On drainage basins near maturity, erosion occurs only sporadically and in local 
patches where local conditions provide requisite length of overland flow, runoff 
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yolume, slope, and low local resistivity to erosion, or where an epicycle of erosion 
occurs. 

As regards advanced stages of erosional development of a drainage basin the prin- 
cipal differences between the Davis concept and the hydrophysical concept seem to 
be that the former does not provide any definite end point of erosional development, 
whereas the hydrophysical concept does provide a definite end point for both stream 
and valley development. 

The physical and mathematical treatment of the subject establishes rational 
quantitative relationships between the interpretation of observed phenomena ac- 
curately and with confidence in the correctness of the results. For example, it has 
long been known that there is a relation between surface erosion and slope, but 
without knowing the physical basis of this relationship gross errors could easily be 
made in interpreting valley-slope erosion. The equations may be applied to the 
study of individual cases. Such possible applications are of infinite variety and must 
be left to others. It has been shown that some of the equations are of practical use 
in engineering problems, as in determining flood-crest modulation by channel storage 
and the requisite spacing of soil conservation terraces. 

In conclusion, what has been given is a framework or outline of drainage-basin - 
development along hydrophysical lines rather than the completed picture. It is 
hoped that the reader will find that a new, more definite, and more quantitative 
meaning has been imparted to Playfair’s law and the Davis concept of the stream 
erosion cycle. It is also hoped that the reader will find stimulation to further study 
and research. 
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